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♦The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  =  1  event/s. 

♦♦The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 

A  more  complete  listing  of  conversions  may  be  found  in  "Metric  Practice  Guide  E  380-74, " 
American  Society  for  Testing  and  Materials. 
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SECTION  1 
INTRODUCTION 

The  ROSCOE-IR  computer  program  provides  a  means  for  calcu¬ 
lating  IR  observables  from  nuclear  detonations  (NUDETS).  The  general 
computational  scheme  in  ROSCOE  is  to  calculate  the  radiance  along  a 
line-of-sight  (LOS)  from  the  arbitrary  location  of  the  sensor  looking 
toward  the  NUDET  source  region.  By  projecting  a  number  of  lines-of- 
sight  to  the  source  and  its  vicinity,  it  is  possible  to  construct  an 
image  for  imaging  type  sensors.  For  non-imaging  sensors  with  finite 
f ield-of-view  (FOV),  it  is  possible  to  assess  the  radiance  transmitted 
directly  from  the  source  and  separately  to  assess  the  contribution  at 
the  detector  of  backgrounds  and  scattering  contributions  along  the 
lines-of-sight  not  leading  directly  to  the  source.  In  each  case,  the 
NUDET  source  is  modeled  as  a  function  of  time  and  wavelength  and  the 
radiance  along  the  LOS  is  calculated  by  integrating  along  the  LOS  toward 
the  source.  The  transmission  of  the  medium  and  the  contribution  to 
radiance  by  differential  elements  of  the  LOS  through  emission  and  scat¬ 
tering  into  the  LOS  are  taken  into  account.  Contributors  to  the 
radiance  are  thus  NUDET  source  radiation  at  the  end  of  the  LOS;  emission 
of  the  atmospheric  elements  between  source  and  sensor;  and  scattering 
into  the  LOS  of  solar  radiation,  earth-emitted  radiation,  and  NUDET 
radiation  from  source  points  off  the  LOS. 

The  presence  of  natural  clouds  in  the  atmosphere  influences 

the  observable  IR  radiance  so  strongly  that  the  inclusion  of  a  cloud 

model  in  the  ROSCOE-IR  program  is  essential  for  the  realistic  assessment 

of  IR  observables.  In  most  areas  of  the  world  of  interest  in  connection 

with  nuclear  explosions,  clouds  are  present  about  two  thirds  of  the  time 

and,  in  fact,  clouds  will  lie  in  a  LOS  about  one  half  the  time.  Clouds 

produce  several  kinds  of  effects  which  are  important  both  to  imaging 

and  non-imaging  sensor  systems.  The  most  obvious  effect  of  clouds  is 

to  obscure  the  source  or  to  so  reduce  its  intensity  that  it  is  not 

detectable  by  the  sensor.  Thick  clouds  provide  direct  attenuation 

4  5 

factors  as  high  as  10  or  10  even  at  the  most  transparent  wavelengths. 


In  wavelength  bands  where  there  is  strong  absorption,  such  as  at  2.7 
microns,  attenuation  factors  of  this  sort  occur  even  for  rather  thin 
clouds.  This  attenuation  by  clouds  has  led  to  the  well  known  rule  of 
thumb  that  at  many  wavelengths  NUDETS  are  not  visible  from  above  until 
the  nuclear  cloud  has  risen  through  the  cloud  tops. 

A  second  effect  of  clouds  even  at  moderately  transparent 
wavelengths  is  to  destroy  the  image  of  the  source  by  multiple  scattering 
so  that  the  apparent  source,  when  viewed  from  above,  is  a  circular 
diffuse  spot  on  the  cloud  top  which  may  be  many  kilometers  in  radius. 

As  the  nuclear  cloud  rises,  this  diffuse  source  becomes  smaller  in 
radius,  turning  gradually  into  the  nuclear  cloud  itself  as  it  emerges 
from  the  cloud  top.  Many  seconds  typically  are  involved  in  the  rise  of 
the  nuclear  cloud  from  low  altitudes  to  cloud-top  altitudes. 

A  third  effect  of  clouds  arises  in  scattered  or  broken  clouds 
so  arranged  that  one  cloud  obscures  the  direct  LOS  to  the  source,  while 
other  LOS's  lead  to  cloud  edges  which  are  illuminated  by  the  source. 

Since  cloud-edge  reflections  may  contribute  greater  radiance  than  that 
produced  along  an  obscured  LOS,  this  effect  gives  rise  to  a  number  of 
lighted  clouds  in  the  vicinity  of  the  source  brighter  than  the  (obscured) 
source.  Other  effects  of  clouds  include  a  higher  background  radiance 
than  in  the  absence  of  clouds  because  of  the  high  albedo  of  cloud  tops 
for  solar  radiation  and  a  spottiness  in  the  background  when  the  clouds 
are  scattered  or  broken. 

It  has  been  the  purpose  of  the  ROSCOE  natural  cloud  project 
to  provide  a  cloud  model  and  the  capability  in  ROSCOE-IR  to  calculate 
cloud  effects  on  systems  which  look  at  NUDET  environments.  Satellite 
systems,  particularly,  have  been  kept  in  mind  during  the  development  of 
this  natural  cloud  capability  and  the  spectral  dependence  of  effects 
has  been  sufficiently  treated  in  the  2-  to  5-micron  region  in  order  that 
realistic  results  will  be  provided  for  sensors  whose  spectral  sensitivity 
lies  within  this  range. 

The  complexity  of  clouds  and  the  variability  of  cloudiness  in 
a  given  region  poses  an  immediate  problem  in  cloud  modeling.  It  is 
typical  that  any  choice  of  cloud  configurations  used  for  transmission 
calculations  leads  to  criticism  by  some  reviewer.  Such  questions  arise 
as,  why  does  the  cloud  obscure  the  source  when  it  is  quite  likely  that 
the  source  will  be  visible  through  breaks?  What  would  be  the  effect  if 
the  source  were  half  obscured  and  half  visible?  Why  worry  about  the 
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brightness  of  illuminated  clouds  when  the  direct  view  of  the  source  will 
be  so  much  brighter?  Why  worry  about  obscuration  of  the  fireball  when 
waiting  awhile  will  allow  it  to  break  through  the  clouds?  Each  user  who 
contemplates  cloud  effects  tends  to  think  of  configurations  which  lead 
to  maximum  or  minimum  impact  on  a  system  and  tends  to  focus  on  worst 
cases  or  best  cases.  Each  possible  cloud  configuration  has  some  validity 
in  the  problem,  but  a  measure  of  the  likelihood  of  the  case  is  important 
for  overall  assessment. 

For  the  ROSCOE-IR  program  a  cloud-effects  capability  has  been 
provided  which  allows  the  user  to  specify,  within  limits,  the  cloud 
scene  of  his  choice  so  that  the  effects  of  scenes  which  he  considers 
likely  to  be  important  can  be  assessed.  This  capability  is  called  the 
deterministic  model. 

In  addition  to  the  deterministic  model,  a  statistical  cloud 
model,  based  on  historical  observation,  is  provided  for  ROSCOE.  Fre¬ 
quency  of  occurrence  of  cloud  types  and  average  descriptions  of  the 
properties  of  each  of  the  cloud  types  are  obtained  from  a  global  cloud 
study  and  data  bank  prepared  for  NASA  (Reference  1).  By  the  use  of  this 
statistical  natural  cloud  model,  the  ROSCOE  program  calculates  the 
probability  distribution  for  radiance  along  a  LOS  and  thus  determines 
the  extremes  in  radiance  which  may  be  exhibited  in  the  presence  or 
absence  of  clouds.  For  a  problem  which  has  been  done  by  the  deter¬ 
ministic  mode,  the  likelihood  of  such  a  result  based  on  natural  cloud 
occurrences  can  be  judged. 

The  physical  and  geometric  treatment  of  these  two  cloud 
models  is  described  in  the  following  section. 
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SECTION  2 
CLOUD  MODELS 

The  ROSCOE-IP.  program  performs  very  detailed  calculations  to 
determine  radiance  along  the  line-of-sight  to  a  sensor.  The  exact 
geometric  location  of  the  sensor  and  sources  including  the  sun,  the 
earth's  surface,  and  the  composition  of  the  atmosphere  as  a  function  of 
location  are  taken  into  account.  The  radiance  may  be  calculated  at  any 
wavelength  selected  in  the  2-  to  5-micron  band.  Source  spectral  powers, 
absorption  and  scattering  cross-sections,  reflectivities,  and  other 
physical  properties  affecting  radiation  transport  must  be  known  in 
detail  as  a  function  of  wavelength.  A  cloud  model  for  use  in  conjunction 
with  ROSCOE-IR  must  thus  have  the  same  detailed  physical  and  geometric 
description  to  enable  these  calculations  to  be  performed  with  clouds  in 
the  scene.  The  deterministic  and  statistical  cloud  models  have  been  so 
constructed  that  each  calculation  is  made  in  the  presence  of  clouds 
sufficiently  described  in  physical  properties  and  geometry  that  such 
detailed  transmission  calculations  can  be  performed.  For  the  statistical 
clouds  a  whole  set  of  such  calculations  is  done  within  the  cloud  model 
and  selected  fractiles  of  the  integral  distributions  of  results  so 
obtained  are  provided  as  output. 

2-1  DETERMINISTIC  CLOUD  MODEL 

The  introduction  of  deterministic  clouds  into  the  NUDET 
environment  involves  a  specification  sufficiently  complete  to  allow 
radiation  transmission  calculations.  Various  deterministic  cloud 
shapes  have  been  used  in  other  radiation  transmission  programs.  These 
include  modeling  of  clouds  by  horizontal  slabs;  spheres;  hemispheres; 
round,  rectangular,  and  hexagonal  prisms;  etc.  For  the  ROSCOE  natural 
cloud  routines  we  have  chosen  ellipsoids.  By  the  specification  of  the 
three  axes,  these  may  vary  from  spherical  shapes  to  elliptical  pancake 
shapes  to  very  elongated  rod-like  shapes.  The  three  axes  are  required 
to  lie  in  the  vertical  direction  and  in  the  horizontal  plane  and  the 
orientation  of  the  longer  axis  in  the  horizontal  plane  must  be  specified. 
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Location  in  three  dimensions  of  the  center  of  the  ellipsoid  then  comple¬ 
tes  the  specification  of  the  geometric  cloud  boundary.  Any  number  of 
clouds  up  to  20  are  permitted  in  the  scene  so  long  as  they  are  not 
intersecting.  This  deterministic  placement  allows  the  investigation  of 
the  effects  of  obscuration  of  the  source  by  clouds,  the  brightness  and 
extent  of  illuminated  cloud  edges,  and  the  brightness  of  edges  surround¬ 
ing  a  hole  through  which  the  source  is  viewed.  The  cloud  type  which  is 
being  represented  is  also  identified  as  one  of  the  naturally  occurring 
types  so  that  the  water  droplet  density  and  size  distribution  may  be 
specified.  A  sample  obscuration  case  is  described  in  Section  4. 

Transmission  calculations  for  both  deterministic  and  statistical 
cloud  models  have  been  done  by  diffusion.  The  absorption  and  scattering 
cross-sections  for  the  water  droplets  and  the  scattering  phase  functions 
have  been  calculated  (Reference  2)  by  the  Mie  theory  of  scattering  of 
photons  by  spherical  droplet  distributions.  The  droplet  size  parameters 
from  Reference  1  (Table  6-1)  were  used  and  are  given  in  Table  2-1.  These 
mean  radius  and  variance  values  were  used  to  generate  a  gamma  distribution 
in  droplet  size  for  each  cloud  type. 

Such  calculations  were  made  for  droplets  in  the  range  of  radii 
for  4.5  to  10  pm  for  which  Mie  phase-function  parameters  were  readily 
calculable.  These  sizes  are  characteristic  of  droplet  sizes  for  cloud- 
types  1  through  10  in  Table  2-1.  The  Mie-theory  parameters  for  40-pm 
radius  droplets  are  extremely  difficult  to  calculate  and  were  regarded 
as  being  beyond  the  scope  of  the  investigation. 

The  so-called  40-gm  droplets  are  in  the  usual  situation  actually 
ice  crystals.  The  approximation  of  simulating  such  large  crystals  by 
spherical  droplets  would  be  dubious,  since  the  dimensions  are  large 
enough  that  irregularities  are  several  wavelengths  in  size. 

The  application  of  the  usual  scattering  and/or  diffusion  con¬ 
cepts  would  be  seriously  complicated  by  the  fact  that  such  crystals  are 
typically  spatially  oriented,  i.e.,  their  axes  tend  to  be  aligned.  For 
example,  in  diffusion  this  would  imply  an  anisotropic  medium,  not  just 
anisotropic  scattering. 

For  these  reasons  Mie-theory  parameters  were  not  obtained  for 
cloud-types  11  through  14.  Instead,  we  idealized  the  droplets  as  being 
40-pm  radius  water  droplets  and  employed  small-wavelength  approximations 
as  described  at  the  end  of  Section  3-2.16. 


The  phase  functions  computed  from  Mie  theory  have  been  used 
to  determine  the  average  cosine  of  the  scattering  angle,  TT.  From  the 
cross  sections  and  TT  were  constructed  the  macroscopic  transport  param¬ 
eters  which  are  used  in  the  diffusion  calculations.  Ten  wavelength 
values  in  the  2-  to  5-micron  range  have  been  selected  in  such  a  way  that 
the  principal  variations  in  cross  section  are  retained  when  linear 
interpolation  between  these  wavelength  values  is  made.  The  chosen  wave¬ 
lengths  are  2.0,  2.5,  2.55,  2.65,  2.75,  2.85,  3.2,  3.5,  4.0  and  5.0.  For 
wavelengths  between  these  values,  linear  interpolation  is  used  on  the 
macroscopic  total  cross-section,  E,  the  macroscopic  scattering  cross- 
section,  Es,  and  the  average  cosine  of  the  scattering  angle,  y.  Since 
the  absorption  cross-section,  E  is  the  difference  between  the  total 
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cross-section  and  the  scattering  cross-section,  the  necessary  quantities 
are  available  to  determine  the  macroscopic  transport  cross-section,  Etr, 
the  diffusion  coefficient,  D,  and  the  relaxation  rate,  K,  which  are  the 
constants  required  for  solution  of  the  diffusion  equation.  According 
to  the  "standard"  evaluation  of  diffusion  parameters  in  terms  of  macro¬ 
scopic  cross-sections,  we  have  the  relations 

-1  2s  /  ^  a\  £  -1 

Etr  =  E  -  TIEs(km  1 ) ,  D  =  gg°  ■  (km)  and  K  =  (-^1  (km  . 

However,  that  evaluation  is  strictly  applicable  only  in  the  limit  of 
small  absorption.  A  more  generally  applicable  evaluation,  which  we 
adopted  for  numerical  results  in  this  report,  is  described  in  Section 
A-1.2. 

For  the  deterministic  case  the  ellipsoidal  boundary  leads  to 
great  mathematical  complexity  in  applying  boundary  conditions.  The 
problem  therefore  has  been  treated  by  means  of  a  Green's  function  and 
application  of  Green's  theorem  so  that  the  problem  is  reduced  to  a  two- 
dimensional  surface  problem.  The  surface  integrals  can  conveniently 
be  treated  numerically  by  partition  of  the  surface  into  a  number  of 
facets  and  summing.  The  general  problem  is  treated  in  Section  A-3. 

A  typical  geometry  for  the  deterministic  cloud  problem  is  shown 
in  Figure  2-1.  A  point  source  S  illuminates  an  ellipsoidal  cloud  and 
gives  rise  to  radiance  in  the  direction  of  a  detector  D.  The  cloud 
surface  is  naturally  divided  into  four  regions  by  the  terminators 
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corresponding  to  the  source  and  detector  points.  Thus,  there  is  an 
illuminated  and  a  shadowed  portion  of  the  cloud  and  a  visible  and  an 
invisible  portion  of  the  cloud.  These  combine  to  give  four  categories 
of  cloud  surface,  only  two  of  which  are  visible  to  the  detector.  Three 
lines-of-sight  are  shown  in  Figure  2-1,  the  lowest  of  which  intersects 
an  illuminated  visible  region  while  the  two  upper  ones  intersect  a 
visible  but  non-illuminated  region.  Radiance  from  this  latter  region 
arises  entirely  from  emergence  of  photons  which  have  diffused  through 
the  cloud  from  the  lower  illuminated  region.  Although  diffused  photons 
contribute  to  the  lowest  line-of-sight  also,  it  is  in  fact  dominated  by 
photons  reflected  from  the  surface.  This  sector  of  the  cloud  is  there¬ 
fore  treated  by  bidirectional  reflectance,  taking  into  account  the 
incident  and  emergent  directions  and  the  local  normal  to  the  cloud 
surface.  Bidirectional  reflectance  coefficients  for  this  purpose  were 
calculated  by  Monte  Carlo  for  a  thick,  flat  cloud  of  stratocumulus  (Sc) 
droplet  density  and  size  distribution.  (This  set  of  bidirectional 
reflectance  coefficients  is  used  as  an  approximation  to  the  set  for 
each  of  the  other  13  cloud  types.)  The  emergent  photons  at  45  directions 
in  exit  zenith  angle  and  azimuth  were  calculated  for  four  incident 
zenith  angles.  These  tabular  data  of  bidirectional  reflectance  coef¬ 
ficients  are  interpolated  for  use  in  the  deterministic  cloud  model. 

The  angular  distribution  for  the  photons  emerging  from  diffusion  trans¬ 
mission  to  the  non-illuminated  surfaces  is  taken  to  be  diffuse  (or 
Lambertian)  in  angular  distribution. 

For  evaluation  of  the  photons  into  and  out  of  the  deterministic 
cloud,  the  cloud  is  divided  into  facets  as  shown  in  Figure  2-2.  The 
facets  are  constructed  so  as  to  lie  between  lines  of  constant  0  and 
constant  4> ,  with  more  facets  of  smaller  area  being  concentrated  near 
the  regions  of  maximum  curvature  in  the  ellipsoid.  In  the  cloud 
coordinate  system,  if  a  is  the  semi-major  axis  lying  in  the  horizontal 
plane  and  b  is  the  semi-minor  axis,  then  the  division  of  the  azimuthal 
angle  <)> ,  where  <p  is  measured  from  the  direction  of  the  long  axis,  is 
according  to  the  prescription 

^6-m  =  “§  ~  ^m  ’  A<t>  =  5 

/b\(n-2)S 

^n  =  %-l  +  -  *1  =  0;  =  2  1  2  <n  <N  =  5 
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ellipsoidal  cloud 


Given  8,  this  can  be  solved  for  A(j>,  the  angular  dimension  of  the  first 
facet,  and  thus  for  each  of  the  <f>  in  a  quadrant.  Similarly  the  polar 
angle  9,  measured  from  the  zenith,  is  divided  in  a  like  manner  according 
to  the  formula 


e  =  e  ,  + 

n  n-1 


«© 


(n-2)ot 


Bj_  =  0; 


2  <  n  <  N  =  5 


The  explicit  formulas  for  9^  and  cf> ^  are  given  in  Section  3  with  the 
description  of  subroutine  CLDGEOM.  Here  a  and  9  are  arbitrary  exponents 
chosen  to  accentuate  the  grouping  of  facets  at  regions  of  high  curvature. 
A  value  of  a  =  3  =  0.8  has  been  used.  These  facets  are  then  reflected 
through  the  principal  planes  into  the  other  seven  octants  of  the 
ellipsoid. 


where 


The  area  of  facet  i,j  is  given  by  the  expression 

f%+ 1  /’9i+l  2  2  2  * 

Aij  =  /  d<t>  /  d9  r  sine  (1  +  S  +  T  )5 


2 

cos  9 


r“2  =  sin2e(^4  +  +  £2^ 

\  a  o'  c 

Tr2  2  r2  2 
S  =  sine  cos0  — g  cos  <p  +  —g  sin  <p 

La  b 

r  r2  r2l 

T  =  sin0  sin<f>  cosij)|-  — ^  +  — g 

L  a  b  J 


Tr2  2  r2  2 
S  =  sine  cosQ  — ^  cos  <p  +  — g  sin  <j 

La  b 

T  .2  2-i 


This  integration  is  performed  numerically. 

As  discussed  in  Section  A-3,  the  solution  for  the  ellipsoid 

illuminated  by  a  point  source  is  in  the  form  of  the  outgoing  photon  flux 

2 

(radiant  exitance),  i.e.,  watts/kilometer  ,  per  watt  emitted  by  an 
isotropic  point  source.  Illumination  by  the  source  is  converted  to 
ingoing  photon  flux  (irradiance)  at  each  of  the  illuminated  facets  by 
computing 

S  cos  9 . 

J_  =  - — L 

4ttR 
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for  a  unit  source.  Specifically,  if  J+(r)  and  J_(r)  are  the  outward 
and  inward  fluxes,  respectively,  then 


J+( r ' )  =  S(r 


S(r '  ) 


)  -  f  J+(r)G+(r,r '  )  do(r)  | 
-J_(r')  -  J"  J_( r )G_( r , r ' )  da(r)  I 


(55)*&(59)* 


G+(r,r ' ) 


A  ■+ 

n(r) 

2tt 


grad  G(r,r’  )  +  G^r  >r '  ) 


(55) 


G(r , r ' )  =  e_KR/R  (53) 

grad  G(r.r')  =  (VG)£  =  G(r.r')  (K+R-1) 

R  =  | 


where  a  and  b  are  given  in  Table  A-2.  Note  that  the  gradient  of  G(r,r') 
is  evaluated  at  the  inward-flux  position  r.  In  the  application  of 
this,  we  used,  in  our  early  work,  the  approximate  solution,  J2  +(r'): 


Jj  +(r ' )=S(r ' )/ 


(r , r ' )  da( r ) 


(60a)*+ 


J2 ,  +( r '  )  =  'J1>+(r'  ) 


J  (Jj  +(?)-J1  +(r ' ) )G+(r,r' )do(r) 


1+ 


( 60b )  *  +  + 


r '  )do(r ) 


In  the  computation,  the  integrals  are  approximated  by  sums, 

C(r')  =  f  f(r)G+(r,r’ )do  =  ^  f  (fj  )G+(r  j  ,  r^A.. 

+  Vt  (-  ±  l) 

*Equation  numbers  in  Appendix  A. 

'The  form  of  Ji  +  is  obtained  by  ignoring  the  last  term  on  the  right- 
,,hand  side  of  Equation  (59)  in  Appendix  A. 

' 'j  is  obtained  by  combining  Equations  (59)  and  (60a). 

^  J  + 
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where  is  the  area  associated  with  the  point  r ^ is  the  average 
radius  of  curvature  at  that  point,  and  f(r)  represents  the  multiplier 
of  G+  in  the  integrals. 

Eventually  we  found  an  example  of  J0  being  negative. 

Accordingly,  after  considerable  review  of  the  situation,  we  have  settled 

for  using  J-^  +  in  the  computer  code  until  a  better  solution  is  found. 

The  geometry  of  the  computation  is  summarized  in  Figure  2-3. 

A  line-of-sight  (LOS)  of  interest  from  detector  to  cloud  surface  has 

been  established.  The  quantity  of  interest  is  the  cloud  radiance  in  the 

(antiparallel)  direction  of  the  LOS  i>  the  intersection  of  the  LOS  with 

the  cloud.  An  irradiance-to-radiance  transfer  coefficient,  TCOEF,  of 

dimension  sr”  ,  is  provided  by  the  cloud  routines  which,  when  multiplied 

by  the  normal  spectral  irradiance  at  the  transfer  point  on  the  cloud, 

* 

gives  the  required  spectral  radiance.  Transfer  point  t-^  in  Figure  2-3 
is  on  an  element  of  the  surface  which  is  visible  to  the  detector  D  but 
is  not  directly  illuminated  by  the  source  S.  It  is  thus  radiant  because 
of  an  outward  flux  of  photons  supplied  by  diffusion.  Assuming  diffuse 
radiance,  the  spectral  radiance  L^  is  given  by 

J+x(ti> 


where  J , 


is  a  function  of  all  the  irradiances  J 


's  over  the  illuminated 


facets  (k),  given  by 


cos  0 


The  transfer  coefficient  then  is 


TCOEFa  =  — g - -  Lx  (sr-1) 

X 


Similarly,  for  points  on  cloud  facets  which  are  both  illumi¬ 
nated  and  visible,  such  as  tg  in  Figure  2-3,  the  radiance  is  dominated 
by  directional  reflectivity.  By  definition  of  bidirectional  reflectance 
coefficient,  p',  the  required  spectral  radiance  is 


In  the  deterministic  model  (in  contrast  with  the  statistical  model) 
attenuation  is  not  included  along  the  air  path.  Such  attenuation  will 
be  computed  by  the  sight-path  integration  model. 


SXCOS0i  -2  -1 

L,  =  p'(0.,0  „,$)  - 5 -  (watt  km  sr  pm  1) 

A  Aif  4jrRJ 

st2 

The  transfer  coefficient  then  is  again  the  quantity  by  which  the  normal 
* 

irradiance  at  the  transfer  point  must  be  multiplied  to  give  radiance 
in  the  direction  of  the  detector  and  is  given  by 

^R2st 

TC°EFX  =  — g - -  Lx  =  (0.,0f,<t))  cos  0i(sr"1)  . 

X 

It  should  be  noted  that  the  bidirectional  reflectance  coefficient  is 
dependent  on  the  angle  of  incidence,  the  angle  of  reflection,  and  the 
azimuthal  angle  <t>  between  the  incidence  and  exit  planes.  This  is  the 
quantity  which  has  been  separately  calculated  by  Monte  Carlo  and 
tabulated  for  use  in  the  deterministic  cloud  problem. 

The  f ield-of-view  of  the  detector  is  one  of  the  specified 
quantities  in  the  problem.  The  transfer  coefficients  are  properly 
averaged  over  this  f ield-of-view  to  give  an  average  radiance  along  the 
line-of-sight .  If  the  f ield-of-view  is  less  than  the  facet  area 
projected  normal  to  the  line-of-sight,  then  the  facec  radiance  is  the 
required  average  radiance.  If,  however,  another  facet  lies  within  a 
distance  equal  to  0.7  times  the  square  root  of  the  area  of  the  f ield- 
of-view,  this  facet  also  contributes  to  the  average  radiance  and  its 
radiance,  weighted  by  the  appropriate  projected  facet  area  in  the  field- 
of-view,  is  included  in  the  calculation  of  the  average  radiance.  When 
part  of  the  f ield-of-view  does  not  intersect  the  cloud,  i.e.,  the  cloud 
fills  only  a  portion  of  the  f ield-of-view,  then  the  contribution  of  the 
off-cloud  area  to  the  average  radiance  is  zero  according  to  the  cloud 
model.  More  properly,  the  line-of-sight  integration  module  may  arrange 
to  provide  the  earth  upwelling  radiance  in  this  off-cloud  fraction. 

In  addition  to  the  transfer  coefficient  which  gives  radiance 
resulting  from  the  source,  the  cloud  as  an  emitter  of  radiation  must  be 
considered.  This  self-emission  radiance  is  calculated  by  assuming  the 
cloud  to  be  at  the  ambient  atmosphere  temperature  at  the  altitude  at 
which  the  LOS  intersects  the  cloud  and  to  emit  a  blackbody  spectrum  of 
that  temperature  with  an  appropriate  emissivity  in  the  direction  of  the 
* 

See  footnote  on  page  22. 


line-of-sight .  The  directional  emissivity  is  taken  to  be  one  minus  the 


directional  reflectance,  p 


dA' 


The  directional  reflectance  is  computed 


by  integrating  the  bidirectional  reflectance  coefficient,  over  the 

2ir  solid  angle,  holding  constant  the  incident  direction  equal  to  the 
angle  between  the  LOS  and  the  normal  to  the  surface.  That  is,  the 
emission  spectral  radiance  at  the  transfer  point  in  the  direction  of 


the  detector , 


Lex>  is  given  by 


LeX  ’  0-»dX>  3tt  V5 


C2^XT-1 


-1 


(watt/(km  sr  pm)) 


where 


c^  =  8irch,  c2  =  ch/k,  and 


PdA(9- 


± )  =y"  pj^  ( 6  i  ,  6 f ,  ^  )  cos0fdflf  . 


Since  this  radiance  is  independent  of  sources  in  the  problem,  it  is  not 
given  as  a  transfer  coefficient  but  is  calculated  as  a  spectral  radiance 
(watt/ (km“  sr  pm) ) .  Note  that  the  emissivity  of  each  of  the  14  cloud 
types  is  taken  to  be  that  computed  for  cloud-type  4. 

If  the  problem  being  examined  is  a  daytime  problem  so  that 
the  sun  also  illuminates  the  cloud,  a  third  source  of  cloud  radiance 
exists.  The  solar  reflected  radiance  is  calculated  exactly  as  though 
the  sun  were  another  point  source  and  the  transfer  coefficient  for  that 
geometry  is  the  output. 

2-2  STATISTICAL  CLOUD  MODEL 

The  ROSCOE-IR  statistical  cloud  model  is  based  on  the  cloud 
descriptions  and  cloud  occurrence  data  bank  of  Fowler  et  al. 

(Reference  1).  The  clouds  are  described  as  horizontal  plane  parallel 
slabs  with  a  base  altitude  and  a  top  altitude  given  for  each  cloud  type, 
together  with  water-droplet  average  densities  and  droplet-distribution 
average  radii.  The  occurrence  data  provide  the  probability  that  each 
such  cloud  type  will  occur  singly  or  in  multiple  layers  with  other 
types.  The  data  are  further  divided  into  coverage  categories  ranging 
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from  clear  (cover  category  1)  to  full  cloud  cover  (cover  category  5) 
approximately  by  quarters  (compare  definitions  of  CCOVER  in  Section  4-2), 
so  that  the  likelihood  of  a  LOS  being  cloudy  or  cloud-free  may  be 
determined . 

The  data  base  is  divided  into  the  30  geographic  regions  shown 
in  Figure  2-4.  Two  regions  of  interest  were  selected.  One  of  these  is 
region  11  of  the  global  data  bank  which  is  a  temperate  continental  region 
and  characterizes  the  northeast  United  States  and  western  Europe.  The 
second  region  is  region  4  which  is  a  tropical  maritime  region  and  includes, 
for  example,  the  Hawaiian  Islands.  For  each  region,  five  subsets  of  the 
total  data  base  have  been  selected  for  use  in  ROSCOE  to  provide  a  total 
of  10  user  options.  They  are,  for  each  region,  the  annual  average,  the 
winter  day  average,  the  winter  night  average,  the  summer  day  average,  and 
the  summer  night  average.  For  each  of  these  options  there  is  available 
the  probability  of  occurrence  of  each  of  the  14  cloud  types  singly  or  in 
multiple  layers,  for  each  of  the  four  cloud-coverage  categories,  and  a 
probability  for  cloudless  skies.  An  eleventh  option  permits  the  user  to 
construct  his  own  cloud-type  occurrence  probability  table,  it  being 
required  only  that  he  utilize  the  same  14  cloud  types  for  which  physical 
descriptions  exist.  This  option  permits  the  user  to  examine  the  effect 
of  any  one  of  the  cloud  types  and  permits  construction  of  cloud  scenes 
characteristic  of  the  information  contained  in  weather  reports  or  other 
cloud  descriptions. 

In  the  statistical  cloud  data  base  there  are  seven  types  of 
low-altitude  clouds,  three  middle  types,  and  four  high  types.  Allowing 
only  one  cloud  type  in  each  altitude  regime,  the  possible  configurations 
include  clear,  the  14  types  alone  as  a  single  layer,  61  combinations  of 
two  types  in  double  layers,  and  84  triple-layer  combinations  of  three 
types.  Thus  160  configurations  are  possible,  neglecting  the  possibility 
of  considering  the  four  fractional  cloud-coverage  categories  as  the 
basis  of  a  further  factor  of  four. 

This  last  consideration  is  avoided  by  treating  only  two  aspects 
of  each  configuration,  namely  the  possibilities  of  a  cloud-free  line-of- 
sight  (CFLOS)  along  the  ROSCOE  user's  LOS  and  of  a  LOS  obscured  by  clouds. 
In  the  latter  case,  the  clouds  are  assumed  to  be  extensive  in  horizontal 
dimension  so  that  the  radiation  transfer  is  the  same  as  for  plane  parallel 
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Figure  2-4.  Geographic  regions  of  the  NASA  worldwide  cloud  data  bank 


layers  of  infinite  extent.  This  is  consistent  physically  with  most 
fractional  coverage  cases,  i.e.,  the  cloud  patches  are  large  and  the 
gaps  are  large.  Thus  only  160  possible  configurations  must  be  addressed. 

For  the  ROSCOE-IR  statistical  cloud  model  it  was  required  to 
provide  spectral  radiance  quantities  along  the  selected  line-of-sight , 
namely  radiance  quantities  arising  from  NUDET  sources,  solar  reflection, 
and  cloud  emission.  In  order  to  retain  an  interface  with  the  main  pro¬ 
gram  similar  to  that  of  the  deterministic  cloud  model,  the  effects  of 
NUDET  and  solar  sources  were  to  be  returned  separately  as  transfer  coef¬ 
ficients,  with  units  of  sr  ^ ,  and  the  emission  was  to  be  returned 

-2  -1  —1 

separately  as  a  spectral  radiance  in  units  of  watt  km  sr  pm 

Since  the  number  of  possible  configurations  of  clouds  is  so 
large,  it  was  decided  not  to  transfer  from  the  main  program  to  the 
statistical  cloud  routine  and  back  again  for  each  possible  cloud  con¬ 
figuration  but  rather  to  loop  through  the  configurations,  each  having  its 
own  probability,  within  the  cloud  routine  and  to  return  to  the  main 
program  the  probability  distribution  of  radiance  results. 

The  probability  distribution  is  the  result  of  the  fact  that 
each  of  the  159  cloudy  configuration  cases  has  a  probability,  or  fre¬ 
quency  weight,  which  can  be  determined  from  the  ratio  of  the  frequency 
of  observed  cases  of  this  conf iguration  to  the  total  number  of  obser¬ 
vations.  These  frequencies  are  in  the  data  base. 

The  weights  of  each  configuration  are  adjusted  to  take  account 
of  the  probability  of  a  CFLOS.  The  weight  W(i)  for  the  LOS  intersecting 
the  i-th  cloudy  configuration  is  thus 

5 

W(i)  =  ^  W(i)r>cc  [l-CFL0Scc(9)] 
cc=2 


where  W(i)r  cc  is  the  raw  weight  of  the  i-th  configuration  by  type, 
layers,  and  coverage-category  cc  and  CFL0Scc(9)  is  the  probability  for 
a  CFLOS  at  the  zenith  angle  0  and  fractional  coverage  index  cc.  The 
CFLOS  probabilities  are  derived  from  Lund's  tables  of  CFLOS  probability 
for  all  cloud  types  combined  as  a  function  of  zenith  angle  and  fractional 
coverage  (Reference  3).  The  weight  assigned  to  the  clear  line-of-sight 
is  then 
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159 

W(clear)  =  P(cc=l)  +  [ 1-W( i ) ] 

i=l 

where  P(cc=l)  is  the  weight  of  the  clear  cloud-coverage  category  (no 
clouds)  and  the  sum  represents  the  gap  contributions  from  the  partial 
coverage  configurations.  A  more  detailed  treatment  of  these  equations 
is  given  in  the  discussion  of  subroutine  CLDWT  in  Section  3-2. 

The  statistical  cloud  routines  then  perform  the  appropriate 
calculation  of  radiation  transmission  for  159  configurations  of  cloud 
cover  and  for  clear  conditions  and  assign  to  each  the  probability 
weights,  W(160)  =  W(clear)  being  the  clear  1 ine-of-sight  weight.  The 
result  is  a  distribution  of  160  values  of  the  radiance-related  quantity 
calculated,  each  with  its  weight,  or  probability.  The  distribution  is 
summed  to  provide  an  integral  or  cumulative  distribution  and  radiance- 
related  values  from  this  cumulative  distribution  and  provided  as  output 
at  designated  optional  fractile  values  of  probability  provided  as  input 
in  the  call.  W( clear),  known  in  the  code  as  PCFLOS ,  is  also  provided 
as  an  output,  to  be  used  as  an  approximation  to  the  probability  of  a 
CFLOS  to  a  NUDET  source  within  the  detector's  f ield-of-view  and  below 
12-km  altitude. 

The  radiance-related  quantities  follow  the  same  system  as 
those  of  the  deterministic  cloud  model,  except  that  distributions  rather 
than  single  values  are  the  output.  For  NUDET  sources  and  for  the  sun, 
the  output  is  a  distribution  in  transfer  coefficients  which,  when 
multiplied  by  the  (unattenuated)  spectral  irradiance  at  the  transfer 
point  from  the  source,  results  in  a  distribution  of  spectral  radiances 
along  the  line-of-sight .  The  required  spectral  radiance  and  the 
transfer  coefficient  are  related  by  the  expressions 
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is  the 


where  S^  is  the  isotropic  point-source  spectral  power,  R 

source-to-transf er-point  distance,  and  E^  is  the  solar  irradiance.  The 

cloud  thermal  emission  output,  however,  is  an  actual  radiance  in  units  of 
-2  -1  -1 

watt  km  urn  sr  rather  than  a  coefficient,  because  the  thermal  emittance 
is  not  related  to  NUDET  or  solar  source  strengths  but  depends  solely  on 
cloud  temperature  and  emission  properties.  Again,  as  for  the  deter¬ 
ministic  case,  the  NUDET,  solar,  and  thermal  emission  outputs  are  pro¬ 
vided  separately. 

Figure  2-5  is  indicative  of  the  geometry  of  the  statistical 
cloud  problem.  One  or  more  cloud  layers  for  each  of  the  159  configu¬ 
rations  will  lie  somewhere  between  0.15-kilometer  and  12-kilometer 
altitude.  These  are  the  base  of  the  stratus  and  the  top  of  the  cumu- 
lonimbus-3  cloud  models,  respectively.  In  addition,  a  point  source  will 
be  located  somewhere  in  the  figure  such  as  at  SI,  S2,  or  S3.  The  source 
location  will  be  either  below  the  base  of  the  lowest  cloud,  somewhere 
in  the  region  occupied  by  one  or  more  cloud  types,  or  above  the  top  of 
the  highest  cloud.  A  method  has  been  developed  for  treating  each  of 
these  cases.  As  a  matter  of  fact,  during  the  progression  upward  of  a 
fireball  from  a  low-altitude  explosion,  the  source  will  progressively 
occupy  each  sort  of  location  relative  to  the  cloud  region. 

The  statistical  cloud  program  calculates  the  irradiance-to- 
radiance  transfer  coefficient  along  the  reverse  LOS  at  a  transfer  point, 
tp,  which  is  chosen  to  be  the  intersection  of  the  LOS  with  the  12-kilometer 
altitude  plane,  because  all  cloud  tops  lie  at  or  below  this  point.  By 
treating  the  atmospheric  transmission  problem  in  the  cloud  region  within 
the  cloud  subroutines,  as  discussed  later,  it  is  possible  to  avoid 
having  the  main  program  integrate  along  the  LOS  to  each  of  the  159  cloud- 
top  altitudes  in  the  entire  statistical  problem. 

In  addition  to  NUDET  sources  the  sun  as  a  radiation  source  is 
also  indicated  in  Figure  2-5.  For  each  of  the  159  cloud  configurations 
the  sun  illuminates  the  highest  cloud  top  of  the  configuration  and  by 
reflectance  contributes  to  the  radiance  along  the  LOS.  The  transfer 
coefficient  corresponding  to  this  solar  reflected  radiance  is  again  a 
distribution  of  160  points  and  is  provided  as  a  separate  output  distri¬ 
bution  from  that  of  a  NUDET  source.  Note,  however,  in  the  160-point 
distribution,  there  are  only  10  different  values  (nine  for  clouds,  one 
for  Earth)  because  (a)  there  are  only  nine  different  cloud-top  altitudes 
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tp( 12-km) 


Figure  2-5. 


Statistical  cloud  geometry 


(cf.  Table  2-1)  and  (b)  the  bidirectional  reflectance  of  any  cloud  is 
assumed  to  be  approximated  by  that  for  cloud-type  4  (stratocumulous) . 

In  addition  to  the  LOS  intersecting  cloud  tops  in  the  159 
cloud  configurations,  there  is  the  probability  of  a  CFLOS  to  the  Earth's 
surface,  which  acts  as  a  solar  reflector  and  as  an  emitter.  The  surface 
may  be  sunlit  or  not,  the  probabilities  depending  on  the  probability  of 
a  CFLOS  to  the  surface  point  from  the  sun.  The  contribution  of  the 
surface  reflectance  is  weighted  by  the  product  probability  of  a  two-leg 
CFLOS . 

Finally,  as  in  the  case  of  the  deterministic  clouds,  the 
statistical  cloud  tops  are  emitters  which  contribute  to  the  radiance 
along  the  LOS.  The  cloud  top  for  each  configuration  is  treated  as  a 
blackbody-spectrum  emitter  at  the  temperature  of  the  cloud  top  and  with 
directional  emissivity  equal  to  (1-p^)  as  in  the  deterministic  case. 

(The  directional  emissivity  computed  for  cloud-type  4  (stratocumulous) 
is  used  for  all  cloud  types.)  This  emission  radiance  distribution  is 
again  carried  separately  from  the  NUDET  source  and  the  solar  distributions. 
The  surface  is  again  the  160th  point  in  this  distribution.  As  in  the 
solar  reflectance  case,  there  are  only  10  different  values  in  the  160- 
point  distribution. 

A  summary  guide  to  the  distributions  available  from  the 
statistical  cloud  model  is  given  in  Table  2-2. 

Figure  2-6  is  a  schematic  of  one  of  the  calculations.  A  point 
source  lies  below  a  plane  parallel  cloud.  By  the  solution  of  the  dif¬ 
fusion  problem  in  cylindrical  geometry,  described  in  Section  A-2.1.3J 
the  spectral  flux  leaving  the  cloud  top  is  calculated  and  has  a  bell¬ 
shaped  radial  distribution  as  indicated  in  (a).  The  outgoing  flux 
(radiant  exitance  or  emittance)  at  the  cloud  top,  J+,  is  given  by  a  sum 
of  weighted  Bessel  functions  (Equation  (30a),  with  Q  =  S/2)  and  the 
radiance,  L,  is  just  this  flux  divided  by  ir,  since  the  angular  distri¬ 
bution  is  taken  to  be  diffuse.  The  radiance  at  the  point  where  the  LOS 
intersects  the  cloud  top  is  calculated  and  divided  by  the  cloud-free 
normal  irradiance  at  the  transfer  point,  TP,  to  give  a  transfer  coef¬ 
ficient  with  units  of  sr”1  as  in  the  deterministic  case.  That  is, 
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Table  2-2.  Guide  to  distributions  from  the  statistical 
cloud  model  . 


Thermal 

Solar 

Artificial 

LOS 

Emission , 

Source 

Source 

Terminates 

Radiance 

Transfer  . 

Transfer 

IDX 

On 

W/(km2  sr  pm) 

Coef f . ,  sr" 

Coef f . ,  sr- 

1  to  159 

Clouds 

a 

c 

e 

160 

Earth 

b 

d 

d 

aEach  emission  in  the  set  of  159  is  from  the  top  of  the  uppermost 
cloud  layer  and  is  attenuated  along  the  path  to  the  transfer  point 
at  12-km  altitude.  There  are  nine  distinct  values  among  the  159 
values  of  such  radiance  corresponding  to  the  nine  different  cloud- 
top  altitudes. 

bThe  emission  is  from  the  Earth's  surface  (the  nature  of  which  can 
be  specified)  and  is  attenuated  along  the  path  to  the  transfer  point 
at  12-km  altitude. 

Each  of  the  transfer  coefficients  in  the  set  of  159  corresponds  to 
the  bidirectional  reflectance  of  the  uppermost  cloud  layer  for  the 
incident-  and  ref lected-photon  geometry  and  includes  air  transmittance 
along  those  portions  of  the  incoming  and  reflected  paths  below 
12-km  altitude.  The  nine  distinct  values  among  the  159  values  of 
such  transfer  coefficients  correspond  to  the  nine  different  cloud- 
top  altitudes. 

^The  transfer  coefficient  corresponds  to  the  bidirectional  reflectance 
of  the  Earth's  surface  for  the  incident-  and  ref lected-photon 
geometry  and  includes  air  transmittance  along  the  incoming-  and 
ref lected-photon  paths  below  12-km  altitude.  The  weight  for  this 
case  includes  the  (product)  probability  of  a  two-leg  CFLOS  instead 
of  a  one-leg  CFLOS  as  for  the  thermal  emission  case. 

A  given  transfer  coefficient  may  correspond  to  either  (a)  trans¬ 
mittance  through  cloud  layers  between  the  source  and  the  transfer 
point  or  (b)  bidirectional  reflectance  from  the  uppermost  cloud 
layer  below  the  source.  Thus  the  distribution  will  include  a 
combination  of  the  above  two  types  of  transfer  coefficients  for 
a  source  above  the  lowest  cloud-top  altitude  (1.0  km)  and  below 
the  highest  cloud-top  altitude  (12.0  km).  For  each  transfer 
coefficient  one  includes  air  transmittance  along  those  portions 
of  (nominal)  paths  below  12-km  altitude. 


Figure  2-6. 


Typical  statistical  cloud  geometry.  (a)  shows 
the  bell-shaped  radial  distribution  of  flux, 
J+,  leaving  the  cloud  top  and  (b)  shows  the' 
source-to-transfer-point  distance  Rst  and  the 
air  transmittance  path  S-CT-TP. 


Since,  for  statistical  clouds  in  contrast  with  deterministic 
clouds,  the  main-program  integration  of  radiance  along  the  LOS  does  not 
extend  below  12  kilometers,  the  air  transmission  reduction  of  the 
source  radiance  is  approximated  within  the  cloud  model  and  reduces  the 
cloud-top  radiance,  ,  prior  to  calculation  of  the  transfer  coefficient. 
The  air  path  from  the  source  to  the  point  of  intersection  of  the  LOS 
with  the  cloud  top,  CT,  and  from  this  point  to  the  transfer  point  is 
used  for  this  air-effect  calculation  as  shown  in  Figure  2-6(b).  If  the 
source  is  above  the  cloud  top,  the  air-transmittance  path  is  again  from 
the  source  S  to  the  LOS  intersection  point  with  the  cloud  top  (or  that 
portion  of  the  path  below  12-km  altitude)  and  then  to  the  transfer 
point  TP.  For  emission,  the  air  path  along  the  LOS  is  from  the  point 
CT  on  the  cloud  top  to  the  transfer  point  TP.  The  air  transmittance 
=  Tsct  Tcttp^  has  t)een  calculated  by  use  of  the  LOWTRAN  code 
(Reference  4 ) . 

The  geometries  used  in  the  diffusion  calculations  are  shown 
in  Figure  2-7.  Diffusion  solutions  have  been  obtained  for  the  case  of 
a  point  source  beneath  a  cloud  (a),  a  point  source  within  a  single¬ 
layer  cloud  (b),  and  for  a  point  source  at  or  below  the  base  of  a 
multi-layered  configuration  of  any  number  of  layers  (c).  These  three 
cases  handle  the  necessary  configurations  with  two  exceptions.  First, 
the  cloud-type  descriptions  do  not  necessarily  result  in  the  base  of 
an  upper  cloud  coinciding  with  the  top  of  a  lower  cloud  so  that  gaps 
or  overlaps  exist  between  multiple  layers.  When  this  occurs  the  base 
of  the  higher  layer  is  adjusted  to  coincide  with  the  top  of  the  layer 
below,  the  optical  thickness  of  the  upper  layer  being  retained  by 
adjustment  of  the  particle  number  density.  In  this  manner  a  continuous 
multi-layer  medium  is  created  which  is  compatible  with  the  available 
diffusion  solutions.  Second,  when  the  source  lies  within  a  cloud 
member  of  a  multiple  layer  as  shown  in  Figure  2-7(d),  the  cloud  region 
below  the  source  is  neglected  (truncated)  so  that  the  problem  is 
reduced  to  that  of  Figure  2-7(c)  with  the  source  on  the  bottom  surface. 
(The  two  adjustments  are  not  commutative.  In  the  code,  we  have  imple¬ 
mented  the  second  adjustment  before  the  first.)  The  truncation  modi¬ 
fication  is  expected  to  cause  about  a  10%  reduction  in  the  cloud-top 
radiance,  this  being  the  approximate  magnitude  of  the  albedo*  for  these 
♦ 

See  function  RHOEPS  for  a  tabulation  at  10  wavelengths  of  the 
hemispherical-directional  reflectance,  or  equivalently,  directional- 
hemispherical  reflectance. 
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Figure  2-7.  Diffusion  and  reflectance  geometries  for  statistical 
clouds.  Diffusion  solutions  have  been  obtained  for 
Cases  (a),  (b),  and  (c).  Case  (d)  is  treated  approxi¬ 
mately  by  truncating  the  cloud  region  below  the 
source,  thereby  converting  it  to  Case  (c)  with  source 
on  surface.  The  solution  for  the  transmitted  photon 
flux  is  given  by  Equation  (30a)  for  Case  (a),  by 
Equations  (26b)  and  (40a)  for  Case  (b),  and  by 
Equations  (69)  and  (70)  for  Case  (c).  Note  that 
Cases  (a)  and  (b)  become  identical  for  the  source  on 
the  surface  (h=0).  In  this  regard,  one  can  verify 
that  for  h=0 ,  Equation  (40a)  becomes  identical  to 
Equation  (29a)  provided  Q  in  Equation  (29a)  is  set 
equal  to  one-half  of  the  isotropic  source  strength,  S. 
In  Equation  (40a),  S  is  tacitly  unity. 
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clouds  in  the  two-  to  five-micron  region.  A  reflecting  layer  of  albedo 
w  below  the  source  increases  the  cloud-top  radiance  by  a  factor  of 
about  l/(l-u>) ,  as  shown  at  the  end  of  Appendix  D. 

Figure  2-8  shows  the  assumed  altitude  range  of  each  of  the 
14  cloud  types  for  single-layer  configurations.  The  dashed  lines  at 
the  lower  edges  of  the  mid-  and  high-altitude  clouds  depict  the 
boundaries  that  are  moved  to  achieve  layer  contiguity  in  multilayer 
conf igurations . 

For  the  source  above  the  cloud  top  as  shown  in  Figure  2-7(e), 
bidirectional  reflectance  dominates  the  radiance  along  the  LOS  and  is 
calculated  from  the  bidirectional  reflectance  table  (Function  CLDBDR) 
as  was  done  for  the  deterministic  case.  The  transfer  coefficient  is 
then 
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Here  the  air  transmittance  is  taken  to  be  the  product  of  the  trans- 
mittances  and  along  those  portions  of  the  incoming  and  reflected 
paths  below  12-km  altitude. 

Examples  of  statistical  cases  are  given  in  Section  4  and  show 
both  the  radial  distribution  of  radiance  across  the  cloud  tops  and 
examples  of  the  probabilistic  distributions  in  radiance  arising  from 
the  many  possible  cloud  configurations. 

Through  use  of  the  statistical  cloud  model  the  user  may  obtain 
three  separate  distributions  for  spectral  radiance  from  NUDET  sources, 
from  the  solar  source,  and  from  thermal  emission  backgrounds.  It  is 
rather  obvious  but  perhaps  worthy  of  note  that  these  distribution  data 
are  more  complex  quantities  than  would  be  the  three  simple  radiance 
values  from  a  completely  deterministic  case.  The  three  simple  radiances, 
for  example,  can  be  added  to  give  the  total  radiance  which  would  be 
observed  by  a  sensor.  They  can  be  divided  to  give  ratios  of  signal- 
to-noise,  signal-to-solar  background,  etc.  Such  quantities  are  often 
important  measures  of  sensor  effectiveness. 
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Figure  2-8.  Altitude  ranges  of  single-layer  cloud  types. 


Such  simple  sets  of  three  radiance  values  are  obtained  from 
the  deterministic  cloud  model  and  also  may  be  obtained  from  the  sta¬ 
tistical  model  if  used  in  the  Option-11  mode  so  that  a  single  cloud 
configuration  is  specified. 

Each  of  the  distributions  from  the  statistical  cloud  model 
results  from  160  such  deterministic  configurations.  However,  once  the 
distribution  has  been  constructed,  the  contributions  of  the  individual 
configurations  are  obscured  and  the  origin  of  any  part  of  the  distri¬ 
bution  loses  its  identity,  unless  special  steps  are  taken. 

It  is  not  true,  for  example,  that  the  lowest  or  highest  ten 
percent  of  NUDET  source  radiances  result  from  the  same  cloud  configu¬ 
rations  as  those  that  result  in  the  lowest  or  highest  solar  background 
radiances.  As  a  matter  of  fact,  even  the  solar  background  and  the 
thermal  emission  background  distributions  cannot  be  correlated  in  such 
a  simple  manner. 

Some  kinds  of  simple  generalizations  may  be  made  when  the 
radiance  values  in  the  distributions  permit.  For  example,  if  the 
highest  background  radiance  value  corresponds  to  the  tenth  percentile 
value  for  the  source,  one  can  conclude  that  in  90  percent  of  the  cases 
the  signal-to-noise  ratio  is  greater  than  unity.  If  the  50th  percentile 
value  for  the  source  is  six  times  the  maximum  background  radiance  value, 
one  knows  that  in  at  least  50%  of  the  cases  the  signal-to-background 
ratio  is  greater  than  six.  Generally  speaking,  however,  the  user 
cannot  simply  add  or  divide  fractiles  of  the  distributions  to  obtain 
meaningful  quantities. 

If  the  distribution  in  such  quantities  as  signal-tornoise 
ratio  of  solar-to-thermal  ratio  are  of  prime  interest,  and  if  the  dis¬ 
tributions  of  signal  and  background  overlap  appreciably  in  radiance 
values,  special  treatment  of  the  statistical  cloud  model  may  be 
warranted.  It  may  be  most  convenient,  for  example,  to  do  a  program 
modification  to  identify  and  store  the  individual  elements  of  the  three 
distributions  so  that  they  may  be  processed  configuration  by  config¬ 
uration  to  obtain  the  probability  distribution  in  the  quantity  of 
interest . 
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SECTION  3 

COMPUTER  PROGRAM  DESCRIPTION 


3-1  INTRODUCTION 

The  main  (or  driver)  program  ROSCOER  is  intended  to  demon¬ 
strate  the  utility  of  the  ROSCOE-IR  cloud  computation  procedure  and  to 
facilitate  integration  of  the  subroutines  into  ROSCOE-IR.  All  require¬ 
ments  of  the  computational  procedure  are  met  in  ROSCOER  to  provide  a 
stand-alone  code  and  to  provide  a  guide  to  the  intended  use  of  the  sub¬ 
routines.  All  overall  program  flow  is  illustrated  in  Figure  3-1. 

Data  provided  by  ROSCOE-IR  (such  as  source,  detector,  and  sun 
location)  are  read  by  the  driver  program  ROSCOER  or  are  entered  via  data 
statements.  Calls  to  the  deterministic  and  statistical  cloud  sub- 
modules  are  illustrated  via  ROSCOER.  There  is  no  preferential  looping 
on  lines-of-sight ,  wavelength,  or  cloud  number.  The  calling  program, 
ROSCOER,  simply  demonstrates  how  the  code  can  be  used.  The  determin¬ 
istic  check  case  (Problem  6  in  Section  4)  loops  over  the  64  possible 
lines-of-sight*  at  the  second  wavelength  in  the  wavelength  table,  i.e., 

X  =  2.5  ym.  A  check  case  for  the  statistical  cloud  submodule  (Problems 
1  to  4 ,  Section  4)  loops  over  as  many  as  161  lines-of-sight  at  one- 
quarter  kilometer  intervals  from  the  source  normal  to  the  lines-of- 
sight.  The  transfer  coefficient  is  integrated  with  proper  conversions 
to  give  a  total  transmission 

Several  parameters  are  initialized  by  statements  in  the  main 
program  only  as  a  means  of  checking  the  code  calculations  and  demon¬ 
strating  its  utility.  The  user  will  necessarily  obtain  this  information 
by  a  replacement  of  this  main  program.  A  call  to  subroutine  CLOUDO 
obtains  data  that  are  required  only  by  the  cloud  package  and  establishes 
the  mode  to  be  that  of  either  a  statistical  cloud  data  base  or  a  set  of 
deterministic  clouds.  Deterministic  cloud  properties  are  still  obtained 
from  the  preprogrammed  data  base;  however,  sizes  and  shapes,  within  the 
ellipsoidal  constraint,  and  locations  are  to  be  provided  by  the  user. 

*The  number  of  lines-of-sight  returned  by  subroutine  LOS  has  been 
reduced  from  a  maximum  of  64  to  a  maximum  of  16  by  incrementing  the 
facet-center  indexes  by  2’s,  a  procedure  that  appears  to  be  adequate. 
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The  transfer  coefficients  for  the  calculations  will,  of 
course,  be  multiplied  by  a  source  irradiance  to  obtain  a  radiance  in 
the  direction  of  the  detector  at  the  transfer  altitude  of  12  kilometers 
(for  statistical  clouds).  This  coefficient,  as  well  as  the  emission, 
has  included  in  it  the  atmospheric  transmission,  in  both  directions, 
if  appropriate,  from  the  1962  U.S.  Standard  Atmosphere  as  calculated 
from  LOWTRAN  3  (Reference  4).  No  atmospheric  attenuation  is  included 
in  the  deterministic  case  because  it  will  be  computed  by  the  line-of- 
sight  integration  module  in  ROSCOE-IR. 

Parameters  such  as  lines-of-sight  and  wavelength  are  also 
user  defined.  A  table  of  wavelengths  is  provided  in  common  CDDATA  at 
which  cross  sections  are  tabulated.  This  list  provides  a  set  of  possible 
wavelengths  to  check  out  the  code  and  to  make  specific  calculations. 

The  only  restriction  on  wavelength  is  that  it  be  between  two  and  five 
microns . 

A  line-of-sight  must  be  provided  to  both  the  deterministic 
and  statistical  subroutines.  Possible  lines-of-sight  are  provided  for 
deterministic  clouds  by  a  call  to  subroutine  LOS.  These  are,  in  turn, 
available  in  common  LOS.  A  maximum  of  64  lines-of-sight  are  returned. 
This  table  is  terminated  by  zeros  in  the  altitude  array  HLOS .  Any  or 
all  of  these  lines-of-sight  are  appropriate  for  subsequent  deter¬ 
ministic  calculations. 

Since  a  statistical  cloud  is  infinite  in  lateral  extent, 
there  is  no  mechanism  to  relate  the  detector  line-of-sight  to  the  cloud. 
One  usually  has  a  particular  scenario  which  directs  which  way  the 
detector  is  looking  and  it  is  this  consideration  which  gives  rise  to 
the  input  for  statistical  clouds.  There  is,  however,  little  reason  to 
make  transmission  calculations  for  horizontal  distances  in  excess  of 
five  times  the  source  to  cloud-top  distance.  This  situation  may  occur 
in  statistical  cases  where  the  transmission  is  negligible;  however,  the 
code  doesn't  bother  to  eliminate  the  calculation. 

The  routines  called  to  get  an  answer  are  CL0UD3  (UL,VL,WL,K, 
ALAM , TCOEF , EMISS , SINCLD)  for  deterministic  clouds  and  SCLOUD(ALAM,UL, 
VL,WL,ISUN)  for  statistical  clouds,  where  UL,VL,WL  are  the  direction 
cosines  of  the  LOS,  ALAM  is  the  wavelength  in  microns,  TCOEF  is  the 
irradiance-to-radiance  transfer  coefficient  at  the  intersection  point 
of  the  LOS  and  cloud-K,  EMISS  is  the  thermal  emission  from  the  cloud 
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at  the  intersection  point,  and  ISUN  is  a  flag  to  distinguish  between  an 
artificial  source  and  the  sun.  Common  PTE  variables  must  be  initialized 
and  utilized  as  follows: 


( Input ) 

P1,P2,P3,P4  -  four  probability  levels  (fractiles)  for 
which  transfer  coefficients  and  thermal 
emission  are  required. 

(Output) 

T1,T2,T3,T4  -  four  transfer  coefficients  corresponding 
to  the  P1,P2,P3,P4  above. 

E1,E2,E3,E4  -  thermal  emission  corresponding  to  PI, 
P2,P3,P4  above. 


S1,S2,S3,S4  -  four  transfer  coefficients  corresponding 
to  P1,P2,P3,P4  above  for  the  sun  as  a 
source . 


PCFLOS  -  probability  of  the  detector  seeing  the 

artificial  source,  either  with  or  without 
clouds  being  present,  provided  the  source 
is  within  the  f ield-of-view. 


3-2  DESCRIPTION  OF  SUBROUTINES 


A  summary  of  the  information  necessary  to  execute  each  of  the 
subroutines  is  provided  in  this  section.  Some  of  the  subroutines  are 
documented  to  a  greater  extent  than  otners  as  the  use  of  many  is  self- 
evident  . 


3-2 . 1  Program  ROSCOER 

The  driver  program  ROSCOER  is  used  to  provide  test  cases  for 
the  deterministic  and  statistical  cloud  submodules.  Data  are  read  for 
the  locations  of  the  source  and  detector.  A  call  to  CLOUDO  reads  the 
data  which  are  required  by  the  cloud  module.  Branching  is  directed  by 
a  mode  flag  to  the  submodules. 

The  deterministic  mode  is  illustrated  by  a  call  to  subroutine 
LOS  and  subsequent  calls  to  subroutine  CL0UD3  for  the  lines-of-sight 
determined  by  LOS.  The  statistical  mode  is  illustrated  by  a  series  of 
lines-of-sight  providing  a  function  which  is  subsequently  integrated 
across  the  cloud  top  to  give  total  transmission. 
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3-2.2  Subroutine  ATMOSU( JJ.ZH) 

This  subroutine,  a  dummy  routine  for  the  real  one  (see  the 

3 

ROSCOE  Manual,  Volume  14a-l),  provides  the  atmospheric  density  (g/cm  ) 
and  temperature  (degrees  K)  as  a  function  of  altitude  (km).  The  common 
ATMOUP  is  utilized  with  the  fifth  and  sixth  locations  being  the  density 
and  temperature.  The  calling  sequence  is  (2,ZH),  ZH  being  the  altitude. 
The  value  of  2  for  JJ  denotes  an  operational  call,  allowed  after  ATMOSU 
has  been  initialized  earlier  with  JJ=1 .  The  density  and  temperature 
are  provided  by  simple  expressions: 

RHO  =  1.225E-3  *  EXP(-ZH/9. ) 

T  =  288.  -ZH  *  6.5 

3-2.3  Function  BESSO(X) 

This  routine  provides  the  zero-order  cylindrical  Bessel 
function  per  formulas  9.4.1  and  9.4.3  in  the  NBS  Handbook  (Reference  5). 

3-2.4  Block  Data 

This  BLOCK  DATA  set  contains  the  extinction  (SIGT( 10 , 14) )  and 

scattering  (SIGS( 10 , 14) )  cross-sections  and  the  average  cosine  of  the 

scattering  angle  (XMUBAR( 10 , 14 ) )  for  the  14  cloud  types.  Units  for  the 

cross  sections  are  square  microns;  cloud  droplet  number  density  is  in 
3  -1 

reciprocal  cm  .  The  code  uses  km  as  macroscopic  cross-section  units 
with  the  conversion  from  microscopic  cross-sections  being 

ax(km_1)  =  ox(ym2)  x  N(cnT3)  x  ±9-  cm2  x  lo5  . 

ym 

SIGT,  SIGS,  and  XMUBAR  are  tabulated  at  10  wavelengths  from  2  to  5 
microns.  The  entries  were  calculated  by  Mie  theorya  using  cloud  pro¬ 
perties  from  Fowler  (Reference  1)  and  given  in  Table  2-1. 

Also  tabulated  are  cloud  coverage,  CCOVER(5 , 11 ) ,  and  fre¬ 
quency,  CFREQ( 17 , 4 , 11 ) ,  data  for  the  statistical  module. 

aSee  the  comment  at  the  end  of  Section  3-2.16  for  an  exception  to 
this  statement. 
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CCOVER( 1 ,K)  is  frequency  of  clear  sky 

(2,K)  is  frequency  of  cloud-coverage  category  2 

(3,K)  is  frequency  of  cloud-coverage  category  3 

(4,K)  is  frequency  of  cloud-coverage  category  4 

(5,K)  is  frequency  of  cloud-coverage  category  5 

where  K  is  the  model  (i.e.,  time,  region)  discussed  in 
Section  2-2, 

CFREQ( I , J , K) ,  1=1,17;  J=l,4:  K=l,ll 


CFREQ(1 , J,K) 

is 

frequency 

of 

Cul 

( 2 ,  J ,  K) 

is 

frequency 

of 

Cu2 

( 3  ,  J  ,  K) 

is 

frequency 

of 

Cu3 

(4,J,K) 

is 

frequency 

of 

Sc 

( 5 , J , K) 

is 

frequency 

of 

St 

(6,J,K) 

is 

frequency 

of 

Cbl 

( 7  ,  J ,  K) 

is 

frequency 

of 

Ns 

(8 , J ,K) 

is 

frequency 

of 

Ac 

(9,J,K) 

is 

frequency 

of 

Asl 

( 10, J ,K) 

is 

frequency 

of 

Cb2 

(11 , J,K) 

is 

frequency 

of 

Cb3 

( 12 , J , K) 

is 

frequency 

of 

As2 

( 13 , J , K) 

is 

frequency 

of 

Ci 

( 14 , J , K) 

is 

frequency 

of 

Cs 

14 

where  ^  CFREQ(I,J,K)  =  1 
1=1 


CFREQ( 15 , J , K) 

is 

frequency 

of 

1 

layer 

(16, J,K) 

is 

frequency 

of 

2 

layers 

( 17, J,K) 

is 

frequency 

of 

3 

layers 

17 

where  ^  CFREQ(I,J,K)  =  1 
1=15 


J  is  the  cloud-coverage  index  and  K  is  the  model  number. 


3-2.5 


Function  CFLOSF( ICC , ANGLE) 


This  function  obtains  a  value  for  the  probability  of  a  cloud- 
free  line-of-sight  from  a  table  by  Lund  and  Shanklin  (Reference  3). 

The  array  is  tabulated  for  cloud  coverage  (in  lOths)  and  for  zenith 
angles  in  10-degree  increments  from  0  to  90°.  The  table  is  for  all 
cloud  types  combined. 

The  calling  sequence  requires  the  cloud-coverage  index  ICC  to 
be  an  integer  in  lOths  and  the  angle  in  degrees. 

INPUTS: 

ICC-1  -  cloud  coverage  (lOtns) 

ANGLE  -  zenith  angle  in  degrees 

OUTPUT: 

CFLOSF  -  probability  of  a  cloud-free  line-of-sight 
at  a  given  zenith  angle  for  the  requested 
cloud-coverage  index  ICC 

Cloud  coverages  reported  in  tenths  were  grouped  into  five 
caregories  by  the  authors  of  Reference  1  (p.3): 


Coverage  category  1 
Coverage  category  2 
Coverage  category  3 
Coverage  category  4 
Coverage  category  5 


0  tenths 
1,2,3  tenths 
4,5  tenths 
6, 7, 8, 9  tenths 
10  tenths 


In  our  cloud  model  the  decimal  cloud-coverage  fractions  of  0.0,  0.3, 
0.5,  0.8,  and  1.0  are  assigned  to  the  cloud-coverage  categories  1,2, 3, 4 
and  5,  respectively. 

3-2.6  Function  CLDBDR  ( CTHIN , CTHOUT , CPHOUT , ALAM , KCLOUD) 

This  function  contains  tables  of  bidirectional  reflectances 
for  a  stratocumulus  cloud.  The  data  are  tabulated  in  forty-five  equal 
solid  angle  bins  at  10  wavelengths  from  2  to  5  microns.  The  bi¬ 
directional  reflectance  has  units  of  sr’1.  In  the  calling  sequence 
the  variables  are: 
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CTHIN  -  cosine  of  incident  zenith  angle 

CTHOUT  -  cosine  of  exit  zenith  angle 

CPHOUT  -  cosine  of  the  azimuthal  angle  between  the 
incident  and  exit  planes 

ALAM  -  wavelength,  pm 

KCLOUD  -  not  used. 

A  linear  interpolation  is  performed  between  wavelengths  but 
no  angular  interpolation  is  utilized. 

3-2.7  Subroutine  CLDGEOM 


Required  bookkeeping  for  cloud  "facets"  is  accomplished  in 
this  routine.  See  Figure  3-2  for  a  flow  chart.  A  maximum  of  20  deter¬ 
ministic  clouds  can  be  considered  in  a  given  job.  Areas  of  the  facets, 
radii  to  their  centers,  coordinates  of  the  boundaries,  and  coordinates 
of  the  centers  are  evaluated  and  stored  in  common  CDCORD  for  use  by 
other  routines.  There  is  no  requirement  for  user  interaction  with  this 
subroutine,  unless  he  wants  to  consider  the  irradiance  on  a  particular 
facet  and  facet  orientation  to  the  source  and  detector.  To  do  so  would 
require  a  code  modification. 

Coordinates  of  the  facets  are  computed  according  to  the 
prescription : 


9  _  JL  _ l _ 

2  2  1  +  (Rl)a  +  (Rl)2a  +  (R1 )3a 
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<t>3  =  |  -  [!  +  (R2)3]  A 


where  R1  is  c/a,  R2  is  b/a,  and  a,  b,  and  c  are  the  semi-axis  lengths. 
These  coordinates  are  reflected  in  the  eight  octants;  consequently,  they 
are  calculated  only  for  the  first  octant  and  are  mapped  to  the  other 
seven  by  the  function  IPRIME.  This  is  also  true  of  areas  and  radii. 

Areas  are  evaluated  by  a  numerical  quadrature  (Simpson's  rule) 
using  twenty  intervals  in  both  0  and  <|>  for  each  facet.  Differential 
area  elements  are  calculated  in  DSURF .  All  data  obtained  from,  and 
returned  to,  common  arrays. 

3-2.8  Subroutine  CLDWT(ALAM, ANG) 

Weightings  for  all  possible  cloud  and  cloud-layer  combinations 
are  evaluated  in  this  routine.  These  are  subsequently  assigned  to  the 
transfer  coefficients  and  values  of  the  emission  as  calculated  by  calls 
to  TRANSF  and  EMISSF. 

The  correspondence  between  the  notation  used  in  Section  2-2 
and  that  used  in  the  following  equations  is  given  below: 


Section  2-2 


Here 


W( clear ) 
P(cc=l) 

1  'SO 

72  tl-W(i)] 
i=l 

W(i),  i=l , 159 

CFLOS  ( e ) 
cc 


WTnewdBO) 

CCOVER( 1 ) 


WT  .  .(160 
old 


) 


WT(IDX),  IDX=1 , 159 
CFLOS 

c 
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For  any  of  the  10  location-season  averaged  cloud  models 
( KMODEL  =  1,10),  WT(IDX)  is  the  probability  that  (a)  the  cloud- 
configuration  set  indicated  by  the  index  IDX  occurs  and  (b)  the  LOS 
at  zenith  angle  6e  intersects  the  cloud-configuration  set.  The  total 
probability  of  the  LOS  intersecting  clouds  is 

159 

^  WT(IDX)  . 

IDX  =  1 

The  probability  that  the  LOS  does  not  intersect  the  clouds  but  "sees 
through"  the  clouds  is 

WT0id  (160)  . 

The  probability  of  seeing  the  ground  without  any  clouds  being  present 
is 

CCOVER  ( 1 )  . 

(For  brevity  we  have  suppressed  the  index  KMODEL  on  which  the  array 
CCOVER  depends.)  The  probability  of  seeing  the  ground  either  with  or 
without  clouds  being  present  is 

WTnew  (160)  =  WTold  (160)  +  CC0VER  (!)  • 

The  probability  of  seeing  either  clouds  or  the  ground  (with  or  without 
clouds  being  present)  is 

159 

Y,  WT(IDX)  +  WTnew  (160)  =  1  . 

IDX  =  1 


•19 


WT( I )  =  ^  CCOVER( C+l )  (l-CFLOSc) 
C=1 


CF(L,C) 
1L  TL 


L=1 , 7  1=1,7 


CF(15,C)«  {  TM  •  CF(7t+m-,C1 


M=  1 , 3  1  =  8,10 


TH 


CF(10+H,C) 

TI1 


H=1 , 4  1=11,14 


TL  •  TM  •  CF(L,C)  »  CF(7+M,C) 
DEN  TL  •  TM 


L=1 , 7  1=15,35 

M=1 ,3 


CF(16,C)«  (  TL  *  TH  CF(L,C)  »  CF(10+H,C) 
DEN  TL  •  TH 


L=1 , 7  1  =  36,63 

H=1 ,4 


TM  •  TH  CF( 7+M,C)  •  CF(10+H,C) 
DEN  TM  •  TH 


M=1 , 3  1=64,75 

H=1 ,4 


CF( 1 7 , C)  •  1  • 


CF( L , C)  >  CF( 7+M ,C)  «  CF(10+H,C)  L=l,7  1  =  76,159 


TL 


TM 


TH 


M=1 , 3 
H=1 ,4 


where 


CF  =  CFREQ 
TL  =  TOTL 
TM  =  TOTM 
TH  =  TOTH 

DEN  =  TL  •  TM  +  TL  •  T1I  +  TM  •  TH 


Here,  U(j | £=1 ,C)  is  the  probability  of  layer  j,  given  one  layer 
and  cloud-coverage  index  C 

V(j,r|*-2,C)  is  the  probability  of  layers  j  and  j',  given 
two  layers  and  index  C 

W( 1 , 2 , 3 | £=3 ,C)  is  the  probability  of  layers  1,2,3,  given 
three  layers  and  index  C 
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These  probabilities  satisfy  the  following  relations: 


3  3  KMAXj 


U(j|M,C)  =  T,  = 


£  £ 


CF( K | j , C)  =  1 


j=l  j=l  K=KMINj 


££t  T . 

,j'\  1=2  ,C)  =  j - -1-  =  1;  W(1,2,3|A-3,C)  =  1 


2  3  >  > 

2  2  v(j’j'iz=2’c)  = 
j  =  l  j'>j  7  j7  4 


T  .  T  .  . 
J  J 


The  probability  of  seeing  clouds  =  ^  ]  WT(I) 


=  ^  CCOVER  (C+l)  •  (l-CFLOSc)  ^  CF(14  +  H,C) 


KMAX  . 

A  CF(K|  j  ,C) 

2^  U(j|A-l,C)  2^  Tj  *=1 

j=l  K=KMIN  . 

tJ 


2  3 

^2  ^2  v(j  j"u=2.c) 

j=i  j'>j 


KMAXj  KMAX . . 

^  CF(K|.j,C)  CF(K^ir.C) 

K=KMIN  .  K '=KMIN . -  j  J 

J  J 


7  10  14 

W(l,2,3|Jl=3,C)  £  £  — (  Kt1-X-  ’ C)-  ;*  CF(  K  ^  1  2  ,  C  )  »  CF(  K  -  [.3 ,  C  ) 


K=1  K'=8  K"  =  ll 
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4 

=  ^  CC0VER(C+1 )  •(l-CFLOSc) 
C=1 


The  probability  of  seeing  either  clouds  or  the  ground  (with  or  without 
clouds) 

159 

=  WT(  I )  +  WTold  (160)  +  OCOVER(l) 

1  =  1 

4 

=  CC0VER(C+1 )  [1  -  CFL0Sc  +  CFL0Xc]  +  CC0VER(1)=1  . 

C=1 

3-2.9  Subroutine  CLOUDO(MODE) 

For  a  detailed  description  of  CLOUDO  see  Section  4-2. 

3-2.10  Subroutine  CLOUD 1 

Coordinates  of  the  source  (XI ,ETA,ZETA)  and  detector  (XID, 
ETAD,ZETAD)  are  calculated  in  the  Cartesian  coordinate  system  of  each 
cloud  with  the  quantities  stored  in  common  CDCORD.  Entry  CLD  E  TO  C 
in  subroutine  CL0UD9  is  used  to  transform  from  earth  to  cloud  coordinates. 

3-2.11  Subroutine  CL0UD2(UL ,VL,WL) 

The  altitude  (ALTLK(K))  and  geographic  coordinates  (THETALK(K) , 
PHILK(K))  of  the  intersection  of  the  line-of -sight  and  the  Kth  deter¬ 
ministic  cloud  are  calculated  in  this  routine  and  outputted  in  common 
INTRCPT.  These  results  are  not  used  by  other  cloud-model  routines. 
However,  this  routine  is  a  convenience  to  the  user.  The  sight-path 
integration  module  can  use  results  of  this  routine  to  avoid  calling  the 
cloud  transfer-coefficient  routine  (CL0UD3)  until  the  sight-path 
integration  has  reached  the  cloud  top. 

The  input  to  this  subroutine  are  the  direction  cosines 
(UL,VL,WL)  of  the  line-of-sight .  The  point  of  intersection  is  calculated 
for  all  clouds  which  fall  in  the  path.  A  flow  chart  is  shown  in 
Figure  3-3. 
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3-2.12  Subroutine  CL0UD3(UL,VL,WL,K,ALAM,TC0EF,EMISS,SINCLD) 

When  all  of  the  bookkeeping  has  been  accomplished,  CL0UD3 
is  called  to  calculate  the  transfer  coefficient  (TCOEF)  and  the  emission 
(EMISS).  Branching  is  required  depending  upon  whether  or  not  the  point 
of  intersection  of  the  LOS  with  the  cloud  is  visible  to  the  source. 

For  those  LOS-intersection  points  which  are  also  illuminated,  the 
bidirectional  reflectance  is  obtained  from  tables  in  function  CLDBDR. 

The  tabulated  values  are  functions  of  incident  and  exit  zenith  angles 
and  the  exit  azimuthal  angle.  The  bidirectional  reflectance  is  picked 
from  a  bin  with  no  interpolation.  The  tables  are,  however,  interpolated 
linearly  as  a  function  of  wavelength. 

The  irradiance  to  all  illuminated  facets  is  evaluated.  A 
diffusion  calculation  is  then  performed  in  subroutine  DIFFUS  to  get  the 
radiance  exiting  the  facet  of  interest. 

The  thermal  emission  is  calculated  from  the  Planck  distri¬ 
bution  for  the  temperature  at  the  point  of  intersection  of  the  cloud 
top  and  line-of-sight . 

Before  the  source  irradiation  of  a  cloud  is  computed,  a  check 
is  made  to  determine  whether  or  not  the  source  is  either  inside  the 
cloud  or  on  its  surface.  If  the  source  is  not  outside  the  cloud,  a 
warning  message  is  printed  and,  without  further  calculation,  the 
routine  returns  a  logical  variable  (SINCLD  =  .TRUE.)  to  the  main  pro¬ 
gram  (ROSCOER)  which  then  proceeds  to  the  next  case. 

We  have  also  devised  a  prescription  which  is  used  if  the 
artificial  (point)  source  is  closer  to  the  cloud  surface  than 
(approximately)  the  radius  Rg  that  is  assigned  to  the  finite-size  arti¬ 
ficial  source.  Let  R?  be  the  position  of  the  source  in  the  cloud 
coordinate  system  and  RS  be  the  radial  distance  to  the  intersection 
point  S  of  the  ray  RP  with  the  cloud  surface.  At  point  S  the  outward 
unit  normal  is  ft(R§).  Then,  for  the  purpose  of  this  routine,  if 

(RP-R§)  •  n(R§)  <  Rg  , 

the  source  is  moved  in  the  direction  of  n(R§)  an  amount  <SP  given  by 
6P  =  0.5  [Rs  -  (R?-r3)  •  n(R§)]  . 
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A  message  is  printed  with  the  original  and  altered  coordinates  of  the 
source . 

3-2.13  Subroutine  CL0UD9(RTP1 , RTP2 , RTP3 , EC1 , EC2 , EC3 , K) 

Transformations  are  made  between  the  earth-centered  polar 
coordinate  system  and  the  cloud-centered  Cartesian  system  of  the  Kth 
cloud.  This  routine  has  two  entry  points.  Entries  CLD  E  TO  C  and  CLD  C 
TO  E  transform  from  earth  to  cloud  and  cloud  to  earth,  respectively. 

The  cloud  geometry  is  most  easily  described  in  the  cloud  .  .  . 
principal-axis  system  (C,n,S)  in  which  the  equation  of  the  cloud  surface 
is 


The  location  of  the  cloud  center  (earth-centered  polar  coordinates)  is 
(Rc,ec,<t’c)  where 

R  =  R  .  +  H  ,  H  =  HCLOUD  (km) 
c  earth  c’  c 

@c  =  THETCD  (degrees) 
ij>c  =  PHICD  (degrees). 

Looking  from  above  the  cloud  toward  the  earth  center: 


The  transformation  from  the  earth-centered  Cartesian  system 
(x,y,z)  to  the  non-rotated  (ipc=0)  Cartesian  system  (X,Y,Z)  in  the  cloud 
can  be  performed  with  two  rotations  and  one  translation. 

The  first  intermediate  coordinate  system  (x',y',z')  is 
obtained  by  making  a  counterclockwise  (positive)  rotation 

<i>  =  <f>c  +  */2 


about  the  z-axis  such  that  the  x'-axis  is  normal  to  the  vector  S  to 

c 

the  cloud  center.  The  (ir/2)-term  is  introduced  to  make  the  positive 
X-axis  of  the  cloud  point  east.  The  transformation  is 


=  RjC*) 


where  the  transformation  matrix  R^  is  given  by 


cos 

RjC*)))  =  -  sin 


V  ° 


The  second  intermediate  coordinate  system  (z",y",z")  is 

obtained  by  making  a  counterclockwise  (positive)  rotation  0c  about  the 

x'-axis  such  that  the  z'-axis  lies  along  the  vector  R  .  The  trans- 

c 

formation  is 


y" J  =  R2(0)  I  y' 


where 
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Finally,  the  non-rotated  cloud  system  (X,Y,Z)  is  related  to 
the  ( x" , y" , z" )  system  by  a  translation  along  the  z"-axis: 


Z  =  z"  -  R„ 
c 

where 

Rc  =  Rearth  +  Hc 
Collating  our  results  we  have 


=  R2(0  )R1(«(») 


(2a) 


with 


R2(0)R1(^) 


sin 

*c 

cos 

*c 

0 

cos 

6 

c 

cos 

e 

c 

-  cos 

e 

C 

sin 

$ 

c 

sin 

sin 

e 

c 

cos 

<f> 

c 

sin 

0 

C 

sin 

<|> 

c 

cos 

(2b) 


Thus,  an  arbitrary  point  in  earth-centered  polar  coordinates 
(R,  ©,$)  has  earth-centered  Cartesian  coordinates  (x,y,z)  and  cloud 
principal-axis  coordinates  (5,n,C)  where 


( x ,  y ,  z)  =  (R  sin  0  cos  4> ,  R  sin  0  sin  $,  R  cos  0) 

and  the  coordinates  (X,Y,Z)  and  (x,y,z)  are  related  by  Equations  (2a) 
and  (2b).  Finally  (£,n,?)  and  (X,Y,Z)  are  related  by  Equation  (1),  or 
equivalently 

=  Ri^c> 

The  transformation  from  the  cloud  system  to  the  earth  system 
may  be  done  by  doing  the  same  steps  in  reverse  order  and  using  the 
transposed  transformation  matrices. 
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3-2.14 


Subroutine  CONVERT( R , TH , PH , X , Y , Z ) 

This  subroutine  has  entry  points  to  transform  from  polar  to 
Cartesian  and  Cartesian  to  polar  coordinates.  The  entry  points  are 
POL  TO  CT  and  CT  TO  POL  with  arguments  as  follows: 

Y  -  Cartesian  X-coordinate 

Y  -  Cartesian  Y-coordinate 
Z  -  Cartesian  Z-coordinate 

R  -  Radius  in  spherical  coordinate  system 
TH  -  Polar  angle  in  spherical  coordinate  system 
PH  -  Azimuthal  angle  in  spherical  coordinate  system. 

3-2.15  Subroutine  DIFFUS(K, ALAM) 

The  computations  involved  in  subroutine  DIFFUS  are  those  from 
Appendix  A  and  as  described  in  Section  2-1 .  Communication  with  this 
routine  is  from  CL0UD3 . 

3-2.16  Subroutine  DIFPRM(ALAM) 

Given  a  wavelength  between  two  and  five  microns,  this  sub¬ 
routine  interpolates  microscopic  scattering  and  extinction  cross- 

2  -1 
sections  in  urn  ,  coverts  to  macroscopic  cross-sections  in  units  of  km 

and  records  the  modified  diffusion  parameters,  as  defined  in  Section 

A-1.2,  for  each  of  the  14  cloud  types  of  the  data  base. 

INPUT: 

ALAM:  wavelength  of  photon  (pm) 

OUTPUT :  (common  DIFFUS,  arrays  of  14-word  length) 

DALFA  =  a  in  Appendix  A  Equation  (2)  (dimensionless) 

DBETA  =  b  in  Appendix  A  Equation  (2)  (km) 

DKAPPA  i  K  in  Appendix  A  Equation  (19)  (km-1) 

For  large-particle-size  clouds  (i.e.,  cloud-types  11,12,13 

_2 

and  14),  the  small  wavelength  approximations  =  2vr  and  oc  = 

_ 2  *  ^ 

(1+W)ur  are  used.  Here,  1  and  W  are,  respectively,  the  diffracted 

and  specular-ref lected  components  of  the  scattering. 
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3-2.17  Function  DSURF( THETA, PHI ,K) 

Differential  area  elements  on  the  surface  of  an  ellipsoid  are 
calculated  in  this  routine  as  a  function  of  0  and  <p .  The  cloud  number 
K  must  also  be  transmitted  to  the  function  as  it  obtains  the  ellipsoidal 
parameters  from  common  CDCORD. 

3-2.18  Function  EM I SSF ( ALAM , LL , MM , HH ) 

2 

This  function  provides  the  thermal  emission  (W/(km  sr  pm)) 
from  clouds,  attenuated  by  clear  air  to  an  altitude  of  12  km. 

3-2.19  Subrout ine  ESURF ( THI , THR , PS I , ZKM , MSM , DD , SPCULR , ZLAM , IDAY , 
IFIRES , ESURFI , SFR , EPSD , TKS ) 

This  subroutine,  a  dummy  routine  for  the  real  one,  provides 
the  surface  upwelling  radiation  and  earth  bidirectional  reflectance. 

The  description  is  given  in  Reference  6. 

3-2.20  Function  IPRIME(II) 

This  function  provides  a  mapping  of  the  facet  areas  and  radii 
onto  those  of  the  first  octant.  This  symmetry  is  utilized  to  obviate 
storing  redundant  data. 

3-2.21  Subroutine  LINEAR(XD , FXD , XX , FXX , NX) 

LINEAR  performs  a  linear  interpolation  with  an  expeditious 
technique  to  find  the  pair  of  tabular  points  which  bound  the  point  of 
interest.  The  interpolation  is  then 


Y 


Y 


i-1 


+ 


Y.-Y 

l 

X.-X 

l 


i-1 

i-1 


(X-X._ 


) 


INPUTS: 

XX:  array  of  points  in  independent  variable 

FXX:  array  of  points  in  dependent  variable 

NX:  number  of  points  in  the  table 

X0:  value  of  independent  variable  for  which  dependent 

variable  is  desired. 

OUTPUTS : 

FXO:  dependent  variable  at  value  of  X0  in  independent 

variable . 
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3-2.22 


Subroutine  LOS(K) 

Possible  lines-of-sight  from  which  the  user  can  select 
directions  are  computed  and  stored  in  common  LOS.  There  are  selection 
rules  based  on  several  criteria.  A  line-of-sight  must  obviously  inter¬ 
cept  a  cloud  facet;  however,  the  facet  is  not  necessarily  illuminated. 
Facets  which  would  fall  within  the  same  f ield-of-view  are  not  dupli¬ 
cated.  See  Figure  3-4  for  a  flow  chart. 

The  first  criterion  selects  first  the  facets  which  are  both 
visible  and  illuminated  and  orders  these  from  smallest  to  largest 
scattering  angle.  Those  facets  which  are  visible  but  not  illuminated 
are  then  ordered  according  to  their  being  closest  to  the  edge  of  the 
cloud.  The  rationale  for  this  ordering  is  that  the  cloud  will  be 
thinnest  along  the  line-of-sight  near  the  edge. 

This  routine  can  produce  a  maximum  of  16  possible  lines-of- 
sight.  For  a  small  cloud  and  large  f ield-of-view ,  there  could  possibly 
be  only  one. 

3-2.23  Function  PS ( XI , ET1 , ZET1 , THETA , PHI , K ) 

This  function  evaluates  a  quantity  which  is  proportional  to 
the  dot  product  of  the  normal  to  the  ellipsoid  and  a  vector  from  the 
surface  point  to  the  external  point.  If  the  quantity  is  negative,  the 
point  is  shadowed  from  the  line-of-sight;  if  positive,  the  point  is 
visible . 

3-2.24  Function  RHOEPS ( ALAM , CTHETA ) 

This  function  evaluates  the  hemispherical-directional 
reflectance  from  the  (horizontal)  surface  of  cloud-type  4  (strato- 
cumulus)  for  radiation  of  wavelength  ALAM  reflected  into  zenith  angle 
cos-1  (CTHETA). 

Tabulated  in  an  array  called  RHO  is  a  set  of  hemispherical- 
directional  reflectances  for  the  stratocumulus  cloud  model  at  wave¬ 
lengths  between  2  and  5  microns  and  for  exit  zenith  angles  whose 
cosines  are  1.0,  0.866,  0.5,  and  0.0.  The  value  of  0.0  is  a  fictitious 
endpoint  to  the  table  to  accommodate  interpolation.  Interpolation  is 
linear  in  wavelength  and  cosine  theta. 

INPUT: 

ALAM:  wavelength  of  photon  (pm) 

CTHETA:  cosine  of  zenith  angle  for  which  emissivity 

is  desired. 
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Figure  3-4.  Flow  chart  of  subroutine  LOS. 
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Figure  3- 
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Flow  chart  of  subroutine  LOS  (Cont. ) 
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OUTPUT : 


3-2.25 


RHOEPS:  hemispherical-directional  reflectance 


Subroutine  SCLOUD(ALAM,U,V,W, ISUN) 


Subroutine  SCLOUD,  for  which  a  flow  chart  is  given  in  Figure 
3-5,  calls  in  the  required  subroutines  to  compute  transfer  coefficients 
and  thermal  emission  for  statistical  clouds.  A  source  position  must  be 
specified,  just  as  for  deterministic  clouds.  It  is  also  required  that 
the  solar  longitude  and  latitude  be  specified  when  the  solar  position  is 
above  the  local  horizon  at  the  intersection  of  the  line-of-sight  with 
the  earth.  Calculations  are  made  for  the  source  (burst)  and  solar 
transfer  coefficients  and  the  thermal  emission.  Calculations  are  by¬ 
passed  for  the  solar  transfer  coefficient  when  the  sun  is  below  the 
local  horizon. 

By  interrupting  this  subroutine,  one  has  access  to  the  entire 
cumulative  arrays  (TRANS(I),  EMISS(I),  I  =  1,160)  in  addition  to  the  four 
fractiles  returned  to  the  calling  program.  This  could  easily  become  a 
normal  mode  of  operation  for  system-type  studies.  Also,  this  (non- 
implemented)  mode  is  appropriate  if  one  desires  calculations  for  a 
specific  realization  of  the  statistical  cloud  configurations  without 
going  to  the  deterministic  mode  of  operating  the  code. 

Subroutine  SCLOUD  calls  subroutine  CLDWT  to  obtain  TRANS(I), 
EMISS(I),  and  WT( I )  for  1=1,159.  For  1=160,  when  the  LOS  is  cloud-free 
and  terminates  on  the  ground,  these  variables  are  evaluated  in  subroutine 
SCLOUD.  In  computing  the  transfer  coefficient  for  1=160,  one  could  get 
arbitrarily  large  values  if  he  let  the  (point)  source  get  arbitrarily 
close  to  the  ground  (cf.  the  equation  for  the  reflectance  transfer 
coefficient  for  the  statistical  cloud  submodel  near  the  end  of  Section  2). 
To  prevent  getting  such  unrealistic  values  for  the  transfer  coefficient, 
we  have  imposed  the  requirement  that  the  main  program  assign  a  radius  R 

s 

to  the  finite-size  artificial  source.  Then,  if  the  altitude  h  of  the 

s 

source  is  less  than  R  ,  we  temporarily  let  the  effective  altitude  of 

s 

the  finite-size  source  be 

heff  =  °-5  <Rs  +  hs>  • 

Thus,  if  the  source  is  initially  placed  on  the  ground  so  that  hg  =  0,  our 

prescription  places  the  effective  altitude  at  0.5  R  . 

s 


. i"1  .  “'T 


Figure  3-5.  Flow  chart  of  subroutine  SCLOUD. 
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INPUTS : 


ALAM:  wavelength  of  incident  photon  (pm) 

U,V,W:  direction  cosines  of  line-of-sight  in  the  earth- 

centered  system 

ISUN:  A  flag  to  denote  nature  of  the  source,  being  zero 

for  an  artificial  source  and  unity  for  the  sun. 

See  Section  3-3  for  definitions  of  variables  in  the  following  name 
commons : 


Section  No. 

Section  No 

CLDFREQ 

3-3.4 

PARAMS 

3-3.15 

CLDWT 

3-3.5 

PTE 

3-3.16 

DETECT 

3-3.8 

SANDD 

3-3.17 

FLAGS 

3  '  10 

SOLARP 

3-3.18 

MATERL 

3-  .  j.4 

SORCE 

3-3.19 

Common  PTE  is  used  only  in  SCLOUD  and  the  calling  program.  The  variable 
definitions  are 

( input ) 

P1,P2,P3,P4  -  four  probability  levels  (fractiles) 
for  which  transfer  coefficients  and 
thermal  emission  are  required. 

(output) 

T1,T2,T3,T4  -  four  transfer  coefficients  for  a 
NUDET  source  corresponding  to  the 
P1,P2,P3,P4  above 

E1,E2,E3,E4  -  thermal  emission  radiance  corresponding 
to  P1,P2,P3,P4  above 

S1,S2,S3,S4  -  four  transfer  coefficients  for  the 
solar  source  corresponding  to  PI , 

P2,P3,P4  above 

PCFLOS  -  probability  of  the  detector  seeing  the 

artificial  source,  either  with  or  without 
clouds  being  present,  provided  the  source 
is  within  the  f ield-of-view. 
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3-2.26 


Subroutine  SGEOM(UL , VL , WL , CLDTOP , XD , YD , ZD , XS , YS , ZS , RDIST , 
SDIST , THETI , THETE , PHIE ) 

Subroutine  SGEOM,  for  which  a  flow  chart  is  given  in  Figure 
3-6,  computes  the  incident  and  exit  zenith  angles  and  the  azimuthal 
scattering  angle  for  a  specified  source  and  detector  geometry  and  a 
given  line-of-sight .  Also  evaluated  is  the  distance  (SDIST)  from  the 
source  to  the  cloud-LOS  intersection  point  anl  the  horizontal  component 
thereof  (RDIST). 


INPUTS : 

UL,VL,WL:  direction  cosines  of  the  line-of-sight  in 

the  earth-centered  coordinate  system 

CLDTOP:  altitude  of  top  of  cloud  or  other  specified 

altitude  (km) 


XD,YD,ZD:  rectangular  components  of  detector  position 

vector  in  earth-centered  system  (km) 


XS,YS,ZS:  rectangular  components  of  source  position 

vector  in  earth-centered  system  (km) 


OUTPUTS : 


SDIST:  distance  from  source  to  intersection  of 

line-of-sight  with  cloud  top  (km) 

RDIST:  horizontal  component  of  SDIST  (km) 

THETI :  zenith  angle  to  source  position  at 

intersection  point  of  LOS  with  CLDTOP- 
altitude  (degrees) 

THETE:  zenith  angle  to  detector  position  (degrees) 

PHIE:  azimuthal  angle  of  scattering  of  photon 

(degrees) 


3-2.27  Subroutine  SORTLJ ( A , B , C , N , LOHI ) 

This  subroutine  sorts  an  array  in  either  ascending  or  des¬ 
cending  order  and  carries  two  additional  arrays  with  the  ordering.  The 
argument-list  variables  are  as  follows: 

A  is  the  array  to  be  sorted 
B  array  carried  along  with  the  sorting  on  A 
C  array  carried  along  with  the  sorting  on  A 
N  is  the  number  of  elements  in  arrays  A,B,  and  C 
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Figure  3- 
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in 


-1 


(RS-RI ) *RI 
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6.  Flow  chart  of  subroutine  SGEOM 
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LOHI  is  the  ordering  flag 

>0  sorts  from  high  to  low 
£0  sorts  from  low  to  high. 

3-2.28  Function  TA I R ( ALAM , THETA , ALT ) 

Air  transmittance,  T,  for  radiation  of  wavelength  X  from 
altitude  h  to  12  km  along  a  path  at  zenith  angle  0  was  calculated 
using  LOWTRAN  3  for  the  U.S.  1962  Standard  Atmosphere.  The  data  are 
cataloged  in  this  function. 

The  independent  variables  are: 
wavelength:  2  <  X  <5  microns 

zenith  angle:  0  £  0  £85  degrees 
altitude:  1  <  h  <12  kilometers. 


The  interpolation  (and  extrapolation  for  0  >85  degrees  and  h  <1  km) 
used  is  linear  in  (T,X)-,  (In  T,  sec  0)-,  and  (In  T,h)-space.  More 
explicit  formulas  follow: 


Altitude 


T 

b 

a 

which  implies 


In  T  =  In  T1  -  (In  Tj 


in  V 


h-h 

h2~h 


Angle 


T 


ab 


sec© 


b  =  (Tg/T^ 


(sec  e2-sec  9^ 


-1 


a  =  T. /bsec6l  , 
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which  implies 


In  T  =  In  Tj  +  (In  T2 


In  Tj) 


secS-sece^ 
sec0  2~sec6  ^ 


For  0  >85  degrees,  we  redefine  sec0  by  a  parabolic  approximation  to  the 
Chapman  function, 

sec0  =  sec  (85  deg)  +1.1  (0-85)2 


Wavelength 


X-Xn 

T  =  T1  +  (T2-ti> 


3-2.29  Function  TRANSF ( ALAM , LL , MM , HH ) 

This  function  evaluates  the  transfer  coefficient  for  the  set 
of  statistical  clouds  whose  indexes  are  in  the  calling  sequence:  Con¬ 
siderable  bookkeeping  is  accomplished  to  assign  the  calculation  to  the 
proper  prescription.  The  prescriptions  are  adequately  documented  in 
Appendix  A  and  discussed  in  Section  2-2. 

The  bookkeeping,  in  general,  evaluates  the  diffusion  parameters 
for  the  appropriate  layer  number  and  stores  them  according  to  the 
required  variable  name  for  the  prescription.  The  first  steps  are  to 
determine  how  many  layers  are  involved  in  the  problem,  whether  the 
source  is  inside  or  outside,  and  whether  the  cloud  layers  are  continuous. 
These  questions  are  answered  and  control  transferred  to  the  equation  for 

reflection , 

1-layer  transmission  for  source  inside, 

1- layer  transmission  for  source  below, 

2- layer  transmission  for  source  below, 
or  3-layer  transmission  for  source  below. 

The  results  are  then  normalized  to  get  a  transfer  coefficient,  TRANSF, 
which  is  returned  to  subroutine  CLDWT. 

If  the  altitude  ALT  of  an  artificial  source  above  the  top  of 
the  topmost  layer  is  less  than  the  radius  R  assigned  to  the  finite-size 
source,  the  effective  altitude  of  the  source  in  computing  the  reflectance 
transfer  coefficient  is  reset  to 


heff  =  CLDTOP  +0.5  (ALT  +  Rg)  . 

3-2.30  Subroutine  VSN0RM(X , Y , Z , EN1 , EN2 , EN3 , K) 

This  subroutine  computes  a  unit  vector  n  normal  to  the  sur¬ 
face  of  the  ellipsoid  at  the  Cartesian  coordinates  x,y,z  in  cloud  K. 


3-3  LABELED  COMMON  VARIABLE  DEFINITIONS 

Variables  in  labeled  common  are  described  in  this  section. 
Since  the  subroutines  which  use  these  variables  are  designed  to  be 
utilized  as  part  of  other  codes,  i.e.,  the  subroutines  are  submodules, 
names  are  located  in  labeled  commons  which  will  tend  to  minimize  their 
conflict  with  other  code  variable  names.  A  list  of  the  names  of 
labeled  commons  and  the  subroutines  in  which  they  are  found  is  located 
in  Table  3-1 . 

The  name  of  each  common  and  the  variables  which  are  con¬ 
tained  therein  will  be  listed  just  as  they  are  found  in  the  subroutines. 
A  brief  description  and  the  variable  units  follow  the  statement. 

3-3.1  COMMON/ ATM0UP/DUM1( 4 ) , RHO , TT ,DUM2( 22) 

3 

RHO:  RH0(g/cm  ),  the  air  density  at  altitude  z(km), 

is  used  by  the  main  program  ROSCOER  to  scale 
the  radius  of  the  artificial  source  with  altitude. 

TT:  TT  (degrees  Kelvin),  the  air  temperature  at 

altitude  z(km),  is  used  by  the  cloud  routines. 

(RHO  and  TT,  the  only  two  variables  from  this 
common,  are  obtained  by  a  call  to  ATMOSU ( 2 , ZH ) , 
where  ZH  is  the  altitude  of  interest,  (km)). 
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Table  3-1.  List  of  subroutines  and  labeled  commons. 


3-3.2  C0MM0N/CDC0RD/XI(1 ,20) ,ETA(1 ,20) , ZETA( 1,20) ,XID(20 ) , ETAD( 20) , 
ZETAD( 20 ) , ASQ( 20 ) , BSQ( 20) , CSQ( 20) , CPH( 5 , 20) , CTH( 5 , 20) , CPHBAR 
(16,20), CTHBAR( 8,20), CRBAR( 4,4,20), CAREA( 4,4,20) 

Variables  in  this  common  provide  communication  between  the 
geometry  routines  of  the  deterministic  cloud  calculations.  Calculations 
relative  to  a  given  cloud  are  made  in  a  coordinate  system  centered  on 
that  cloud  and  results,  if  required,  are  transformed  back  to  the  earth- 
centered  system.  Ellipsoidal  elemental  areas  used  in  integration  are 
referred  to  as  facets. 

XI,ETA,ZETA:  Cartesian  coordinates  of  the  Nth  source  in 

the  Kth  cloud-centered  system.  N  currently 
is  restricted  to  one.  (km) 

XID , ETAD , ZETAD :  Cartesian  coordinates  of  the  detector  in  the 
Kth  cloud-centered  system,  (kin) 


ASQ,BSQ,CSQ:  Values  of  the  squares  of  the  semi-axes  of 

an  ellipsoid,  (km^) 


CPH:  Cloud-centered  azimuthal  angle  bounding  a 

facet  or  area  on  the  surface  of  the  cloud, 
(radians) 

CTH:  Cloud-centered  polar  angle  bounding  a  facet 

on  the  surface  of  the  cloud,  (radians) 

CPHIBAR:  Cloud-centered  azimuthal  angle  locating  a 

facet  center  in  the  Kth  cloud.  Symmetry  is 
taken  into  account  in  the  eight  octants, 
(degrees) 

CTHBAR:  Cloud-centered  polar  angle  locating  a  facet 

center  in  the  Kth  cloud,  (degrees) 

CRBAR:  Distance  from  the  cloud  center  to  the  facet 

center,  utilized  in  calculating  fluxes.  All 
symmetry  is  taken  into  account  and  the  function 
IPRIME  locates  the  element  in  the  CRBAR  array 
which  has  the  proper  distance  to  one  of  the 
128  facets  when  there  are  16  azimuthal  angles 
and  8  polar  angles,  (km) 

CAREA:  Areas  of  the  facets  are  calculated  by  a 

2-dimensional  Simpson's  rule  of  integration 
in  subroutine  CLDGEOM  with  the  differential 
area  element  being  calculated  in  DSURF.  Again, 
function  IPRIME  returns  the  proper  index  to 
obtain  the  desired  area  element,  (km2) 
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3-3.3  COMMON/CDDATA/ALAMT( 10) ,SIGT (10,14) ,SIGS(10,14) ,W(10)  , 

XMUBAR( 10,14) ,DN0( 14 ) , RAD( 14 ) , CLDBASE( 14 ) , CLDTHK( 14 ) 

Data  in  this  common  are  physical  and  optical  properties  of 
the  cloud  data  base.  Physical  properties  are  taken  from  the  NASA  cloud 
data  base  and  optical  properties  for  cloud  types  1  through  10  are  from 
Mie  theory  calculations.  There  are  14  types  of  clouds  in  this  data  base 
for  which  data  are  tabulated  at  ten  wavelengths. 

ALAMT:  ALAMT  is  a  table  of  wavelengths  for  which  micro¬ 

scopic  cross-sections  and  mean  cosines  of 
scattering  angles  are  tabulated.  Values  are 
2.0,  2.5,  2.55,  2.65,  2.75,  2.85,  3.2,  3.5, 

4.0,  and  5.0.  (ym) 

SIGT:  Microscopic  extinction  cross-sections  for  the 

14  cloud  types  are  tabulated  at  the  10  wave¬ 
lengths  of  the  ALAMT  table,  (ym2) 

SIGS:  Scattering  cross-sections  for  the  14  cloud 

types  are  tabulated  at  the  10  wavelengths  of 
the  ALAMT  table,  (ym2) 

W:  (Normalized)  contribution  of  specularly-reflected 

radiation  to  scattered  radiation,  with  the 
diffracted  radiation  contributing  unity;  used 
only  for  the  large-particle-size  clouds,  types 
11  through  14.  (dimensionless) 

XMUBAR:  Cosine  of  the  scattering  angle  averaged  over 

the  particle  size  distribution  of  the  cloud 
type  and  tabulated  at  the  10  wavelengths, 
(dimensionless) 

DNO:  Droplet  density  of  the  given  cloud  type. 

(droplets/cm2) 

RAD:  Mean  droplet  radius  for  the  cloud  type,  (ym) 

CLDBASE:  Cloud-base  altitudes  for  the  14  cloud  types,  (km) 

CLDTHK:  Cloud  thickness,  (km) 

3-3.4  COMMON / CLDFREQ / KMODEL , CCOVER (5,11), CLDFREQ (17,4 ,11) 

Common  CLDFREQ  has  its  origin  in  block  data.  The  index 
KMODEL  selects  one  of  the  ten  sets  of  data  present  in  data  statements 
or  allows  a  user-supplied  set  of  data  to  be  input  as  the  eleventh  set 
of  data. 
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KMODEL:  A  variable  which  selects  the  cloud  data  base 

to  be  used  in  the  calculations. 

=1  for  dry  annual  average 
=2  for  dry  winter  average  night 
=3  for  dry  winter  average  day 
=4  for  dry  summer  average  night 
=5  for  dry  summer  average  day 
=6  for  moist  annual  average 
=7  for  moist  winter  average  night 
=8  for  moist  winter  average  day 
=9  for  moist  summer  average  night 
=10  for  moist  summer  average  day 
=11  for  array  supplied  by  the  user. 


See  Section  4  for  a  detailed  illustration  of  using 
KMODEL  =11. 

CCOVER:  This  array  provides  the  probability  of  the 

Lth  cloud-coverage  category  for  the  Kth 
model,  K  being  defined  above,  where 


5 

^  CCOVER ( L , K )  =  1 
L=1 


For  the  K  =  11  base,  the  input  instructions  in  Section  4 
illustrate  the  use  of  the  table.  CC0VER(L,K)  is  a  frequency  of 
occurrence  of  the  Lth  coverage  category. 

CLDFREQ :  With  the  indexes  (I,L,K),  CLDFREQ( I , L , K)  is 

the  frequency  of  the  lth  cloud  type  for  I 
less  than  fifteen.  For  I  =  15,16,17  it  is 
the  frequency  of  occurrence  of  1-layer, 
2-layers,  or  3-layers,  respectively.  For 
K-values  between  1  and  10  (inclusive)  the 
array  elements  are  from  the  model  as  defined 
under  KMODEL  above.  For  K  =  11 ,  the  use 
of  the  table  is  described  in  Section  4. 
Basically,  however,  each  element  is  the 
frequency  of  the  lth  type  of  cloud: 


1=1  for  Cul 
1=2  for  Cu2 
1=3  for  Cu3 
1=4  for  Sc 
1=5  for  St 
1=6  for  Cbl 
1=7  for  Ns 


1=8  for  Ac 
1=9  for  Asl 
1=10  for  Cb2 
1=11  for  Cb3 
1=12  for  As2 
1=13  for  Ci 
1=14  for  Cs 
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3-3.5  COMMON/CLDWT/ IDX ,  WT( 160) , TRANS ( 160 ) , EMISS( 160) 

The  common  array  WT  contains  the  probability  of  a  given 
statistical  cloud  configuration,  i.e.,  the  product  of  the  probability  of 
a  given  cloud-coverage  category,  the  probability  of  the  number  of  layers 
in  that  coverage  category,  and  the  probability  of  an  occurrence  of  a 
cloud  type  in  that  layer  and  category  configuration.  This  configuration 
is  ultimately  weighted  by  one  minus  the  cloud-free  line-of-sight  for 
this  cloud  coverage  and  zenith  angle. 

There  are  160  combinations  of  layer-type-coverage  configurations. 
Subsets  of  restrictions  on  any  of  the  components,  e.g.,  coverage  cate¬ 
gory,  reduces  the  160  to  a  fewer  number  of  cases. 

IDX:  This  is  the  total  number  of  configurations 

which  is  always  less  than  160.  With  a  maximum 
constraint,  it  can  be  as  few  as  two,  i.e.  the 
cloud  type  (with  coverage  category  specified) 
and  the  cloud-free  line-of-sight  component. 

WT:  WT(IDX)  is  the  probability  of  the  IDXth  con¬ 

figuration.  IDX  is  less  than  or  equal  to 
160. 

TRANS:  TRANS(IDX)  is  the  transfer  coefficient  for 

the  IDXth  configuration.  It  corresponds  to 
either  a  transmission  or  a  reflection  and 
includes  attenuation  by  the  atmosphere  (1962 
U.S.  standard  via  LOWTRAN-3)  to  an  altitude 
of  12  km.  (sr-l) 

EMISS:  EMISS(IDX)  is  the  thermal  emission  spectral 

radiance  attenuated  to  12-km  altitude. 

(W/(km2  sr  ym)) 

3-3.6  COMMON / CLOUD / NCLOUD , HCLOUD (20), THETCD (20) ,PHICD(20) , 

CL0UDA( 20 ) , CL0UDB( 20) ,CLOUDC(20) ,CLDPSI(20) ,KCLOUD(20) 

Input  deterministic  cloud  descriptions  are  read  into  these 
arrays  by  subroutine  CLOUDO.  As  many  as  twenty  clouds  can  be  part  of 
the  calculations;  however,  there  are  no  interacting  clouds,  i.e., 
reflections  of  one  cloud  to  another  are  not  considered. 

NCLOUD:  Number  of  deterministic  clouds  for  which 

altitude,  orientation,  and  semiaxes  are  to 
be  read  into  the  arrays.  This  variable  is 
used  as  the  upper  limit  on  DO-loops  throughout 
the  code. 

THETCD:  Colatitude  of  the  cloud,  (degrees) 
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PHICD:  East  longitude  of  the  cloud,  (degrees) 

CLOUDA:  Semiaxis  of  the  longer  cloud  axis  in  horizontal 

plane,  (km) 

CLOUDB:  Semiaxis  of  the  shorter  cloud  axis  in  the 

horizontal  plane,  (km) 

CLOUDC:  Semiaxis  of  the  cloud  in  vertical  direction, 

(km) 

The  values  of  a,b,  and  c  are  those  in  the  equation  of  an 

ellipsoid : 

(if  ♦  (§)2  *  (if  ■  * 


CLDPSI:  Azimuthal  angle,  measured  counterclockwise 

looking  down,  orienting  cloud  about  (vertical) 
z-axis.  Clouds  are  normally  oriented  with 
(+a)-axis  eastward  (psi  =  0).  (degrees) 

KCLOUD:  Index  for  Kth  cloud,  representing  one  of  the 

14  cloud  types  from  the  data  bank. 

3-3.7  COMMON / CONF I G/ CTYPL ( 8 ) , CTYPM ( 4 ) , CTYPH ( 5 ) , CFLAG , Cl , C2 

This  common  is  a  bit  unwieldy  but  is  used  in  the  input 
routine  CLOUDO .  The  entries  are  the  Hollerith  names  of  the  cloud 
types  in  the  low,  middle,  or  high  layer. 

CTYPL ( LL )  :  Type  of  cloud  LL  (1  £LL  £7)  or  no  cloud 

(LL=8 )  in  the  low  layer.  The  variable  is 
lef t-justif ied  in  5-column  fields. 

CTYPM(MM):  Type  of  cloud  MM  (1  <MM  £3)  or  no  cloud 

(MM=4)  in  the  middle  layer. 

CTYPH (HH) :  Type  of  cloud  HH  (1  <HH  £4)  or  no  cloud 
( HH=5 )  in  the  high  layer. 

CFLAG:  CFLAG  is  a  coverage  category  flag.  (Not 

utilized) 

Cl:  Cl  is  the  (integer)  index  for  the  first 

cloud-coverage  category  being  utilized. 
Normally  Cl  is  2. 

C2:  C2  is  the  upper  limit  counterpart  of  Cl. 

Normally  C2  is  5. 
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3-3 . 8  COMMON /DETECT/ HD , THETAD , PHID , FOV 

This  common  contains  the  location  of  the  detector  in  an  earth- 
centered  polar  coordinate  system. 

HD:  Altitude  of  the  detector,  (km) 

THETAD:  Colatitude  of  detector,  (degrees) 

PHID:  East  longitude  of  detector,  (degrees) 

FOV:  Field-of-view  (FOV)  of  the  detector  (km*2). 

Radiances  in  the  deterministic  cloud  calculations 
are  averaged  over  the  field-of-view.  The  FOV 
is  also  used  in  subroutine  LOS,  as  potential 
lines-of-sight  are  separated  by  the  order  of  the 
square  root  of  the  FOV. 

3-3.9  COMMON/DIFFUS/DALFA( 14 ) , DBETA( 14 ) ,DKAPPA(14) 

This  common  contains  the  coefficients  in  the  diffusion  equation 
for  photon  flux  transmitted  through  clouds.  They  are  generated  for  all 
fourteen  clouds  by  a  single  call  to  DIFPRM  for-  a  'specified  wavelength. 

DALFA:  Diffusion  parameter  defined  in  Appendix 

A-1.2.  (unitless) 

DBETA:  Diffusion  parameter  defined  in  Appendix 

A-1.2.  (km) 

DKAPPA:  Diffusion  parameter  defined  in  Appendix 

A-1.2.  (km-1) 

3-3.10  COMMON/ FLAGS / I TFL AG 

This  common  is  currently  used  only  to  avoid  calculation  of  the 
transfer  coefficient  in  function  TRANSF  in  the  statistical  calculations 
when  only  the  thermal  emission  is  needed. 

ITFLAG:  ITFLAG  turns  off  the  transfer  coefficient 

computation  as  described  above. 

3-3.11  COMMON /FLUX/FIN( 8,16), FOUT( 8 , 16 ) 

This  common  provides  communication  between  subroutine  CL0UD3 , 
which  generates  the  input  FIN  to  subroutine  DIFFUS  and  uses  the  output 
FOUT,  and  the  subroutine  DIFFUS.  DIFFUS  effects  the  diffusion  solution 
to  get  the  photon  flux  output  at  each  of  the  facets,  given  the  input 
flux. 
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FIN:  This  is  47t  times  the  photons  km  sec  that 

are  incident  normally  on  each  of  the  cloud 
facets  which  are  visible  to  the  unit  strength, 
isotropic  source. 

FOUT:  This  is  4tr  times  the  photon  flux  (radiant 

emittance)  that  is  transmitted  by  each  facet, 
according  to  the  diffusion  theory  approximation. 

( 4  tt  photons/ (km2  sec  )  ,  for  unit-strength  source). 

3-3.12  COMMON/ I NTRCPT / ALTLK ( 20 ) , THETALK( 20 ) , PH ILK ( 20 ) 

This  common  provides  communication  between  subroutine  CL0UD2 
and  the  calling  program  which  is  requesting  the  geographic  location 
of  the  intersection  of  the  line-of-sight  and  the  cloud. 

ALTLK:  The  altitude  at  which  the  line-of-sight 

intercepts  the  Kth  cloud,  (km) 

THETALK:  The  colatitude  of  the  Kth  cloud  intercept. 

( degrees ) 

PHILK:  The  east  longitude  of  the  Kth  cloud 

intercept,  (degrees) 

3-3.13  COMMON/ LOS /HLOS (1,64), ULOS (1,64), VLOS (1,64), WLOS (1,64) 

This  common  contains  lines-of-sight  to  all  facets  of  deter¬ 
ministic  clouds  which  provide  reflectance  and  to  those  facets  which 
have  the  highest  transmission.  These  are  provided  following  a  call  to 
subroutine  LOS  with  the  argument  K  for  the  Kth  cloud. 

HLOS:  Altitude  at  which  the  line-of-sight  from 

the  detector  intercepts  the  Jth  facet  in 
the  Kth  cloud,  (km,  above  mean  sea  level) 

ULOS , VLOS , WLOS :  Direction  cosines  in  the  earth-centered 
Cartesian  system,  of  the  above-described 
lines-of-sight . 

Each  of  the  vectors  ( ULOS , VLOS , WLOS )  provides  to  the  code 
user  a  potential  line-of-sight  which  is  separated  by  approximately  the 
square  root  of  the  f ield-of-view  from  another  possible  line-of-sight. 

3-3.14  COMMON / MATERL/MSM ,DDM( 7 ) , SRFALT 

MSM :  Integer  defining  susface  type. 

=1  for  Lambertian  diffuse  surface  with  spectrally- 
independent  reflectance  set  by  DD(1)  and 
emissivity  by  (l.-DD(l)). 

=2  for  water 

=3  for  snow 
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=4  for  sand 
=5  for  soil 
=6  for  foilage 
=7  for  urban  material 

DDM(M) :  Additional  descriptor  for  selected  surface  material. 

See  comments  in  subroutine  ESURF,  a  routine  pro¬ 
vided  by  the  model  documented  in  Reference  6. 

SRFALT:  Surface  altitude,  (km) 

3-3.15  COMMON /PAR AMS/IN, I OUT , P I , RE , DTR , RTD 

This  common  contains  program  constants,  set  in  program  ROSCOER, 
which  are  frequently  used  in  numerous  routines.  The  constants  can  be 
replaced  by  equivalencing  or  renaming  variables. 

IN:  The  integer  variable  name  used  to  identify 

the  input  file  unit.  It  is  not  used 
consistently. 

IOUT:  The  integer  variable  name  which  identifies  the 

output  file  unit.  It  is  not  used  consistently. 

PI :  The  value  of  PI  as  used  throughout  the  program. 

It  has  the  numerical  value  of  3.14159265. 


RE:  The  value  used  for  the  Earth's  radius  and  has 

a  numerical  value  of  6371.03  km. 

DTR:  This  is  the  conversion  factor  to  convert  from 

degrees  to  radians  and  has  a  numerical  value 
of  0.01745329251994. 

RTD:  This  is  the  conversion  factor  from  radians  to 

degrees  and  has  a  numerical  value  of  57.2957795131. 

3-3 . 16  COMMON /PTE/PI , P2 , P3 , P4 , T1 , T2 , T3 , T4 , SI , S2 , S3 , S4 , El , E2 , E3 , E4 


P1...P4:  Probability  fractiles  at  which  transfer 

coefficients  and  emission  radiances  are  desired. 
Set  in  program  ROSCOER. 

T1...T4:  Transfer  coefficients  at  12-km  altitude  for 

point  source.  This  is  the  ratio  of  a  radiance 
to  an  irrad'iance.  (sr-l) 


S1...S4:  Transfer  coefficients  at  12-km  altitude  for 

sun  as  a  source.  (sr~l) 


E1...E2:  Thermal  emission  (spectral  radiance)  from  top 

layer  of  clouds  or  the  ground  (for  the  cloud- 
free  line-of-sight) .  (W/(km^  sr  pm)) 
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PCFLOS:  Probability  of  the  detector  seeing  the 

artificial  source,  either  with  or  without 
clouds  being  present,  provided  the  source 
is  within  the  f ield-of-view. 

3-3 . 17  COMMON/SANDD/XS , YS , ZS , XD , YD , ZD , UL , VL , WL 

The  common  SANDD  contains  the  Cartesian  coordinates  of  the 
location  of  both  the  source  and  detector,  and  the  direction  cosines  of 
the  line-of-sight . 

XS,YS,ZS:  Cartesian  coordinates  of  the  source  location 

in  the  earth-centered  system.  The  x-direction 
points  to  zero  degrees  longitude  and  ninety 
degrees  colatitude.  The  z-direction  is  through 
the  north  pole,  (km) 

XD,YD,ZD:  Cartesian  coordinates  of  the  detector  location 

in  the  earth-centered  system,  (km) 

UL,VL,WL:  Direction  cosines  of  the  line-of-sight 

vector  in  the  earth-centered  system. 

3-3.18  COMMON/SOLARP/SOLLAT, SOLLON 

SOLLAT:  Solar  latitude  in  earth-centered  polar 

coordinates  with  positive  being  north  of 
the  equator,  (radians) 

SOLLON:  Solar  (east)  longitude  in  earth-centered 

polar  coordinates.  (radians) 

3-3.19  COMMON /SORCE/NSORCE , HSORCE( 1 ) , THETAS ( 1 ) ,PHIS(1) ,RSORCE(l) 

NSORCE:  This  variable  is  not  utilized  in  this  version 

of  the  code,  and  it  is  always  set  to  one. 

HSORCE:  Altitude  of  source,  (km) 

THETAS:  Colatitude  of  source,  (degrees) 

PHIS:  East  longitude  of  source,  (degrees) 

RSORCE:  Radius  of  source,  (km) 

3-3 . 20  COMMON /XXllll/KCLDK , CLDLAM 

KCLDK:  A  test-flag  used  in  subroutine  CLOUD3  to 

determine  whether  CLOUD3  is  now  being  called 
with  the  same  cloud-number  K  as  for  its  last 
call . 

CLDLAM:  A  test  quantity  used  in  subroutine  CL0UD3  to 

determine  whether  CLOUD3  is  now  being  called 
with  the  same  wavelength  ALAM  as  for  its  last 
call . 
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Note:  If  K  and  ALAM  inputted  to  subroutine  CL0UD3  have  the  same  values 

as  in  the  last  call,  the  subroutine  bypasses  its  call  to  subroutine 
DIFFUS(K, ALAM) .  Both  KCLDK  and  CLDLAM  are  initialized  in  the  main  program 
to  the  arbitrary,  out-of-range  values  of  99  and  -99.,  respectively. 
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SECTION  4 

DATA  REQUIREMENTS  AND  SAMPLE  PROBLEMS 

4-1  DATA  REQUIRED  FOR  STATISTICAL  AND  DETERMINISTIC  CLOUDS 

Required  user-supplied  data  for  transmission  of  IR  signatures 
through  clouds  are  read  by  CLOUDO .  Many  of  the  data,  as  outlined  in 
the  following  paragraphs,  are  stored  in  the  program  and  the  code  need 
only  be  directed  by  an  appropriate  flag  to  use  them.  Consequently, 
only  a  few  directives  need  be  inputted  to  run  standard  cases.  Data 
cataloged  in  the  program  are  described  first,  followed  by  descriptions 
and  illustrations  for  user-supplied  data. 

Data  which  are  cataloged  are  accessible  to  both  the  statis¬ 
tical  and  deterministic  submodules.  Only  a  limited  number  of  sets  of 
data  are  included  in  the  code  to  demonstrate  the  utility  of  the  code 
without  its  becoming  unwieldy. 

Cloud  descriptions  in  terms  of  mean  water  droplet  radius  (urn), 

_3 

droplet  density  (cm  ),  base  altitude  (km),  and  thickness  (km)  are 

incorporated  into  the  code  via  BLOCK  DATA.  These  data  are  tabulated 

in  Table  2-1.  Microscopic  extinction  and  scattering  cross-sections 
2 

(qm  )  and  average  cosine  of  the  scattering  angle,  averaged  (as 
described  in  Reference  2)  over  droplet  size  distributions  from  Table 
2-1,  are  also  included  in  BLOCK  DATA  for  14  cloud  types  and  10  wave¬ 
lengths  in  the  2-  to  5-qm  range.  No  provision  has  been  made  to  replace 
or  adjust  these  parameters  for  user-supplied  data.  A  modification 
could  be  implemented  to  provide  such  an  option. 

The  only  bidirectional  reflectance  data  that  are  in  the  code 
are  those  for  a  stratocumulus  cloud,  a  frequently  occurring  cloud  type 
according  to  the  NASA  cloud-data  bank.  The  tables  in  function  CLDBDR 
are  rather  lengthy.  Some  considerable  effort  has  been  expended  to  fit 
the  bidirectional  reflectance  data  for  inclusion  in  the  code;  however, 
the  code  contains  only  tabulated  data. 

The  cloud-frequency  data  are  stored  in  an  array  (CFREQ)  in 
BLOCK  DATA  dimensioned  (14  +  3)  by  4  by  (1C  +  1).  For  each  of  the  10 
(or  11  including  the  user-option  case)  region-averaged  data  sets,  there 


are  14  cloud-frequency  numbers  and  three  layer-frequency  numbers  for 
each  of  the  four  nonclear  coverage  categories.  Each  specification 
requires  five  coverage  categories,  of  which  four  are  independent,  rang¬ 
ing  from  completely  clear  to  complete  cloud-cover. 

Nominal  statistical-cloud  data  sets  are  included  in  the  pro¬ 
gram  which  include  averages  over  summer,  winter,  day,  night,  and  annual 
for  a  dry  and  a  moist  region.  They  are  arranged  as  in  Table  4-1. 

4-2.  USER-REQUIRED  AND  USER-OPTIONAL  DATA  FOR  NATURAL  CLOUD 

MODULES 

For  utilizing  the  ROSCOE-IR  natural  cloud  module  in  a  stand¬ 
alone  mode,  several  data  are  required  to  be  read  by  the  driver  program 
ROSCOER.  ROSCOE-IR  will  provide  these  data  to  the  natural  cloud  module 
as  implemented  in  the  ROSCOE-IR  code.  The  requirements  for  the  driver 
program  will  be  described  first,  followed  by  the  description  of  inputs 
read  by  subroutine  CLOUDO.  Sufficient  information  is  included  to 
illustrate  the  options  and  to  provide  test  cases  for  the  stand-alone  mode. 

Detector  and  source  locations  are  read  by  the  driver  program 
ROSCOER.  The  first  input  record  contains  H,  0,  $ ,  and  FOV  for  the 
detector;  the  second  record  contains  the  same  information  (less  the  FOV) 
for  the  source.  H  is  the  altitude  in  kilometers,  4>  is  the  colatitude 

(degrees),  and  is  the  east  longitude  (degrees).  FOV  is  the  field-of- 
2 

view  (km  )  for  the  detector.  The  third  input  record  contains  the  flag 
ISUN  which  is  zero  if  the  only  source  to  be  considered  is  the  artificial 
source  and  which  is  unity  if  the  sun  is  to  be  considered  as  a  source  in 
addition  to  the  artificial  source.  The  solar  location  is  specified  in 
a  data  statement  in  ROSCOER.  The  fourth  input  record  is  the  wavelength 
ALAM. 

Subroutine  CLOUDO,  for  which  a  flow  diagram  is  shown  in 
Figure  4-1,  reads  all  of  the  user-supplied  data  and  sets  the  proper 
flags  for  subsequent  code  execution.  For  the  statistical  case,  only 
data  are  read;  the  deterministic  mode  (MODE=0)  also  causes  the  geometry 
routines  (CLDGEOM,  and,  in  turn,  DSURF)  to  be  called. 

The  flags  are  set  by  Card  1  and  Card  2. 

Card  1:  MODE  (110  format) 

=  0  Deterministic  model 

=  1  Statistical  model 

(branch  here  to  required  submodule) 


Table  4-1.  Cloud  data  bank  contained  in  ROSCOE  natural 
cloud  program . 


Index  for 
Statistical 
Cloud  Set 

Region 

Type  of 

Average 

1 

~  a 
Dry 

Annual 

2 

r.  a 
Dry 

Winter 

night 

3 

Drya 

Winter 

day 

4 

Drya 

Summer 

night 

5 

Drya 

Summer 

day 

6 

Moist*5 

Annual 

7 

Moist*5 

Winter 

night 

8 

Moist*5 

Winter 

day 

9 

Moist*5 

Summer 

night 

10 

Moist*5 

Summer 

day 

11 

Array  reserved  for  user-supplied 
data 

ci  b 

’  Regions  4  and  11  in  Reference  1  nomenclature, 
i.e.,  continental  and  maritime. 


CLOUDO 


Figure  4-1. 


Flow  chart  of  subroutine  CLOUDO 
(Continued  on  next  page.). 


Read  (Card  2) 

KMODEL ,  KCOVER,  KLAYER 
KTYPL ,  KTYPM,  KTYPH 


KMODEL 


Read  (Card  2a) 
CCOVER  (1, 11), 1  =  1, 5 


Read  (Cards  2b) 

< -  CFREQ  (I,J,11) 

1  =  1 ,17;  J=1  ,4 

Options  to  override  program  data 
for  KMODEL  =  1,10 


(Total  of 
4  cards) 


KCOVER  >0  (Card  2a) 
CCOVER  (I, KMODEL)  1=1,5 


KLAYER  >0  (Cards  2c) 
CFREQ  ( I, J, KMODEL) 

1  =  15,17;  J=1 .4 _ 


KTYPL  >0  (Card  2d) 
CTYPL  (I),  1=1,8 


KTYPM  >0  (Card  2d) 
CTYPM  (I),  1  =  1 ,4 


KTYPH  >0  (Card  2d) 
CTYPH  (I),  1=1,5 


(Total  of 
4  cards) 


Figure  4-1.  Flow  chart  of  subroutine  CLOUDO  (Cont . ) 
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4-2.1 


Statistical  Mode  (M0DE=1)  Data 

Card  2:  KMODEL ,  KCOVER,  KLAYER,  KTYPL ,  KTYPM ,  KTYPH 

(615  format) 

KMODEL  =1  to  10,  cloud  model  (see  Table  4-1) 

=  11,  user-supplied  data 

A  nonzero  value  of  one  or  more  of  the  following  five  flags 
allows  the  user  to  override  selected  portions  of  the  program  data  for  a 
given  KMODEL=l , 10 . 

KCOVER  =  0  no  action  required 

=  1  read  cloud-coverage-frequency  data  (Card  2a) 

( ( CCOVER (I, K), 1=1, 5) ,  K=KM0DEL . LE . 1 0 ) 

KLAYER  =  0  no  action  required 

=  1  read  layer  frequency  (Cards  2c) 

( ( (CFREQ( I , J , K) ,  I  =  15,17),  J  =  1,4),  K  =  KMODEL . LE . 10) 

See  description  for  Cards  2d  for  a  nonzero  value  of  any  of 
the  following  three  flags. 

KTYPL  =  0  no  action  required 

=  1  read  low-layer  cloud  type  desired  (Card  2d) 

KTYPM  =  0  no  action  required 

=  1  read  middle-layer  cloud  type  desired 

KTYPH  =  0  no  action  required 

=  1  read  high-layer  cloud  type  desired 

Card  2a:  for  KMODEL  =  K,  K  =  1,11 
CCOVER(l.K),  ...  ,  CCOVER(5,K)  (5F4.3  format) 

CC0VER( 1 , K)  is  the  frequency  of  cloud-coverage  category  1  (clear) 
CC0VER(2,K)  is  the  frequency  of  cloud-coverage  category  2  (0.3  cloud 
cover) 

CC0VER(3,K)  is  the  frequency  of  cloud-coverage  category  3  (0.5  cloud 
cover) 

CC0VER(4,K)  is  the  frequency  of  cloud-coverage  category  4  (0.8  cloud 
cover) 

CC0VER(5,k)  is  the  frequency  of  cloud-coverage  category  5  (complete 
cover)  ^ 

where  ^  CC0VER(I,K)  =1.0  . 


Cards  2b:  for  K.MODEL  =11  (four  cards) 


CFREQ( I , J, 11 ) 

1  J 

c  1  <  i7,  : 

L  < 

J  <  4 

(17F4.3  format) 

CFREQ(1, J,K) 

is 

frequency 

of 

CU1  ^ 

CFREQ(2, J,K) 

is 

frequency 

of 

CU2 

CFREQ(3, J,K) 

is 

frequency 

of 

CU3 

CFREQ(4 , J,K) 

is 

frequency 

of 

sc  ' 

?  low-altitude 

CFREQ(5, J.K) 

is 

frequency 

of 

ST  1 

clouds 

CFREQ( 6 , J , K) 

is 

frequency 

of 

CBl 

CFREQ(7, J,K) 

is 

frequency 

of 

NS 

CFREQ( 8 , J , K) 

is 

frequency 

of 

AC  j 

| 

CFREQ(9 , J , K) 

is 

frequency 

of 

AS1 

mid-altitude 

clouds 

CFREQ(10, J,K) 

is 

frequency 

of 

CB2  ' 

1 

CFREQ( 11 , J , K) 

is 

frequency 

of 

CB3  'j 

CFREQ( 12 , J , K) 

is 

frequency 

of 

AS2  | 

^  high-altitude 

CFREQ( 13 , J , K) 

is 

frequency 

of 

CI  1 

:  clouds 

CFREQ( 14 , J , K) 

is 

frequency 

of 

CS  J 

where 


14 


CFREQ(  I ,  J  ,  K) 

1  =  1 


=  1, 


CFREQ( 15 , J , K) 

is 

frequency 

of 

1 

layer 

CFREQ( 16 , J , K) 

is 

frequency 

of 

2 

layers 

CFREQ( 17 , J , K ) 

is 

f requebcy 

of 

3 

layers 

where 


17 


y;  CFREQ(  I  ,  J  ,K)  =  1. 


1  =  15 


Cards  2c:  for  K  =  KMODEL  =  1,10  (3F4.3  format)  (four  cards) 

Cards  2c  provide  a  mechanism  to  change  layer  frequencies  for  KMODEL  <11, 
i.e.,  (((CFREQ(I ,J, KMODEL) ,  1  =  15 , 17 )  , J=1 , 4 ) ,  KMODEL< 11)  without  chang¬ 
ing  the  individual  cloud-type  frequencies. 
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_ ■  ______ 


CFREQ( I , J , K) ,  1=15,16,17  and  J=l,2,3,4 


1=15  is  the  frequency  of  1  layer 

=16  is  the  frequency  of  2  layers 

=17  is  the  frequency  of  3  layers 

Cards  2d:  cloud-type  specification  (8A5  format) 

(one,  two,  or  three  cards) 

Read  under  8A5  format  (a  blank  field  rejects  that  cloud 
type).  Data  must  be  left-justified. 


For  KTYPL  =  1  (one  card): 

CTYPL(I),  1=1,8 
CTYPL(l)  =  CU1 
CTYPL( 2 )  =  CU2 
CTYPL( 3)  =  CU3 
CTYPL( 4 )  =  SC 
CTYPL( 5)  =  ST 
CTYPL(6 )  =  CB1 
CTYPL( 7 )  =  NS 
CTYPL( 8)  =  CLR  =  clear 


For  KTYPM  =  1  (one  card): 

CTYPM(I),  I  =  1,4 
CTYPM(l)  =  AC 
CTYPM ( 2 )  =  AS1 
CTYPM( 3)  =  CB2 
CTYPM ( 4 )  =  CLR 


For  KTYPH  =  1  (one  card): 

CTYPH(I),  I  =  1,5 
CTYPH(l)  =  CB3 
CTYPH( 2 )  =  AS2 
CTYPH(  3)  =  Cl 
CTYPH( 4 )  =  CS 
CTYPH( 5 )  =  CLR 


By  utilizing  combinations  of  Card  1  through  Cards  2d,  one  can 
accomplish  any  degree  of  specification  desired.  However,  for  a  complete 
specification,  the  problem  reduces  to  that  of  a  deterministic  calculation 
where  the  thicknesses  and  altitudes  come  from  the  program  data  bank  and 
the  cloud  is  infinite  in  lateral  extent. 


4-2.2 


Deterministic  Mode  (MODE=0)  Data 


Card  2:  NCLOUD  (110  format) 

Card  3:  KCLOUD(I),  HCLOUD(I),  THETCD(I),  PHICD(I) 
CLOUDA(I),  CLOUDB(I),  CLOUDC(I),  CLDPSI(I) 
where  1  <_  I  <_  NCLOUD 
(12,  E8.3,  6E10.4  format) 


KCLOUD(I) 
HCLOUD(I) 
THETCD( I ) 
PH I CD ( I ) 
CLOUDA(I) 
CL0UDB( I ) 
CLOUDC(I) 
CLDPSI ( I ) 


Cloud-type  index  of  Ith  cloud  (1,2,..., 14) 

Altitude  of  Ith  cloud  (km) 

Colatitude  of  Ith  cloud  (degrees) 

East  longitude  of  Ith  cloud  (degress) 

Semi-axis  in  x-direction  of  Ith  cloud  (km) 

Semi-axis  in  y-direction  of  Ith  cloud  (km) 

Semi-axis  in  z-direction  of  Ith  cloud  (km) 

Orientation  of  'the  a-axis  with  respect  to 
the  east-west  direction,  i.e.,  ip  orients 
the  longest  (a-axis)  at  an  angle  of 
degrees  north  of  east . 


4-3.  SAMPLE  PROBLEMS 

A  number  of  sample  computer  runs  have  been  made  using  the 
ROSCOE  natural  cloud  program.  Several  of  these  are  included  here  to 
facilitate  subsequent  code  comparisons.  Listings  of  input  cards 
necessary  for  each  problem  are  provided. 

Some  modifications  of  the  code  were  required  to  obtain  the 
data  output  shown  in  the  samples.  Code  calculations  were  affected  only 
as  outlined  below.  Print  statements  were  placed  in  the  code  to 
intercept  calculations  in  a  convenient  manner  for  presentation.  These 
are  given  in  Appendix  C  to  facilitate  reproducing  the  output  for  each 
problem. 


4-3.1  Problems  1  to  4 :  Use  of  the  Statistical  Cloud  Submodel  in 
a  Limited,  Deterministic  Mode  (with  KMODE1=1,10) 

The  purpose  of  these  four  problems,  which  form  a  set,  is 
(1)  to  illustrate  use  of  the  statistical  cloud  submodel  in  a  limited, 
deterministic  mode  (with  KMODEL=1,10)  and  (2)  to  show  that  diffusion 
calculations  performed  by  the  code  for  a  point  source  beneath  either 
(a)  single-layer  (planar)  clouds  or  (b)  multilayer  (planar)  clouds 
are  in  good  agreement  with  exact,  analytic  diffusion  calculations. 


To  help  provide  the  motivation  for  a  portion  of  these  pro¬ 
blems,  we  note  that  the  total  transmittance  T  for  a  planar  slab  is  given 
by  Equation  (31a)  of  Appendix  A,  i.e., 

T  =  Q  [cosh  KL  +  i  +  & 

where  Q,  the  total  flux  into  the  slab,  can  be  set  to  unity.  To  get  a 
computer  code  result  that  can  be  compared  with  this  transmittance,  we 
must  take  several  steps: 

a.  Instruct  the  code  to  treat  a  single  planar  cloud  above 
a  point  source. 

b.  Compute  the  radial  distribution  of  the  transmission 
transfer  coefficients,  without  including  air 
attenuation . 

c:.  Convert  this  distribution  of  transfer  coefficients 
to  a  radial  distribution  of  radiant  exitance. 

d.  Integrate  this  radial  distribution  of  radiant  exitance 
across  the  cloud  top  and  divide  by  the  spatially- 
integrated  flux  into  the  slab  to  obtain  a  total 
transmittance . 

The  quantity  that  the  code  was  actually  instructed  to  compute 

is  one  we  shall  refer  to  as  a  modified  transfer  coefficient,  TC' ,  in 
-2  -1 

units  of  km  sr  .  These  units  are  those  of  the  transfer  coefficient 

(sr-1)  divided  by  the  square  of  the  distance  from  the  source  to  the 

12-km  altitude  transfer  point.  These  units  are  equivalent  to  describing 

-2  -1 

TC'  as  being  4tt  times  the  radiance  (W  km  sr  )  directed  toward  the 
downward-looking  radial-scanning  detector,  for  a  unit  strength  isotopic 
source.  (Alternatively,  TC’  is  4  times  the  radiant  exitance.)  This 
modified  transfer  coefficient,  4vL ,  was  divided  by  4tt  and  multiplied 
by  v  to  obtain  the  radiant  exitance  and  was  subsequently  integrated 
across  the  cloud-top  radial  distribution.  Division  by  the  source 
power  incident  on  the  underside  of  the  cloud  gave  the  total  trans¬ 
mittance  T: 

T  =  JJ. - i - T  r  *)  2irrdr  . 

(i  source  power)  JQ  V  4tt  / 


^  sinh  KL] 
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The  information  on  the  input  cards  read  by  ROSCOER  for  the 
detector,  source,  and  wavelength  and  by  subroutine  CLOUDO  for  Cards  1,2 
and  2d  (cf.  Figure  4-1  and  corresponding  text)  is  given  in  Figure  4-2 
for  the  four  sample  problems.  Each  of  the  first  three  problems  treats 
a  different  single  cloud  type,  namely,  Type  4  (Sc),  Type  8  (Ac),  and 
Type  13  (Ci),  respectively:  the  fourth  problem  treats  the  three-layer 
configuration  of  these  same  cloud  types.  We  shall  discuss  the  input 
cards  later. 

For  each  of  the  four  problems,  and  for  a  range  of  lateral 
displacements  in  0 . 25-km  steps,  the  values  of  the  modified  transfer 
coefficient,  4irL,  are  given  in  Columns  2  to  5  of  Table  4-2.  These 
quantities  are  plotted  in  Figure  4-3.  Columns  6  to  9  of  Table  4-2 
give  the  corresponding  values  of  the  radial  integrals  for  the  trans¬ 
mittance  . 

In  Figure  4-3  the  curve  for  cloud  Sc  is  not  shown  at  larger 
distances  because  it  soon  develops  wiggles  and,  starting  at  35.75  km, 
an  occasional  negative  value.  These  features  result  from  limitations 
of  the  internal  integration  scheme  which  is  intended  to  be  efficient 
in  the  lateral  range  of  importance  but  not  at  very  large  distances. 

The  curve  for  cloud  Ac  begins  to  wiggle  a  bit  at  about  20  km,  that  for 
cloud  Ci  at  about  30  km,  and  that  for  the  combined  set  of  clouds 
(Sc  +  Ac  +  Ci)  at  7  km. 

The  analytic  results  for  the  transmittance  and  reflectance  of 
each  of  the  three  layers  are: 

T(Sc )  =  0.820  x  10-2  R(Sc)  =  0.653 

T( Ac )  =  0.113  x  10_1  R( Ac )  =  0.637 

T(Ci )  =  0.127  x  10"6  R(Ci )  =  0.0761 

The  analytic  result  for  the  transmittance  of  the  three-layer 
configuration  is 

T( Sc  +  Ac  +  Ci)  =  0.211  x  10-10 
which  was  evaluated  from  the  formula,  derived  in  Appendix  D, 
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Figure  4-2.  Input  data  for  Problems  1  to  4 . 
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where 


TV  =  transmittance  of  ith  layer 

R^  =  reflectance  of  ith  layer. 

The  total  reflectance  R  for  a  planar  slab  is  given  by  Equation  (31b)  of 
Appendix  A  as 

R  =  T  •  *  (ke  -  ?)  sinh  KL  • 

We  can  now  compare  our  code  results  with  the  analytical  results. 
This  comparison  is  given  in  Table  4-3  where  we  give,  as  a  function  of 
the  lateral  displacement  from  the  source,  the  ratio  of  the  code-computed 
integrated-transmittance  to  the  analytic  integrated-transmittance.  We 
see  that  for  cloud  Sc  the  code  result  is  less  than  the  analytic  result 
by  only  7%,  3%,  and  2%  at  distances  of  5,  10,  and  15  km.  For  cloud  Ac 
and  cloud  Ci  the  code  results  are  small  by  5%  and  15%  even  at  40  km. 

For  the  three-layer  configuration,  the  code  result  is  small  by  only  15% 
at  5  km  and  2%  at  10  km. 

We  conclude  that,  whereas  the  comparisons  between  the  code 
results  and  the  analytic  results  for  the  integrated  transmittance  are 
not  perfect  for  these  examples,  the  comparison  is  of  sufficiently  high 
quality  for  the  radial  distances  of  importance  that  one  can  have  a  high 
confidence  in  the  diffusion  calculations  performed  in  the  statistical 
cloud  submodel. 

We  return  now  to  the  input  cards  for  Problems  1  to  4.  We 
ha^e  two  general  choices  for  KMODEL  which  can  be  either  (a)  any  value 
in  the  range  from  1  to  10  or  (b)  a  value  of  11. 

If  KMODEL  is  less  than  11,  it  doesn't  matter,  for  Problems  1 
to  4,  which  of  the  allowed  values  we  choose.  We  can  understand  that 
fact  by  reviewing  how  the  code  operates  for,  say,  Problem  1: 

(1)  ROSCOER  reads  the  detector,  source,  and  wavelength 
cards . 

(2)  Subroutine  CLOUDO  reads  the  Type-1  and  Type-2  cards 
and  three  Type-2d  cards. 

(2.1)  On  the  Type-2  card  we  set  KCOVER  and  KLAYER  to  0 

since  we  have  no  need  to  override  the  stored  data  for 
either  the  CCOVER  or  CFREQ  arrays;  in  fact,  we  really 
don't  use  the  stored  data  for  these  arrays. 
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Table  4-3.  Ratio  of  T(code)-  to  T(analytic)-results  for  sample 
Problems  1  to  4 . 


Cloud  Configuration 
Distance  - 


(km)  Sc  Ac  Ci  Sc  +  Ac  +  Ci 


5 

0.934 

0.543 

0.228 

0.850 

10 

0.968 

0.792 

0.485 

0.979 

15 

0.978 

0.866 

0.630 

20 

0.901 

0.714 

25 

0.922 

0.769 

30 

0.936 

0.806 

35 

0.946 

0.834 

40 

0.953 

0.855 

(2.2) 

On  the  Type 

-2  card  we  set  KTYPL 

,  KTYPM ,  and  KTYPH  to 

1  to  allow 

the  code  to 

read  in 

three  Type-2d  cards. 

(3) 

The  three  Type-2d  cards 

specify 

that  the  only  cloud 

configuration  of  interest  is  that  for  which  the  low- 
altitude  layer  has  cloud  Type-4  (Sc)  and  the  mid- 
and  high-altitude  layers  are  clear. 

(4)  ROSCOER  calls  subroutine  SCLOUD. 

(5)  Subroutine  SCLOUD  calls  subroutine  CLDWT. 

(6)  Subroutine  CLDWT' s  calculations  are  very  much 
restricted  compared  with  a  normal  run. 

(6.1)  Consider  Pass-1  through  the  DO-500  loop  (11=2).  In 
the  DO-101  loop,  the  tests  are  satisfied  only  for 
LL=4.  PIJKL  is  computed  and  IDX  is  set  to  1.  TRANS(l) 
and  WT(160)  are  also  computed  (but  never  used).  None 
of  the  loops  after  DO-101  becomes  effective. 

(6.2)  Consider  Pass-2  through  the  D0-500  loop  (11=3).  Again, 
in  DO-101  loop,  tests  are  satisfied  only  for  LL=4. 

PIJKL  is  computed  and  IDX  is  again  set  to  1.  This 
time,  computations  of  TRANS ( 1 )  and  EMISS(l)  are 
skipped  because  they  were  computed  on  Pass-1.  The 
Pass-1  values  of  WT(1)  and  WT(160)  are  augmented  with 
the  Pass-2  values  (but,  again,  will  never  be  used). 

(6.3)  Pass-3  and  Pass-4  through  D0-500  loop  are  just 
like  Pass-2. 

(6.4)  Subroutine  CLDWT  returns  to  subroutine  SCLOUD  values 
for  TRANS ( 1 ) ,  EMISS(l),  WT(1),  and  WT(2)  =  WT(160). 

(7.1)  In  subroutine* SCLOUD  we  exploit  common  PTE  by 

arbitrarily  setting  T4  equal  to  TRANS ( 1 ) ,  the  modified 
transfer  coefficient. 
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(7.2)  The  code  effectively  returns  to  ROSCOER  without 

doing  any  of  the  usual  calculations  or  sorting  in 
subroutine  SCLOUD. 

(8)  ROSCOER  then  performs  the  radial  integration  over 
the  cloud-top  surface. 

We  now  consider  what  would  happen  if  we  inputted  11  for 
KMODEL.  First,  we  note  that  the  code  would  expect  to  read  a  Type-2a 
card  and  four  Type-2b  cards.  The  natural  selection  of  values  for 
these  cards  would  be  as  follows: 

For  Problems  1  to  4 : 


CCOVER  (1,11) 


J 


0.0,  1=1,4 


1.0,  1=5 


CFREQ( I , J , 11 )  =  0.0  (1=1,17),  J=1 , 3 


and 


Problem 


CFREQ( 1,4,11)  = 


1.0  1=4  &  15 
0.0  Otherwise 

1.0  1=8  &  15 
0.0  Otherwise 

1.0  1=13  &  15 
0.0  Otherwise 

1.0  1=4,8,13  &  17 
0.0  Otherwise 


However,  it  really  wouldn't  matter  what  is  read  in  for  CCOVER  and 
CFREQ  since  we  do  not  use  them  in  Problems  1  to  4. 

Subroutine  CLDWT  will  return  the  complete  arrays  for  TRANS 
and  EMISS  but  the  WT  array  will  have  zeros  except  that: 


WT(4)fO 

and 

WT( 160)^0 

for 

Problem 

1 

WT( 8 )/0 

and 

WT( 160)^0 

for 

Problem 

2 

WT( 13)^0 

and 

WT( 160)^0 

for 

Problem 

3 

WT(114)?0 

1  and 

WT( 160)^0 

for 

Problem 

4 

Of  course,  we  are  interested  in  TRANS(IDX)  for  selected  values  of  IDX: 
i.e.,  IDX  =  4,8,13,  and  114  for  Problems  1,2,3,  and  4,  respectively. 

These  members  of  the  TRANS  array  can  and  should  be  intercepted  in 
subroutine  SCLOUD  immediately  after  the  return  from  subroutine  SCLOUD 
immediately  after  the  return  from  subroutine  CLDWT.  Failure  to  intercept 
at  this  point  would  allow  the  TRANS  array  to  be  sorted,  with  the  loss  of 
identity  of  the  particular  member  of  TRANS  that  is  of  interest. 

4-3.2  Problem  5:  Use  of  the  statistical  cloud  submodel  in  the 

optional  mode  with  KM0DEL=11 

The  purpose  of  this  problem  is  the  illustrate  (1)  use  of  the 
statistical  cloud  submodel  in  the  optional  mode  with  KM0DEL=11  and  (2) 
a  detailed  printout  of  the  complete  results  for  a  statistical  cloud 
problem. 

Use  of  the  cloud  model  with  KM0DEL=11  requires  the  user  to 
read  in  the  Type-2a  card  for  (CCOVER( I , 11 ) ,  1=1,5)  and  four  Type-2b  cards 
for  (CFREQ( I , J , 11 ) ,  1=1,17;  J=l,4),  as  seen  from  Figure  4-1.  The  input 
data  deck  is  shown  in  Figure  4-4. 

Selection  of  appropriate  data  for  use  with  KM0DEL=11,  i.e., 
for  the  CCOVER  and  CFREQ  arrays,  is  always  the  user's  responsibility. 

For  the  authors'  convenience,  and  for  illustration,  we  have  arbitrarily 
selected  the  data  from  that  given  for  KM0DEL=1,  as  can  be  seen  by  com¬ 
paring  with  BLOCK  DATA.  (In  this  instance,  because  we  have  used  data 
from  KM0DEL=1,  we  could  get  the  same  results  more  directly  by  not  inputting 
11  for  KMODEL  and  the  associated  Type-2a  and  Type-2b  cards  but  instead 
inputting  1  for  KMODEL. ) 

Because  we  have  inputted  a  1  for  the  source-type  flag  ISUN, 
the  code  will  first  do  the  problem  with  the  artificial  source  as 
specified  and  then  do  the  problem  with  the  sun  as  the  source.  The 
emission  radiance  will  be  calculated  in  each  case. 

The  output  results  for  the  artificial  source  are  shown  in 
Tables  4-4a  through  4-4d  and  those  for  the  solar  source  are  shown  in 
Table  4-5a  thorugh  4-5c. 

The  quantities  in  the  various  columns  in  Tables  4-4a  and  4-5a 

are  described  in  Table  4-6.  The  results  for  IDX=1,159  are  provided  by 

subroutine  CLDWT  for  the  159  different  cloud  configurations.  The 

-51 

entries  for  TRANS iTCOEF  with  values  of  'vlO  occur  when  the  argument 
of  the  sinh  and  cosh  functions  in  the  integral  being  evaluated 
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PROBLEM  5 


HD 

THETAD 

PNID 

FUV 

UETECTUR 

35793. 

90.00001 

•90.0001 

1.0 

HSORCE(l) 

ThCTAS(I) 

PMXSd) 

WKT 

SOURCE 

0.01 

90.000001 

-90.00 

1000. 

ISUN 

1 


ALAM  HAVELENCTh 

2.50 

MUOE  STATISTICAL  CARD  t 

1 

KMUDEL  KLATCR  KTTPM  USER  CLOUD  MODEL  CARO  2 

KCOVER  KTTPL  rfTVPH 
It  0  0  0  0  0 

1«1  Is2  I ■  J  1*4  1*5  CCOVERd*  11)  CARD  24 

.235.150. 070. 201. 343 

CFREQd.l.ll)  CARD  2B 

1«1  I »2  1*3  1*4  1*5  1*6  1*7  1*0  1*9  10  II  12  13  14  IS  16  17 

.046. 046. 046.1 02.006.005.001.217.005.005. 005.005.222.293*760.212.019 

CFREQd.2,11)  CARD  2B 

1*1  1*2  1*3  1*4  1*5  1*6  1*7  1*6  1*9  10  11  12  13  14  15  |6  17 

. 055. 055. 055. 121. 006. 006. 001. 243. 00*. 006. 006. 009d64. 260. 542. 377. 061 

CFREO(I.3,d)  CARD  2B 

1*1  1*2  1*3  1*4  1*5  1*6  1*7  1*6  1*9  10  11  12  13  14  15  16  17 

.036, 030.036, 169.ol5.01 0.000. 276.0 12.01 0.01 0*0 12. 104.266.41 5.435.150 

CFREQ ( 1 , 4, 11)  CARO  2B 

1*1  1*2  1*3  1*4  1*5  1*6  1*7  1*8  1*9  10  11  12  13  14  IS  16  17 

.010.010.016.336.155.010.030,197.023.010.010.023.017.156.448.423.126 


Figure  4-4.  Input  data  for  Problem  5. 
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Table  4-4a.  Initial  arrays  for  artificial  source  transfer-coefficients,  emission  radiances 
and  weights  in  sample  Problem  5.  Wavelength  =  2.5  pm. 
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Table  4-4b.  Sorted  array  for  artificial  source  transfer-coefficients 
(and  carried-along  weights)  for  sample  Problem  5. 


TAAISsMT  - 


0. 

0. 

•3S37E-51 

.56128-05 

•3637E-51 

•8001E-04 

•3637E-51 

.3308E-0S 

.36372-51 

•1092E-03 

.3S37E-51 

•3308E-05 

•3637E-51 

•6080E-04 

• 3637E-51 

.33082-05 

•3637E-51 

• S612E-05 

• 3637E-51 

• 24Q0E-04 

• 3637E-51 

.45182-05 

• 3637E-51 

•4781E-03 

• 36372-51 

.54462-04 

.3637E-51 

• 1 2288-03 

•3637E-51 

•10392-04 

• 3637E-51 

.19622-05 

•3637E-51 

• 8903E-04 

• 36372-51 

.58268-04 

• 3637E-51 

.12552-03 

• 3637E-51 

•1210E-02 

.43652-51 

• 5612E-05 

• 4365E-51 

•11328-03 

•4365E-51 

.18512-02 

•43652-51 

.12828-03 

•4365E-51 

• 1053E-04 

•4365E-51 

• 1 2822-03 

.4  36  52-51 

•  7  298E-04 

• 4365E-51 

•1132E-03 

• 4365E-51 

•7298E-04 

•4365E-51 

•11442-02 

•4365E-51 

•  7  298E-04 

.43658-51 

.56128-05 

•4365E-51 

• 10532-04 

.43658-51 

•56122-05 

•4365E-51 

•1798E-03 

.43652-51 

.17988-03 

.4365E-51 

.1132E-03 

•  4  365E-51 

•1 851E-02 

•  4  365E-51 

•  185  IE-02 

.43652-51 

• 3865E-03 

. 4365E-51 

. 86928-03 

• 4365E-51 

• 1 798E-03 

.43652-51 

.12288-03 

.4365E-61 

•12442-03 

•4365E-51 

•2573E-02 

•4365E-51 

.31402-03 

•4365E-51 

•7760E-02 

.43652-51 

• 2561E-02 

•4365E-51 

•81 14E-03 

• 4365E-51 

.16422-03 

•4365E-51 

.21562-02 

.4365E-51 

.39252-02 

•4365E-51 

•  5446E-0  4 

.4365E-51 

.26232-03 

•4  365E-51 

.10798-02 

.43658-51 

.10538-04 

.43652-51 

•1223E-03 

.4365E-51 

•  2 164E-02 

•4365E-51 

•  7  143E-Q4 

.43652-51 

.4009E-05 

•  43652-51 

• 1409E-02 

•4365E-51 

.26322-01 

.43652-51 

•23892-04 

.43658-51 

.20478-03 

•4365E-51 

.1282E-03 

.43658-51 

•861 5E-05 

.43652-51 

•  8  335E-04 

• 4365E-51 

•  o  +*0£-04 

• 4365E-51 

.3609E-04 

• 4365E-51 

.10392-04 

.43652-51 

•  1 174E-01 

• 4365E-51 

.20852-02 

•4365E-51 

•3133E-01 

•  41652-51 

.21522-01 

.43652-51 

.16478-02 

.43652-51 

.38212-02 

• 4365E-51 

•2547E-02 

•4365E-51 

. 1 255E-03 

•4365E-51 

.27852-03 

•4365E-51 

.10062-02 

•  4851 £-51 

•6359E-03 

• 4851E-51 

.23142-04 

.48512-51 

•30252-04 

.48518-51 

• 4050E-02 

• 48512-51 

•2314E-04 

•  4351 £-5 1 

.30252-04 

•4851E-51 

•9978E-04 

•48512-51 

•  23142-04 

. 4851E-51 

•4683E-04 

• 4851E-51 

.39252-03 

•48512-51 

•9815E-03 

.48512-51 

.92788-03 

• 4851E-51 

.297BE-03 

.45518-51 

•2602E-03 

.48512-51 

.15338-03 

.48512-51 

•7555E-05 

•48512-51 

•  1699E-0  3 

• 48518-51 

.77392-05 

.48512-51 

• 1 772E-04 

•72B0E-51 

.1092E-03 

•7280E-51 

•1786E-02 

•7280E-51 

.17982-03 

•  7  280E-51 

• 1 366E-03 

•  72dOE-5i 

.72688-03 

•35972-38 

.10922-03 

• 4506E-37 

.60802-04 

.45062-37 

•10922-03 

.61358-37 

.19922-02 

. 1 2002-23 

•2561E-02 

• 13262-23 

.40532-01 

•7354E-23 

. 1 144E-02 

.73542-23 

.25618-02 

•  1 339E-12 

•1184E-04 

•1339E-12 

• 3025E-04 

.37052-11 

•40082-05 

.37052-11 

.56128-05 

•8333E-10 

.27268-02 

•2307E-09 

.63592-03 

.32342-09 

•237yE-03 

.73938-39 

.19052-03 

.7393E-09 

• 6359E-03 

•8132E-09 

.68732-02 

•8534E-09 

.14502-02 

.14918-08 

.56978-02 

•25902-08 

.8001E-Q4 

•3062E-08 

.13662-03 

•58352-08 

.72688-03 

.76938-08 

.60012-04 

•7693E-08 

•  42332-0  4 

• 8740E-08 

.13282-02 

•8990E-08 

•4231 E-02 

.10838-07 

.17668-02 

.1384E-07 

•9278E-03 

.2199E-07 

.11672-01 

.31192-07 

•  3 140E-03 

.61928-07 

.92788-03 

•6192E-07 

• 2602E-Q3 

.8970E-07 

.16472-02 

.11602-06 

.19928-02 

• 1482E-06 

•10922-03 

-1553E-06 

•3821E-02 

.4S14E-06 

•6080E-04 

.461 4E-06 

.10922-03 

•8300E-06 

.17988-03 

• 4550E-05 

•3925E-02 

•5278E-05 

.17982-03 

.52782-05 

.12288-03 

.64748-04 

.26952-02 

.99462-04 

.1450E-01 

.2231E-03 

•37362-01 

.19012-02 

•2726E-02 

•4801E-02 

•426SE-02 

.6281E-02 

•2726E-02 

.62818-02 

.13282-02 

•68122-02 

.45192-01 

.46172-01 

.57412-01 

. 115dE+00 
.3798E+01 

•2271E-01 
• 4231E-02 

.12768+00 

.42312-02 

•1616E+01 

•53992+00 

.37988+01 

• 1992E-02 
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Table  4-4c.  Sorted  array  for  emission  radiances  (and  carried-along 
weights)  for  sample  Problem  5. 


EH1SS/VT1  = 


0. 

0. 

.11266*01 

•5612E-05 

•1126E*01 

•8001E-04 

•  U26E*01 

. 3308E-05 

.11266*01 

.80016-04 

.11266*01 

.33086-05 

•1126E+01 

•12106-02 

•  1 1 26E*01 

.33086-05 

•11 266*01 

•56126-05 

.11266*01 

.10926-03 

•1 126E*01 

•1 992E-02 

• 11266*01 

•7268E-03 

.11266*01 

• 1786E-02 

•1126E*01 

•1092E-03 

•  1 126E*01 

• 1366E-03 

• 1126E*01 

. 1092E-03 

.1126E*01 

.60806-04 

•  11 26fi^01 

.122RE-03 

•11266*01 

.13286-02 

.11266*01 

.56126-05 

.1126E*01 

•60806-04 

•11 26E*01 

•2400E-04 

si 126E*01 

.54466-04 

•1126E*01 

.47816-03 

•  1 126E*01 

• 4008E-05 

.11266*01 

.80016-04 

.11266*01 

.10396-04 

•11266*01 

.89036-04 

•1126E*01 

•1962E-05 

.11266*01 

•58266-04 

•11266*01 

.42336-04 

•11266*01 

.1255E-03 

•11 266*01 

•4519E-05 

.52966*01 

.72986-04 

•5296E*01 

•1 132E-03 

.52  96E*01 

.12826-03 

•5296E*01 

•5612E-05 

•5296E*0i 

•1053S-04 

•5296E*01 

•5612E-05 

•5296E*01 

•1132E-03 

•52966*01 

•7298E-04 

.52966*01 

.11322-03 

.52966*01 

.18516-02 

.52966*01 

•  561 2E-05 

•  5296E*  01 

• 1053E-04 

.52966*01 

• 1 282E-03 

•5296E*01 

•1798E-03 

•5296E*0i 

.25476-02 

.52966*01 

•25616-02 

•5296E*01 

.21526-01 

•5296E*01 

•7  298E-04 

•5296E*01 

•25616-02 

.5296E*01 

.25616-02 

•52966*01 

.38216-02 

•5296E*01 

.17986-03 

•5296E*01 

•1798E-03 

•5296E*01 

•1228E-03 

•5296E*01 

.1647E-02 

•5296E*01 

•1 144E-02 

.52966*01 

.11746-01 

• 5296E*0l 

•3133E-01 

•5296E*01 

.39252-02 

•5296E*01 

.21646-02 

.52966*01 

•4053E-01 

•52966*01 

•2573E-02 

•5296E*01 

.3140E-03 

.52966*01 

.26326-01 

•5296E*0i 

• 1228E-03 

. 5296E*01 

• 2623E-03 

.52966*01 

.11446-02 

•5296E*01 

.21566-02 

.52966*01 

•1282E-03 

•5296£*0l 

•10066-02 

•5296E*01 

•1851E-02 

.52966*01 

.10796-02 

• 5296E*01 

.10536-04 

• 5296E»01 

• 7760E-02 

.52966*01 

.54466-04 

.52966*01 

.12446-03 

•5296E*01 

.3865E-03 

.5296E*01 

•1B51E-02 

.52966*01 

.83352-04 

.52966*01 

•69406-04 

.52966*01 

•8114E-03 

•  52966*01 

• 4003E-05 

• 5296E*01 

.86156-05 

•5296E*01 

.36096-04 

•5296E*01 

•2389E-04 

•5296E*01 

•2047E-03 

•52966*01 

•1642E-03 

•52966*01 

.14096-02 

.52966*01 

.27856-03 

•52966*01 

• 1039E-04 

• 5296E*01 

,20856-02 

.52966*01 

.71436-04 

•5296E*01 

.12556-03 

.52966*01 

•8692E-03 

.10736*02 

.92786-03 

• 1073E*02 

• 2314S-04 

.10736*02 

.23I4E-04 

• I073E*02 

•6359E-03 

.10736*02 

• 3025E-04 

.10736*02 

.26026-03 

•1073E*02 

• 4050E-02 

• 10736*02 

•5697E-02 

.1073E*02 

•  231 4E-04 

.10736*02 

.635*6-03 

•1073E*02 

•63596-03 

•1073E»02 

•1450E-02 

•1073E*02 

.39256-03 

.1073E*02 

.30256-04 

•  107  36*02 

.92786-03 

.10736*02 

• 1905E-03 

.10736*02 

.26026-03 

.10736*02 

•9978E-04 

•  1 073E*02 

• 92786-03 

.1Q73E*02 

•2379E-03 

.10736*02 

•2978S-03 

•1073E*02 

.30256-04 

• 1073E*02 

.15336-03 

•  107  3E*02 

■6873E-02 

• 1073E*02 

•16996-03 

.1 07  36*02 

•1 167E-01 

.107  3E*02 

•  1772E-04 

• 1073E*02 

• 1184E-04 

.17736*02 

.98156-03 

•1073E*02 

•46836-04 

•1073E*02 

.77396-05 

• 1073E*02 

•7555E-05 

.63896*02 

•  1992E-02 

•63896*02 

•38216-02 

•63896*02 

• 1798E-03 

.63896*02 

•  1092E-0  3 

.63896*02 

.37386-01 

•6389E*02 

•27266-02 

.63096*02 

•1786E-02 

•63896*02 

• 1450E-01 

.63896*02 

•  451 9E-01 

•6389E*02 

• 3925E-02 

.63896*02 

.1092E-03 

•6389E*02 

• 2726E-02 

.63896*02 

•60806-04 

.63896*02 

•2695E-02 

•6389E*02 

•3t40E-03 

•6389E*02 

.13666-03 

•6789E*02 

.17986-03 

•63896*02 

•1092E-03 

• 6389E*02 

•1647E-02 

•63896*02 

•  7268E-03 

.63896*02 

.1328E-02 

.63896*02 

.12286-03 

•6389£*02 

.17986-03 

.63896*02 

.27266-02 

.12086*03 

•42316-02 

•14736*03 

•  4  265E-02 

• 1473E*03 

.42316-02 

•2227£*03 

•3275E*03 

•5741E-01 

•54356*00 

•2S29E*03 

• 2271E-01 

•2553E*03 

•1 992E-02 

•2553E*03 

•42316-02 

Table  4-4d.  Cumulative,  normalized  weight  arrays  for  artificial 
source  transfer-coefficients  and  emission  radiances 
in  sample  Problem  5. 


ITT,»T1  = 


0. 

0. 

•5640E-05 

• 5620E-05 

•8604E-04 

•8573E-04 

•8936E-04 

•8904E-04 

• 1 991E-03 

•1691E-03 

• 2024E-03 

•1725E-03 

•2635E-03 

•13846-02 

•2669E-03 

• 1 387E-02 

• 2725E-03 

•1393E-02 

.2966E-03 

• 1502E-02 

•3012E-03 

•3497E-02 

•7816E-03 

•4224E-02 

•8364E-03 

.60136-02 

.9598E-03 

.61226-02 

•970  2E-03 

•6259E-02 

• 9722E-03 

•6368E-02 

.1062E-02 

• 6429E-02 

•1120E-02 

.65526-02 

•1 246E-02 

•  7  881E-02 

•2462E-02 

•78876-02 

•2468E-02 

•7948E-02 

•2581E-02 

•7972E-02 

•4442E-02 

•8026E-02 

•4571E-02 

• 8S05E-02 

•4581E-02 

• 8509E-02 

•  47 10E-02 

• 8589E-02 

•4784E-02 

•8600E-02 

• 4897E-02 

•8689E-02 

•4971E-02 

•8691E-02 

•61 20E-02 

•8749E-02 

•6194E-02 

•8791E-02 

•6199E-02 

•8917E-02 

•6210E-02 

•8922E-02 

•6216E-02 

.8995E-02 

•6396E-02 

.91086-02 

•6577E-02 

.92366-02 

•  66916-02 

•9242E-02 

• 85516-02 

•9253E-02 

•1041E-01 

•9258E-02 

•1080E-01 

•9372E-02 

• 1 167E-01 

•9445E-02 

•11 85E-01 

•9558E-02 

•1 198E-01 

•  U41E-01 

•1210E-01 

•  114  2E-01 

.14696-01 

• l 143E-01 

•  15  00E-01 

•  1 156E-01 

•2280E-01 

•  1 174E-01 

•2538E-01 

•14296-01 

•2619E-01 

• 1685E-01 

.2636E-01 

•3839E-01 

•2852E-01 

•3847E-01 

•3247E-01 

•4103E-01 

•32526-01 

• 43596-01 

•3279E-01 

•4742F-01 

•3387E-01 

•4760E-01 

• 3388E-01 

•4778E-01 

• 34006-01 

•4790E-01 

.36186-01 

•4955E-01 

•3625E-01 

•5  070E-01 

•3626E-01 

•6245E-01 

•3767E-01 

•  9  382E-01 

•6412E-01 

•9776E-01 

.64146-01 

•9992E-01 

•6435E-01 

•14056*00 

•6448E-01 

• 1431 E*00 

•6449E-01 

•  143  4E*00 

•6457E-01 

•1 6976*00 

•6464E-01 

• 1699E*00 

•6467E-01 

• 17016*00 

•6469E-01 

.171 36*00 

•7648E-01 

•1734E*00 

•7858E-01 

•1736E*00 

.11016*00 

.17466*00 

•1317E*00 

.1764E*00 

•1 333E*00 

.17756*00 

•1372E*00 

•  1775E  *00 

•1397E*00 

.18536*00 

• 1399E >00 

•  1  853E  *00 

•1401E*00 

.18556*00 

•141 26^00 

•18596*00 

•141 8E*00 

•1877E*00 

• 141 9£*00 

•1 878E*00 

.141 8E*00 

.18796*00 

•  1 459E*00 

•18876*00 

•1459E*00 

• 1887E*00 

• 1460E*00 

.1 887E*00 

.  146 lfc*00 

• 1 887E*00 

. 14616*00 

.18876*00 

•146 1E*00 

.18906*00 

.14652*00 

•18916*00 

•1 475E*00 

•19056*00 

• 1485E*00 

• 1908E*00 

•148flE*00 

•1908E*00 

•1490E*00 

.19296*00 

.'l492E*00 

•1 930E*00 

•14926*00 

•1931E*00 

•1493E*00 

• 1940E*00 

•1494E*00 

• 1 9496*00 

•  149  4E*00 

• 1 9496*00 

•1495E*00 

• 1949E*00 

• 15136*00 

•1956E*00 

• 15156*00 

•1956E*00 

•15166*00 

•19596*00 

• 1523E*00 

•19996*00 

•1524E*00 

.20566*00 

•1S25E*00 

.20576*00 

•  1526E*00 

•2063E*00 

• 1546E*00 

•20696*00 

•1572E*00 

• 2084E*00 

.19796*00 

. 2088E*00 

.19916*00 

. 2088E*00 

•20166*00 

•2097E*00 

•2017E*00 

.2099E*00 

.2317E*00 

•2102E*00 

•20176*00 

•21036*00 

•2017E*00 

•  21 12E*00 

•2044E*00 

•  21156*00 

•20516*00 

.21176*00 

.20536*00 

•2118E*00 

•2055E*00 

•21 19E*00 

•  2061 E*00 

•  21 88E*00 

•21306*00 

•21906*00 

•2145E*00 

•2307E*00 

•2202E*00 

.23076*00 

•  2203E*  00 

• 2307E*00 

•2204E*00 

•23176*00 

.22126*00 

•2317E»00 

•2213E*00 

.23176*00 

•  221 3E*00 

•23176*00 

.22266*00 

.23372*00 

•2269E*00 

•2376E»00 

•  22876*00 

•2377E*00 

.22966*00 

.23796*00 

•24136*00 

•27536*00 

.24176*00 

.  2780E»00 

•2426E*00 

•2798E*00 

•2429E*00 

•29436*00 

• 2445E*00 

•  3  3965 ♦OO 

.24656*00 

• 3435E»00 

•2466E*00 

•34366*00 

• 2505E*00 

•34636*00 

•2505E*00 

•  3  46  46*00 

.25066*00 

.3  491 E^00 

•25086*00 

•3494E*00 

•2548E*00 

•3495E*0  0 

•2S49E*00 

•  34976*00 

.25516*00 

• 3498E»00 

•2578E*00 

•35156*00 

. 27236*00 

.3522E*Q0 

.3099E*00 

•35356*00 

•3 126E*00 

•  35376*00 

•  3 1 69E*00 

• 3539E*00 

•3197E*00 

.3566E*00 

•  32 1 OE*00 

•3608E*00 

•3664E*00 

. 3651E»00 

• 4241E*00 

•3693E*00 

•  44696*00 
• IOOOe+OI 

.42686*00 
• 10006*01 

.4S12E*00 

•4495E*00 

•9937E*00 

•45156*00 

•9957E*00 

•  4558E*  00 

STATISTICAL  CL0U0S,  P1,P2, 

P3,P4  = 

•  10 

.25 

•  50 

90 

SOURCE  TRAMSPER 

CP  T1,T2, 

T3#T4  = 

.4376-51  .1 

L54E-06  • 2806*00  .1526*01 

EMISSION 

£l,F2, 

E3,S4  = 

.S30E*01  .6396*02  .261E*03  .314E*03 

Table  4-5b.  Sorted  array  for  solar  source  transfer-coefficients 
(and  carried-along  weights)  for  sample  Problem  5. 


TBANS,UT  s 


0. 

0. 

• 6451E-03 

.43856*00 

.38876-02 

•1992E-02 

.38876-02 

.42316-02 

•43986-02 

.22716-01 

.10486-01 

•5741E-01 

.2714E-01 

.42656-02 

.27146-01 

.42316-02 

•  3171E-01 

.42316-02 

.50616-01 

.45196-01 

.50616-01 

.  3925E-02 

.50616-01 

.60806-04 

.50616-01 

.26956-02 

•5061E-01 

.31406-03 

.50616-01 

.27266-02 

.50616-01 

.17986-03 

.50616-01 

•1092E-03 

.50616-01 

.27266-02 

.50616-01 

•1798E-03 

.50616-01 

.19926-02 

.50616-01 

.12286-03 

•5061E-01 

.17986-03 

.5061E-01 

.10926-03 

.50616-01 

.37386-01 

.50616-01 

• 3821E-02 

•5061E-01 

.17866-02 

.506  1  6-01 

.14506-01 

.50616-01 

.16476-02 

.5061E-01 

•7268E-03 

.50616-01 

.13286-02 

•5061E-01 

•2726E-02 

.50616-01 

.10926-03 

•  5061 E- 01 

. 1 366E-03 

.72556-01 

•  9278E-03 

.72556-01 

.11676-01 

.72556-01 

.30256-04 

.72556-01 

.23796-03 

.72556-01 

.927  8E-03 

.72556-01 

.14506-02 

.7  255E-01 

.92786-03 

•7255E-01 

.6873E-02 

.72556-01 

• 3025E-04 

•7255E-01 

.26026-03 

.72556-01 

.26026-03 

•72556-01 

•  56  97E-02 

.72556-01 

.23146-04 

.72556-01 

.6  359E-0  3 

•7255E-01 

.63596-03 

•7255E-01 

• 4683E-04 

.72556-01 

.39256-03 

.7255E-01 

.30256-04 

•7255E-01 

.16996-03 

•7255E-01 

.23146-04 

•7255E-01 

.40506-02 

.7  25  5  6-01 

.63596-03 

.72556-01 

.58156-03 

•7255E-01 

. 11846-04 

.72556-01 

.29786-03 

•  7255E-01 

.17726-04 

.72556-01 

.23146-04 

•7255E-01 

.75556-05 

.72556-01 

.77396-05 

.72556-01 

.19056-03 

.72556-01 

.15336-03 

.72556-01 

.9978E-Q4 

.76216-01 

.17986-03 

.76216-01 

.25616-02 

.76216-01 

.25616-02 

•7621 E-01 

• 3821E-02 

.76216-01 

.18516-02 

.76216-01 

.31336-01 

•7621E-01 

.12286-03 

•  7621E-01 

.25616-02 

.76216-01 

.11446-02 

•7621E-01 

.17986-03 

.76216-01 

.31406-03 

•7621E-01 

•2632E-01 

.76216-01 

.16476-02 

•7621E-01 

.26236-03 

.76216-01 

.11746-01 

•7621E-01 

.21566-02 

.76216-01 

.12826-03 

.76216-01 

.21646-02 

•7621E-01 

.18516-02 

.76216-01 

.25736-02 

.76216-01 

.17986-03 

.76216-01 

.72986-04 

.76216-01 

.12286-03 

•7621E-01 

.12826-03 

.76216-01 

.11446-02 

.76216-01 

.86926-03 

.76216-01 

.39256-02 

•7621E-01 

.11326-03 

.76216-01 

.40536-01 

•7621E-01 

•56126-05 

.76216-01 

.12826-03 

•76216-01 

.56126-05 

.76216-01 

•  1 053E-04 

.76216-01 

.72986-04 

.76216-01 

.11326-03 

•7621E-01 

.72986-04 

.76216-01 

•  1 132E-0  3 

.76216-01 

•  185  IE- 02 

.76216-01 

.56126-05 

•  7621E-01 

.10536-04 

.76216-01 

.12556-03 

.76216-01 

.27856-03 

.76216-01 

.25476-02 

•7621E-01 

.20856-02 

.75216-01 

.21526-01 

.76216-01 

• 1 409E-02 

.76216-01 

.77606-02 

•7621E-01 

.10066-02 

.76216-01 

.12446-03 

.76216-01 

.10796-02 

.76216-01 

.86156-05 

•7621E-01 

•8335E-04 

•7621E-01 

.54466-04 

.76216-01 

.81146-03 

.76216-01 

.38656-03 

•7621E-01 

. 1053E-04 

.76216-01 

.36096-04 

.76216-01 

.69406-04 

.76216-01 

.20476-03 

•7621E-01 

.40036-05 

.76216-01 

.23896-04 

.7621E-01 

.10396-04 

• 7621E-01 

.16426-03 

.76216-01 

.71436-04 

.84086-01 

.10926-03 

.64086-01 

.6  0806-04 

.84086-01 

.72686-03 

.84086-01 

.17866-02 

.84086-01 

•  1 092E-03 

.8408E-01 

.13666-03 

.84086-01 

.60016-04 

.84086-01 

.60806-04 

.84086-01 

.12286-03 

.84086-01 

•1 3286-02 

.84086-01 

.10925-03 

•8408E-01 

.33086-05 

•84086-01 

.56126-05 

.84086-01 

.80016-04 

.84086-01 

. 3308E-05 

• 8408E-01 

. 8001 E-04 

.84086-01 

.12556-03 

.84086-01 

.12106-02 

.84086-01 

.33086-05 

• 8408E-01 

.56126-05 

•8408E-01 

.58266-04 

.84086-01 

.42336-04 

.84086-01 

.24006-04 

.84086-01 

•5446E-04 

. 8408E-01 

.47816-03 

.84086-01 

.40086-05 

.84086-01 

.45186-05 

.84086-01 

.84086-01 

.1039E-04 

.19926-02 

.84086-01 

•8903E-04 

.84086-01 

• 1 962E-05 

.84086-01 

.56126-05 

Table  4-6.  Definitions  of  quantities  in  Tables  4-4a  and  4-5a  . 


Column  Heading  Description 

I  IDX  Cloud-configuration  index 

2,3,4  LL,MM,HH  Cloud  indexes  for  low-,  medium-,  and 

high-layer  clouds 

5  TRANS  Transfer  coefficient  (for  artificial  source 

in  Table  4-4a  and  for  solar  source  in 
Table  4-5a) 

6  EMISS  Emission  radiance  from  highest  cloud  in 

set  (independent  of  type  of  source) 

7,8,9,10  CC=2,3,4,5  Weights  or  relative  occurrence  frequencies 

of  given  cloud  set  with  cloud-coverage 
category  CC  (or  cloud-coverage  fractions 
of  0.3,  0.5,  0.8,  and  1.0) 

II  WT  Total  weight  or  relative  occurrence 

frequency  of  given  cloud  set  (sum  of  weights 
for  CC=2 ,3,4,5) 


12 


WT( SUM ) 


Cumulative  sum  of  WT  for  all  cloud- 
configuration  sets 


(Appendix  Equation  (30a))  exceeds  100.  In  this  event  function  TRANSF  uses 
the  initialized  value  of  FLUX,  l.E-50,  an  arbitrarily  small  value.  The 
results  for  IDX=160  are  provided  by  subroutine  SCLOUD  for  the  Earth's 
surface.  The  initial  setting  of  WT(160),  the  probability  of  seeing  the 
ground  with  or  without  clouds  being  present,  occurs  in  subroutine  CLDWT, 
and  the  corresponding  WT(SUM)  would  be  exactly  1.0  except  for  small 
numerical  inconsistencies.  WT(160)  is  reset  in  subroutine  SCLOUD  to  take 
account  of  the  two-leg  CFLOS  required  to  see  tha  ground  lighted  by  the 
source  (either  artificial  or  the  sun). 

Before  sorting  the  TRANS  and  EMISS  arrays,  we  have  augmented 
each  of  them  (as  well  as  the  WT  and  WT1  arrays)  with  one  additional  member 
whose  value  is  set  to  zero.  This  procedure  allows  subroutine  LINEAR  to 
interpolate  within  the  array  given  to  it  when  the  weight  of  the  smallest 
(nonzero)  member  exceeds  the  smallest  fractile  (now  0.10)  for  which  an 
integral-distribution  value  is  desired. 

Tables  4-4b  and  4-5b  give  the  TRANS  array  after  it  has  been 
sorted  (by  a  call  to  subroutine  SORTLJ  from  subroutine  SCLOUD)  in  the  order 
of  increasing  values.  The  WT  array  is  simply  carried  along.  The  TRANS 
and  WT  values  are  presented  in  pairs  in  Tables  4-4b  and  4-5b. 

Table  4-4c  gives  the  EMISS  array  after  it  has  been  sorted  in  the 
order  of  increasing  values.  The  WT1  array  is  simply  carried  along.  The 
EMISS  and  WT1  values  are  presented  in  pairs  in  Table  4-c. 

Tables  4-4d  and  4-5c  give  the  cumulative,  normalized  weight 
arrays  WT  and  WT1 ;  these  quantities  are  also  presented  in  pairs. 

At  the  end  of  Tables  4-4d  and  4-5c  the  values  of  the  transfer 
coefficient  and  emission  radiance  are  given  corresponding  to  the  four 
user-selected  fractiles  (here  taken  to  be  0.10,  0.25,  0.50,  and  0.90) 
determined  from  the  cumulative,  normalized  weight  arrays  WT  (for  transfer 
coefficient)  and  WT1  (for  emission  radiance). 

We  repeat,  as  noted  in  Section  2.2,  that  in  the  160-member  dis¬ 
tributions  for  the  solar  source  transfer-coefficients  and  emission 
radiances,  there  are  only  10  different  values  (nine  for  clouds,  one  for  the 
Earth)  because  (a)  there  are  only  nine  different  cloud-top  altitudes 
(cf.  Table  2-1),  (b)  the  bidirectional  reflectance  of  any  cloud  is  assumed 
to  be  approximated  by  that  for  cloud-type  4  (stratocumulous) ,  and  (c)  the 
directional  emissivity  computed  for  cloud-type  4  is  used  for  all  cloud 
types.  We  have  collated  the  essential  results  for  the  solar  source 
transfer-coefficient  array  in  Table  4-7a  and  those  for  the  emission  radiance 
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array  in  Table  4-7b.  These  latter  two  tables  should  facilitate  under¬ 
standing  the  values  of  the  solar  transfer-coefficient  and  emission 
radiance  for  the  fractiles  0.10,  0.25,  0.50,  and  0.90  printed  at  the  end 
of  Table  4-5c. 

The  variation  in  values  of  the  (attenuated)  solar  source 
transfer-coefficient  (TCOEF)  for  the  nine  cloud  tops  is  due  only  to  the 
variation  in  the  (two-leg)  air  transmittance  between  the  different  cloud- 
top  altitudes  and  12-km  altitude,  as  verified  by  our  including  in  Table 
4-7a  the  cloud-top  altitude  and  the  air  transmittance  from  the  cloud-top 
altitude  to  12-km  altitude.  The  variation  in  values  of  the  (attenuated) 
emission  radiance  (EMISS)  for  the  nine  cloud  tops  is  due  to  the  combined 
effects  of  (a)  the  (one-leg)  air  transmittances  from  the  different  cloud- 
top  altitudes  to  12-km  altitude  and  (b)  the  temperature  entering  the 
Planck  emission  function. 

In  the  distributions  presented  in  Tables  4-7a  and  4-7b  the 
importance  of  the  Earth's  surface  is  apparent  since  it  has  by  far  the 
largest  single  weight  associated  with  it.  This  feature  occurs  not  only  for 
the  data  in  KM0DEL=1  but  for  the  data  given  in  BLOCK  DATA  for  each  of  the 
10  KMODELs.  To  illustrate  this  general  feature  we  have  evaluated  the 
probability  (Pgnd)  for  a  detector,  pointed  straight  down,  of  seeing  the 
ground  either  with  or  without  clouds  being  present : 

Pgnd  =  CCOVER(l.K)  + 

CC0VER(2,K)  x  CFL0SF(4,e)  + 

CC0VER( 3 , K)  x  CFL0SF(6 , 6 )  + 

CC0VER( 4 , K)  x  CFL0SF( 9 , 8 )  + 

CC0VER( 5 , K)  x  CLFOSF(ll,e) 

where  K  =  KM0DEL  and  6=0  degrees . 

The  results  are: 


KMODEL 

Pgnd 

KMODEL 

P  , 
gnd 

1 

0.543 

6 

0.611 

2 

0.543 

7 

0.666 

3 

0.484 

8 

0.591 

4 

0.633 

9 

0.628 

5 

0.513 

10 

0.559 
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It  is  instructive  to  relate  some  of  the  Problem-5  results  to 


those  for  Problems  1  to  3.  The  cloud  configurations  for  IDX  equal  to 
4,  8,  and  13  in  Table  4-4a  are  the  same  as  those  for  Problems  1,  2,  and  3, 
respectively.  Since  the  source  is  below  all  clouds  in  both  cases,  we 
should  expect  the  (transmission)  transfer  coefficients  in  Table  4-4a  to 
correspond  to  those  at  zero  lateral  displacement  in  Problems  1,  2,  and  3, 
except  that  air  attenuation  is  included  in  Problem  5  but  not  in  Problems 
1,  2,  and  3.  However,  the  air-attenuation  path  is  the  same  (cf.  Figure  2-6b) 
for  all  cloud  configurations  in  Problem  5.  Thus,  we  expect  that  the 
TRANS  results  given  in  Table  4-4a  for  IDX  =4,  8,  and  13  will  differ  by 
a  constant  factor  from  those  for  zero  lateral  displacement  in  Problems 
1  to  3.  To  predict  this  factor  we  write: 

TRANS (Pm  5)  =  TC'(Pms  1,2,3)  x  TAIR(2 . 5 , 0 . , 0 . 01 )  x  (SDISTO)2. 


By  using  the  interpolation-extrapolation  scheme  for  air  transmittance  as 
described  in  Section  3-2.28,  we  hand-calculate  the  air  transmittance 
to  be 

TAIR  =  ^  (T1/T2)_(h_hl)/(h2"hl) 

=  (0.2150)  (0. 2150/0. 3247)_(0,01_1)/(2_1) 

=  0.1430 

and 

(SDISTO)2  =  (11. 99)2  . 

Thus  we  predict 

TRANS  =  20.55  TC' 

It  is  satisfying  that  our  prediction  is  correct,  as  seen  from  the  data 
in  Table  4-8. 

The  reason  we  did  not  include  IDX=114  and  Problem  4  in  this 
comparison  is  that  the  integral  to  be  evaluated  for  TRANS  in  Table  4-4a 
was  skipped  owing  to  the  large  argument  of  the  sinh  and  cosh  functions 
and  the  concomitant  smallness  of  the  value  of  TRANS.  (See  the  related 
comment  made  previously  with  respect  to  Table  4-4a. ) 


120 


Table  4-8.  Comparison  of  transfer  coefficients  from  Problem  5  with 
modified  transfer  coefficients  from  Problems  1  to  3. 


IDX 

TRANS a 

-1 

sr 

(TC' )b 
sr"^  km-2 

TRANS /TC’ 
km2 

4 

0.4617"1 

0.2247-2 

20.55 

8 

0.6812-2 

0.3315-3 

20.55 

13 

0.2199"7 

0.1070-8 

20.55 

aFrom 

Table  4-4a. 

. 

bFrom 

Table  4-2  for 

RANGE=0 . 0 . 

4-3.3  Problem  6:  Use  of  the  deterministic  cloud  submodel  for  an 
ellipsoidal  cloud 

Problem  6  illustrates  the  selection  by  subroutine  LOS  of  lines- 

of-sight  and  the  evaluation  by  subroutine  CL0UD3  of  transfer  coefficients 

and  emission  radiances  for  these  same  lines-of-sight . 

For  this  problem  we  have  placed  a  single  ellipsoidal  cloud  at 

10-km  altitude.  The  lengths  of  the  semiaxes  are  a=4  km,  b=2  km,  and 

c=l  km.  The  cloud's  properties  are  those  for  the  Type-4  cloud  (strato- 

cumulous;  see  Table  2-1).  The  artificial  source  is  also  placed  at  10-km 

altitude  and  10  km  to  the  east  of  the  cloud  center.  The  solar  source 

has  a  zenith  angle  of  45  degrees.  The  detector  is  at  synchronous 

altitude,  directly  above  the  cloud.  Wavelength  is  2.5  pm.  We  have  made 

o 

two  runs,  one  (Problem  6a)  with  FOV  =  1  km  and  the  other  (Problem  6b) 

2 

with  FOV  =  20  km  .  The  input  deck  is  shown  in  Figure  4-5. 

Figure  4-6a  is  an  outline  of  the  cloud  in  a  horizontal  plane 

through  its  center,  The  64  different  x,y  coordinates  of  the  facet  centers 

are  shown  by  small  dots.  The  16  x,y  coordinates  of  the  facet  centers 

that  are  considered  in  subroutine  LOS  by  incrementing  the  facet-center 

indexes  by  two's  are  denoted  by  crosses. 

Figure  4-6b  pertains  to  the  cloud  for  the  artificial  source. 

We  show  the  16  above-defined  facet  centers  visible  to  the  detector, 

marked  by  crosses.  Of  these,  the  three  that  are  also  visible  to  the 

source  are  marked  as  RI(R  denotes  reflective,  I  the  sequence  number 

2 

(1  to  3)).  For  FOV  =  1  km  ,  each  of  these  three  reflective  points 


PROBLEM  6A 


HD 

THETAD 

PHID 

fov 

DETECTJR 

35793. 

90.00001 

-90.0001 

1.0 

HSORCE(l) 

THETAS(l) 

PH1S(1) 

MKT 

SOURCE 

10. 

1SUN 

1 

ALAM 

2.50 

NODE 

0 

NCLOUD 

1 

KCLOUD 

90.000001 

-89.90 

1000. 

DETEKMI NX  STIC 

NO.  OF  CLOUDS 

CLOUD  PROPERTIES 

HCLOUD 

THETCD 

PHICU  CLOUDA  CLOUDB 

CLOUDC  CLDPSI 

4  10. 

90. 

-90. 

4.0  2.0 

1.0  0. 

CARD  1 
CARD  2 
CARD  3 


PROBLEM  6B  1 

HD 

THETAD 

PHID 

POV 

DETECTOR 

35793. 

90.00001 

-90.0001 

20. 

HSOHCE(l) 

THETAS ( 1 ) 

PUS(l) 

WKT 

SOURCE 

10. 

90.000001 

-89.90 

1000. 

ISDN 

ALAM 

WAVELENGTH 

2.50 

MODE 

DETERMINISTIC 

CARD 

1 

NCLOUD 

NO.  OF  CLOUDS 

CARD 

2 

KCLOUD 

CLOUD  PROPERTIES 

CARO 

3 

HCLOUD 

THETCD 

PHICD  CLOUDA  CLOUDB 

CLOUDC  CLDPSI 

4  10. 

90. 

-90. 

4.0  2.0 

1.0  0. 

Figure  4-5.  Input  data  for  Problems  6a  and  6b. 


Figure  4-6a.  Problem  6  cloud.  Shown  is  the  horizontal  plane  through  the  center  of 
the  ellipsoidal  cloud  at  10-km  altitude.  Projected  positions  of  all 
facet  centers  are  indicated  by  dots;  those  examined  in  subroutine  LOS 
for  visibility  are  marked  with  crosses. 


the  center  of  the  ellipsoidal  cloud  at  10-km  altitude.  Source  is  at 
coaltitude  and  10-km  eastward.  Subroutine  LOS,  for  F0V=1  km2,  selects 
R1 ,  R2 ,  and  R3  as  the  reflective  facet  centers  and  T1  through  T8  as  the 
transmittive  facet  centers.  For  FOV  =  20  km2 ,  R2  and  R3  are  rejected 
and  all  others  remain. 


survives  the  selection  process  since  none  lines  within  0.7  /FOV  of 

o 

another.  For  FOV  =  20  km  ,  R2  and  R3  do  not  survive  the  selection 

process.  The  remaining  13  points  (16  minus  3)  are  transmittive ;  of 

these,  eight  (the  maximum  number  of  transmittive  points  retained  by  the 

code)  are  marked  as  TI  (T  denotes  transmittive,  I  the  sequence  number 

(1  to  8).  These  same  eight  points  obtain  for  both  values  of  FOV. 

Figures  4-6c  and  4-6d  pertain  to  the  cloud  for  the  solar  source 

2  2 

and  correspond  to  FOV  being  1  km  and  20  km  ,  respectively.  Again  we  show 

the  16  facet  centers  visible  to  the  detector,  marked  by  crosses.  All 

of  these  facet  centers  are  visible  to  the  sun;  there  are  no  transmittive 

2 

points.  For  FOV  =  1  km  ,  11  of  the  points  (marked  as  R1  through  Rll  in 

2 

Figure  4-6c)  survive  the  selection  process;  for  FOV  =  20  km  ,  four  of 
the  points  (marked  as  R1  through  R4  in  Figure  4-6d)  survive  the  selection 
process. 

The  transfer  coefficients,  emission  radiances,  and  cloud 
coordinates  of  the  facet  centers  for  both  sources  and  FOV’s  are  given  in 
Table  4-9,  to  which  a  guide  is  provided  in  Table  10.  The  coordinates  of 
the  first  reflective  point  and  the  eight  transmittive  points  are,  of 
course,  the  same  in  Parts  (a)  and  (b). 


the  center  of  the  ellipsoidal  cloud  at  10-km  altitude.  Source  is  th 
sun,  eastward  at  45-deg  zenith  angle.  Subroutine  LOS,  for  F0V=1  km^ 
selects  R1  through  R] 1  as  the  reflective  facet  centers.  There  are  n 
transmittive  facet  centers. 


Problem  6b  with  solar  source.  Shown  is  the  horizontal  plane  through 
the  center  of  the  ellipsoidal  cloud  at  10-km  altitude.  Source  is  the 
sun,  eastward  at  45-deg  zenith  angle.  Subroutine  LOS,  for  FOV=20  km^ 
selects  R1  through  R4  as  the  reflective  facet  centers.  There  are  no 
transmittive  facet  centers. 


"i 

i 

Table  4-9. 

Results 

for  Problem  6.  Transfer  coefficients  (sr’1) , 

emission 

radiances 

(W/(km^ 

'  sr  pm) ) 

,  and  cloud  coordinates 

(km).  See  Table  4- 

-10  for 

guide . 

TCOEP 

EHISS 

ICC 

tree 

zee 

•3483E-01 

•2492E*01 

.36986*01 

•  2633E*  00 

.3577E*00 

•  1 844E-0 1 

• 15736*01 

• 3314E*01 

-• 90  26E*00 

•33146*00 

•5076E-02 

.11176*01 

.22605*01 

•15636*01 

•  26515*00 

(a) 

.4383E-04 

•4348E*01 

-.7572E*00 

.1923E*01 

• 1994E*00 

\ 

•  6605E-0  3 

• 4352E*0l 

•7572E*00 

-•  1923E*01 

•19945*00 

• 5299E-06 

• 1 11 8R*01 

-• 2260E*0l 

-.  1563E*01 

•26516*00 

• 1079E-06 

•1573E*01 

-•  331 4E*  01 

•  90  26E*0Q 

•  3314E*0  0 

.89405-07 

• 24935*01 

-.3698E*01 

2C  33  E*00 

•  3577  E*00 

•  1 1 51E-03 

• 2635E*01 

-.20806*00 

.5282E*00 

• 9t>3lE*00 

• 2615E-03 

• 2635E  *01 

•2080E*00 

-•5282E*00 

•  9631 E*00 

•3276E-03 

•2603E*01 

•  4747  E*00 

• 3283E*00 

•9793E*00 

TCOEP 

EHISS 

XCC 

YCC 

zee 

•  46  26  E- 0  2 

• 1501E*0l 

•  369BE*  01 

.  26  336*00 

.35776*00 

•7286E-03 

•3561E*01 

-.7572E*00 

• 1923E*01 

•1994E*00 

•  377  IE- 02 

.35652*01 

.75726*00 

-• 1923E*0l 

.19945*00 

(b) 

• 1054E-03 

•3086E*0l 

-.2260E*01 

-• 1S63E*01 

•26515*00 

•  17  22E-04 

.23006*01 

-.3314E*01 

•  9Q26E*  00 

•  3314  E*00 

• 1650E-05 

•1502E*0l 

-• 36986*01 

-• 26336*00 

•  3577  E*00 

• 2307E-02 

•2552E*01 

-• 2080E*00 

. 52  d2E*00 

• 9631E*00 

•  3846E-0  2 

• 25525*01 

• 2080E*00 

-.52826*00 

•  96  31 E*00 

•5055E-02 

•  2515E  >0 1 

• 4747E*00 

.3283E*00 

•  9793  E*00 

TCOEP 

EHISS 

XCC 

YCC 

zee 

• 4833E-0 1 

• 2492E*01 

•  3«j98E*01 

•  2t>  336*00 

.35775*00 

•2992E-01 

•1573E*0l 

•  331 4E*01 

-• 9026E*00 

.33145*00 

• 1930E-01 

■  1 1 1 7E*0l 

.2260E*0l 

•  1563E*01 

•  2651 E*00 

•51 25E-01 

•43S2E*0l 

.7S72E*00 

-.  1923E*01 

.19945*00 

•6405E-01 

• 2581E*01 

• 5816E*00 

•4140E-01 

. 9892  E*00 

•5751E-01 

.  26  35£*  01 

2080E*  00 

• 52826*00 

•  96  31 E*00 

r 

(c) 

•5927E-01 

•2603E*01 

-. 4747  E*00 

-• 3263E*00 

•9793E*00 

•4914E-01 

.43485*01 

-.7572E*00 

.1923E*01 

.19945*00 

• 5535E-02 

•  11 185*01 

-• 2260E*01 

-.1563E*01 

.26515*00 

•6696E-02 

.15735*01 

-• 3314E*01 

•  90  26E*00 

.33146*00 

1 

• 8916E-02 

•2493E*01 

-• 3698E*0l 

-• 26336*00 

•  3577  E*00 

I  1 

TCOEP 

EHISS 

XCC 

YCC 

zee 

(d) 

• 2768E-0 1 

.15016*01 

•3698E*01 

.2633E*00 

•  3577  E*00 

• 4876E-01 

.35652*01 

• 75726*00 

-• 19236*01 

.19946*00 

• 4263E-01 

•  3561E  *01 

-.7572E*00 

.1923E*01 

•  1994  E*00 

• 1612E-0 1 

•1502E*01 

-.3698E*01 

-.26336*00 

•  3577  E*00 

/ 

/ 


Table  4-10. 

Guide  to 

Table  4-9  results. 

Table  4-9 

FOV 

km^ 

Reflective 

Transmittive 

Part 

Source 

Points 

Points 

(a) 

Artificial 

1 

3 

8 

(b) 

Artificial 

20 

1 

8 

(c) 

Solar 

1 

11 

o 

(d) 

Solar 

20 

4 

0 
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APPENDIX  A 

SOME  DIFFUSION  PROBLEM  SOLUTIONS 


A-l  INTRODUCTION 

The  problem  of  determining  the  light  flux  out  of  clouds  in 
various  configurations  of  clouds  and  sources  recurs  frequently.  For 
convenience  the  results  for  a  number  of  such  problems  are  collected  here. 
No  claim  is  made  that  these  are  novel  or  original,  many  of  the  solutions 
being  well-known  or  using  standard  methods.  A  unified  compilation  with 
a  consistent  terminology  may  nevertheless  be  useful. 

A  rigorous  treatment  using  the  transport  equation  or  Monte  Carlo 
methods  to  evaluate  the  effects  of  scattering  and  absorption  would  be, 
at  best,  extremely  time  consuming.  Additionally,  the  specification  of 
the  geometry  of  the  clouds  only  approximates  real  cloud  configurations, 
making  such  a  careful  evaluation  of  flux  overly  elaborate.  Therefore  the 
approach  taken  here  has  been  to  solve  the  problem  in  the  approximation 
of  the  diffusion  equation.  The  basic  assumptions  are  that  there  be  many 
scatterings  and  that  precise  surface  effects  be  unimportant.  An  improve¬ 
ment  over  the  most  elementary  diffusion  equation  solutions  can  be  achieved 
by  using  scattering  parameters  based  on  rigorous  transport  theory  or  by 
using  the  results  of  Monte  Carlo  calculations  to  adjust  the  results. 

A-l.l  Basic  Equations 

The  development  follows  the  lines  indicated  by  Weinberg  and 
Wigner^1^  for  neutrons.  The  basic  equation  to  be  solved  is  the  dif¬ 
fusion  equation: 

D  div  grad  <t>-Ia<J>  +  S  =  i  ||  (1) 

2 

'"here:  $  is  the  light  flux,  photons/ (cm  sec) 

D  is  the  diffusion  coefficient 

is  the  macroscopic  absorption  cross-section 

S  is  the  light  source  per  unit  volume  per  unit  time, 
photons/(cm3  sec) 

v  is  the  speed  of  light 

The  required  boundary  conditions  are  given  in  terms  of  the  flux: 
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J±  -  *♦  *  »  (*  II  -  T!  It)  <2> 

J+  is  the  flux  in  the  +  Z  direction. 

The  constants,  a  and  b,  are  left  arbitrary.  From  the  usual  definition  of 
flux,  b  =  D/2.  In  "standard"  usage,  a  =  Eg/4 Z,  D  =  Zs/3E£tr,  and  = 

£-  where  Z,  I  ,  and  Z^r  are  the  macroscopic  total,  scattering,  and 

transport  cross-sections,  respectively,  and  jj  is  the  average  cosine  of 
the  scattering  angle.  Here  a  and  b  are  left  undermined  to  allow  later 
evaluation  using  better  values  based  on  suitable  exact  calculations.  In 
general  the  physical  boundary  conditions  require  a  specified  inward  flux, 
sometimes  zero.  The  frequently  used  "extrapolated  end-point”  formulation 
relates  0  and  9<p/dZ  on  the  surface.  The  desired  quantities  are  the  emerg¬ 
ing  fluxes,  the  actual  flux  distribution  inside  the  medium  being  of,  at 
most ,  peripheral  interest . 

A-1.2  Diffusion  Parameters  for  an  Absorbing  medium 

In  Section  A-l.l  the  "standard"  evaluation  was  given  for  the 

diffusion  parameters  in  terms  of  the  macroscopic  cross-sections.  That 

evaluation  is  strictly  applicable  only  in  the  limit  of  small  absorption. 

A  more  generally  applicable  evaluation,  which  we  adopt  for  numerical 

results  in  this  report,  is  described  in  this  section. 

(2) 

Case,  de  Hoffmann,  and  °laczekv  '  provide  material  suitable  for 

a  better  evaluation  based  on  exact  solutions  of  some  problems.  The  results 
(2) 

of  Case  et  al .  are  derived  for  isotropic  scattering.  In  the  present 

application,  some  modification  has  been  made  to  account,  at  least  partially 

for  anisotropic  scattering.  A  judicious  utilization  of  the  transport  cross 

section  (which  introduces  the  anisotropy)  must  be  made. 

(21 

Case  et  al .  use  several  parameters  to  which  we  must  refer: 

c  =  average  number  of  emergent  particles  per 
collision  (cf.  Reference  2,  p.  43) 

=  Zs/Z  (3) 

jo  for  pure  absorption 
\  1  for  pure  scattering 
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1/k  (c)  =  diffusion  length  (dimensionless) 

°  (cf.  Reference  2,  p.  55) 

X(c)  =  linear  extrapolation  length  (dimensionless) 

(cf.  Reference  2,  p.  13-7) 

The  parameter  kq(c),  determined  by  c  according  to  Equation  (68)  on 
p.  55  of  Reference  2, 

c  =  Ko/tanh-1Ko  ,  (4) 

is  tabulated  and  plotted  in  Reference  2  (pp.  60,  63;  Figure  16  on  p.  56). 
The  parameter  X(c)  is  also  tabulated  and  plotted  in  Reference  2  (p.  136; 
Figure  35  on  p.  138).  Tabular  values  of  kq(c)  and  A(c)  are  given  here 
in  Table  A-l.  (In  the  computer  code,  we  use  more  closely-spaced  values 

Of  Kq.) 

Since  kq(c)  and  A(c)  are  dimensionless,  we  are  free  to  associate 

scale  lengths  with  them  in  some  reasonable  way.  We  do  so  by  choosing 

a  suitable  effective  macroscopic  cross-section  (equivalent  to  choosing  an 

effective  mean  free  path).  It  is  by  introducing  such  a  scale  length  that 

(2) 

we  adapt  the  results  of  Case  et  al.  to  accommodate  non-isotropic 
scattering  through  use  of  the  transport  cross-section.  We  first  discuss 
the  scale  length  used  with  kq(c)  and  then  that  with  X(c). 

We  assume  that  the  dimensioned  (1/km)  diffusion  length,  K,  and 
kq(c)  are  related  by  the  expression 

K  =  kq(c)  a  (5) 

where  a  is  an  effective  macroscopic  cross-section  (or  inverse  mean  free 
path)  whose  precise  form  will  be  established  below.  We  determine  a  by 
the  requirement  that  the  diffusion  coefficient,  D,  whose  basic  definition 
(Reference  1,  pp.  241,  242)  is 

D  =  £a/K2  ,  (6) 

approach  its  anisotropic-scattering  asymptotic  value  of  (32tr)-1  (cf. 
Reference  1  pp.  196,  197)  in  the  limit  that  2  approaches  zero.  To 

cl 

determine  a,  we  first  introduce  the  usual  relation  between  the  total, 
scattering,  and  absorption  cross-sections, 
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Table  A-l. 

Values  of 

,  .  a 

k  and  X  . 
o 

c 

X 

0.0 

1.000000 

1.0000 

0.1 

1 . 000000 

1 . 0000 

0.2 

0.999909 

0.9993 

0.3 

0.997414 

0.9889 

0.4 

0.985624 

0.9606 

0.5 

0.957504 

0.9201 

0.6 

0.907332 

0.8750 

0.7 

0.828635 

0.8300 

0.8 

0.710412 

0.7871 

0.9 

0.525430 

0.7472 

0.92 

0.474002 

0.94 

0.413976 

0.96 

0.340829 

0.98 

0.242983 

0.99 

0.172511 

1.00 

0.000000 

0.7104 

aFrom  Tables  8  and  23  of 

Case  et  al.^2^ 

Ea  =  1  "  >  (7) 

and  then  use  Equations  (3),  (5),  (6),  and  (7)  to  obtain 


(l-c)E 


Next,  we  use  the  fact  (Reference  2,  p.  58)  that 


■+  3  ( 1-c)  as  c+1 


(8) 


(9) 
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Thus,  in  order  that  we  can  have 


D  —  as  c+1 
tr 


we  must  have 


Z  _  1 

2  Z^. 


o  =  ZZ„ 


By  collating  our  results  thus  far,  we  have,  first  of  all,  the 
general  relation  from  the  usual  definition  of  flux, 

b  =  D/2  .  (12 

Then,  by  combining  Equations  (8)  and  (lib),  we  have 


k  2  Z 
o  tr 


and  from  Equations  (5)  and  (lib)  we  have 


K  =  K 


We  shall  now  relate  a,  given  by  Equation  (lib),  and  the 
parameter  A(c).  But  first,  we  note  that  the  use  of  the  linear  extra¬ 
polation  length  A(c)  implies  that  the  relation  between  <p  and  the  normal 
component  of  the  gradient  of  <f>  shall  be  that  of  the  result  for  a  semi¬ 
infinite  medium,  where  the  value  of  A(c)  is  determined  by  relating  the 
diffusion  solution  to  the  exact  transport  solution. 


length , 


We  now  assume 
d,  and  A(c)  are 


that  the  dimensioned  (km)  linear  extrapolation 
related  by  the  expression 


d  =  A(c)/cr 


(15) 


From  the  definition  of  the  linear  extrapolation  length,  d 
(cf.  Reference  1,  pp.  199,  261),  and  the  relation  for  the  photon  current, 
we  must  have 

d  =  b/a  .  (16) 

Thus,  from  Equations  (16),  (15),  and  (5)  we  have 


bo  _  bK 
a  A  (c )  A(c)ico(c)  ' 

To  summarize  our  parameter  evaluation,  we  have: 


(17) 


1. 

2. 

3. 

4. 

5. 

6. 

An  examination  of  the  values  for  <Q(c)  and  A(c)  in  Table  A-l 

shows  that 

~  =  <0(c)A(c)  11  .  (18) 

Thus,  the  factor  [(a/bK)  -  bK/a] ,  which  occurs  in  the  diffusion  solutions 
to  be  developed  later,  has  the  very  important  feature  of  never  being 
negative . 


c  =  Eg/E,  dimensionless,  Equation  (3) 

k  (c),  A(c)  from  tables 
o 

1-c 

D  =  — 5 -  km,  Equation  (13) 

Ko  <c>2tr 

K  =  k  (c)  /TT—  km-1,  Equation  (14) 
o  « 

b  =  D/2,  km,  Equation  (12) 

a  =  K  (c)x(c)-  ’  dimens:i-on]-ess  >  Equation  (17) 
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The  diffusion  parameters  for  both  the  "standard"  usage  and 
our  adopted  usage,  as  well  as  their  values  in  the  limiting  cases  c  =  0,1, 
are  given  in  Table  A-2.  That  the  limiting  values  are  valid  for  c  -*•  0 
(pure  absorption)  can  be  seen  by  considering  the  pure  absorption  solution 
for  a  semi-infinite  medium,  exp(  -  Z  Z).  The  positive-direction  photon 

Si 

current  is  just  exp(-£  Z),  the  negative-direction  current  is  0.  For 

Si 

c  ->  1  (pure  scattering)  ,  note  that  for  the  parameter  a,  our  limiting 
value  differs  somewhat  from  the  usual  one. 

In  the  code,  the  variables  a,b,  and  K  are  known  as  DALFA,  DBETA, 
an d  DKAPPA ,  respectively. 

A-2  SINGLE-LAYER  PLANAR  SLABS 

The  first  problem  to  be  considered  is  diffusion  in  a  single 
layer  of  finite  thickness.  The  sequence  of  problems  is  from  the  simplest 
to  the  more  complicated.  Time-independent  problems  are  treated  before 
time-dependent  ones,  one-dimensional  problems  before  two-dimensional 
ones.  This  sequence  has  value  both  from  the  standpoint  of  gradually 
increasing  mathematical  complexity  and  from  the  fact  that  the  simpler 
solutions  can  serve  as  limiting-case  checks  on  the  more  complicated 
solutions . 

A-2.1  Steady-State  Problems 

The  first  category  of  problems  is  steady  state,  hence  time- 
independent  .  They  are  of  intrinsic  value  and  also  serve  as  checks  on 
later  time-dependent  problems. 

A-2. 1.1  Plane  Source,  Surface  Deposition.  A  slab  of  thickness  L  has 
unit  flux  per  unit  area  incident  on  the  lower  surface  (Z  =  0),  no  flux 
incident  on  the  upper  surface.  The  transmitted  and  reflected  fluxes 
are  sought.  This  yields  the  equation  for  <j> : 

32<KZ)/8Z2  -  K2  <f>(Z)  =  0  (19) 

2  2  2 
with  K  =  Z  ID  for  the  standard  usage  and  K  =  K  (c)Z£.  for  our 
a  o  t  r 

adopted  usage.  The  boundary  conditions  are: 

J+(0)  =  1  (20a) 


J  (L)  =  0 


(20b) 


Table  A-2.  Summary  of  Diffusion  Parameters. 


The  reflected  flux  is  J_(0),  the  transmitted  flux  is  J+(L).  Solution 
of  this  problem  by  elementary  methods  yields  expressions  for  these 
quantities. 


1 


J+(L)  =  (cosh  KL  +  i  (^  +  x)  sinh  KL)_1  (21a) 

J_(0)  =  *  (fS  "  ?)  sinh  KL  J+(L)  •  (21b) 

There  are  several  comments  to  be  made  about  this  solution. 

The  condition  of  small  absorption  shows  up  here  immediately.  For  J  (0) 
to  be  a  positive  quantity  (as  it  should)  it  is  necessary  that  K  <a/b. 
This  condition  is  met  for  the  "standard"-usage  parameters  if  the 
absorption  is  small.  For  our  adopted  parameters  (see  Section  A-1.2), 
this  condition  is  always  met.  Secondly,  in  the  limit  of  no  absorption 
( K  ■+■  0 ) ,  the  sum  of  the  transmitted  and  reflected  fluxes  is  unity,  as  it 
should  be.  A  small  effort,  and  using  b=D/2,  gives  the  result  that 

J  (L)  +  J  (0)  +  f  Z  (|>(Z)dZ  =  1  . 

Jo  a 


—KL 

Finally,  it  can  easily  be  shown  that  J+(L)  <e  ,  the  defect  being  due 
to  increased  loss  in  multiple  collisions. 


A-2.1.2  Plane  Source,  Surface  Deposition  (alternate  treatment).  In 
Section  A-2.1.1  the  problem  of  flow  through  a  slab  was  treated  purely 
in  the  diffusion  approximation.  An  alternate  approach,  based  on  physical 
considerations,  separates  the  flux  into  two  parts,  that  which  has  had  no 
collisions  and  the  remainder.  This  leads  to  an  equation  for  the 
scattered  flux  with  a  source  term  coming  from  first  collisions.  The 
equation  to  be  solved  is 


3  2d> 
3Z 


Za<t>  +  E_  e~EZ 
a  s 


0 


(22) 


The  source  term  reflects  the  exponential  decrease  in  the  unscattered  flux 
and  the  fact  that  only  the  scattering  (not  the  absorption)  contributes 
to  the  source.  If  J+  is  now  considered  as  applying  to  only  the  scattered 
flux,  the  appropriate  boundary  conditions  are  that  there  be  no  flux 
incident  on  the  slab  surfaces: 
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(23) 


J+(0)  =  J_(L)  =  0 

The  solution  of  this  system  follows  the  usual  lines.  A  parti¬ 
cular  solution  of  the  inhomogeneous  equation  is 

v  C  e~EZ,  C  =  £s  D_1(K2-Z2)-1 

o 

where  K  is  given  as  before.  To  this  is  added  a  general  solution  of  the 
homogeneous  equation. 


I  .  ry  \  (-1  —  £Z  .  KL  Q  —  KL 

<t>(Z)=Ce  +Ae  +Be 


The  values  of  A  and  B  are  determined  from  the  boundary  conditions.  The 
quantities  of  interest,  the  transmitted  flux,  J+(L),  and  the  reflected 
flux,  J_(0),  can  then  be  evaluated.  The  quantities  now  include  unscattered 
flux.  ) 


J+(L) 


=  e 


-EL 


!i  * 


C( a+bE )  + 


C(a-b£) 

2 


/_5_  _ 

\bK  a  / 


sinh  KL 


den 


(24a) 


C( a+bE ) 
den 


J  (0) 


C  ja-b  E- 


_-£L  ,  a+bE 

e  (a-b£)  +  — = — 


( _a_  bK\ 
\bK  a  / 


sinh  KL 


den 


(24b) 


den  =  cosh  KL  + 


*(&♦¥)  ■* 


inh  KL 


(24c) 


An  interesting  check  on  this  solution  can  be  made  by  letting 
the  deposition  occur  in  a  thin  layer  (£+°°,  £s/£  finite,  with  a,  b,  and 
K  unchanged).  The  resulting  limits, 


£  /£ 

VL>  " 

J-(0)  If  j3  +  *  (ek  -  ¥)  sinh  KL/(den) 
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can  be  interpreted  (using  Equations  (21a)  and  (21b)  as  a  surface  source 
of  strength  Es/E,  half  entering  the  medium  and  being  transmitted  and 
reflected  according  to  the  prescription  of  Section  A-2 . 1 ,  and  half 
being  emitted  outward  immediately. 

A-2. 1.3  External  Point  Source,  Surface  Deposition.  A  slab  of  thickness 
L  is  illuminated  by  an  isotropically  emitting  point  source  at  a  distance 
h  below  the  slab.  The  relevant  equations  are: 


1  _3_ 
r  3r 


(r  mrp.} 


3  «Kr.Z) 


3Z 


K2<J>(r  ,Z)  =  0 


(25a) 


inside  the  slab,  with  boundary  conditions 


V'°>  -  ft 

h 

(25b) 

,  2  ,,2.3/2 
(r  +h  )  ' 

J_( r , L)  =  0 

(25c) 

f  J+(r,0) 

Jo 

2irrdr  =  Q 

Thus  Q  is  the  total  flux  into  the  slab.  This  is  just  half  of  the  total 
flux  emitted  by  the  source. 

The  solution  uses  standard  methods,  although  not  as  elementary 

as  the  method  of  Sections  A-2 . 1 . 1  and  A-2. 1.2.  As  a  first  step,  the 

problem  is  reduced  to  a  one-dimensional  problem  by  introducing  the  Hankel 
(33 

transform. 


/•  00 

$(P,Z)  =  /  <Kr,Z)  JQ(pr )rdr 


(26a) 


/•  00 

4>(r,Z)  =  /  4>(p,Z)  JQ(pr)pdp 


(26b) 


Equation  (25a)  is  multiplied  by  rJQ(pr)  and  integrated.  The  term 
involving  radial  derivatives  is  integrated  by  parts  twice.  Assuming 
r2<J >  and  ( 34>/3r)r  being  zero  at  r  =  0  and  r3/2<J>  and  r3^2 3<t> / 9 r  being  zero 
at  r  =  <*>  then  yields  the  relation 
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(27a) 


with 


$ — _  k2  <f>(p,Z)  =  0 
dZ^ 


k2  = 


P2  + 


The  corresponding  boundary  conditions  are,  using  Equations  13.51(5)  and 
3.7(13)  of  Reference  4, 


J+(P,0) 

JI(P,L) 

VP.Z) 


_Q_  p-Ph 

2  IT 

0 


a  <j>(P  >z)  +  b 


3<Kp,Z) 

9Z 


(27b) 

(27c) 

(27d) 


The  solution  of  these  equations  follows  the  pattern  of  Sections  A-2.1.1 
and  A-2.1.2. 


J+<P’Z)  =  |r  e_ph 


cosh  k(L-Z)  +  4  I 

a  kb ' 

kkb  a  ) 

)  sinh  k(L-Z) 

cosh  kL  +  4  | 

[a/kb  + 

(  a  j 

|  sinh  kL 

(28a) 


J  (P,Z) 


„  w  4  (t%  -  — )  sinh  k(L-Z) 

Q.  e-ph  _ s  \kb  a  / _ ' _ _ _ 

2lT  cosh  kL  +  4  sinh  kL 


(28b) 


These  relations  then  give  the  Hankel  transforms  of  the  transmitted  and 
reflected  fluxes: 


I+(P,L)  = 


■  A  e'ph/[ 


cosh  kL  +  4 


(JL  +  kb\ 
\kb  a  / 


sinh  kL 


(29a) 


J  (p,0) 


Q  p-Ph  i  (±.  _ 

2tt  s  Ikb  a  / 

cosh  kL  +  4  (4  + 


sinh  kL 


sinh  kL 


(29b) 
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The  actual  transmitted  and  reflected  fluxes  are  evaluated  as  integrals 
using  Equation  (26b): 


y*  oo 

j  (r  LI  -  ^  / 

e-phjo(pr)  P 

+  1,  ’  '  2tt  /  .  ,  T 

J  o  cosh  kL 

*  *  (ft  *  f) 

ap 

sinh  kL 

J-(r’0)  =  ft  / 

J  o 


e-p"J0(pr)p  i  (ft  -  Tj  31"h  kL 


o  cosh  kL  +  J  sinh  kL 


k2  =  p2  +  K2;  K2  =  1  /D 

ct 


A  number  of  quantities  can  be  derived  from  these  accumulated 
results.  The  total  transmitted  and  reflected  fluxes  are  easily 
obtained. 


i 

r°° 

Q 

J 

o  cosh  KL  +  £  | 

(ft  +  ¥ )  Sinh  KL 

/oo 

J_(r,0)  2irrdr  =  2ir  Jq(0,0) 


Q  Mie  -  -t)  sinh  KL 

cosh  KL  +  i  sinh  KI 


These  are  just  the  results  (adjusted  for  flux  Q  incident)  of  Equations 
(21a)  and  (21b).  A  quantity  of  some  interest  is  the  flux  from  a  strip 
of  width  xQ,  with  one  edge  going  through  r  =  0.  Let  Jg+  be  the  result 
for  the  transmitted  flux. 


s+ 


=  /°dx/  dy  J+  ^x2+y2  ,L^ 


(32a) 


_  Q 


r  dP  ( - — ) 

Jo  \cosh  kL  +  £  /  j-g  +  — 1  sinh  kL  / 


(32b) 


fX°  dx  f  dy  Jo(pV^x2+y2)  . 

J  o  •/  o 


By  Equations  13.47  (5)  and  3.4  (6)  of  Reference  4,  the  y-integral  is 


f 

J  o 


dy  JQ(p  V^2+y2)  =  -  J_j(px) 


=  ^  cos  px 

T  2p  T  irpx 


=  cos  px/p 


The  x-integral  can  be  done  immediately. 


J 


s+ 


q  r  sin  pxo 

11  Jo  p 


-ph 


cosh  kL  +  £ 


/X TW 

|kb  a  / 


sinh  kL 


dp 


(33a) 


The  analogous  result  for  the  reflected  flux  is  Jg 


A-1A 


J 

s- 


Q  f*Sin  PXo 

«'-ph  4  (4 

kb ' 
a  i 

)  sinh  kL 

TT  JQ  P 

cosh  kL  +  4 

rrr 

ikb 

sinh  kL 

dp 


(33b) 


The  limiting  case  as  xq-*“>  can  again  be  compared  with  previous  results. 
The  quantity  sinpxQ/p  behaves  as  a  Dirac  delta-function  (see  Reference 
5,  p.  29,  Equation  (47)). 


lira  sin  p  x 

X  -+oo  - — 

o  p 


it  S(p) 


lim  j  Q  _ 1 _ 

Xo"  S+  2  cosh  KL  +  J  ^  sinh  KL 


4 

^  a 
iKb 

Kb^ 
a  1 

|  sinh 

KL 

cosh  KL 

+  4  | 

a 

Kb 

sinh  KL 

These  are  just  half  the  results  in  (31a)  and  (31b)  as  they  should  be. 

A-2.1.4  Internal  Point  Source.  A  slab  of  thickness  L  is  illuminated 
by  an  isotropic  point  source  located  a  distance  h  (<L)  above  its  lower 
surface  (Z  =  0).  The  relevant  system  of  equations  to  be  solved  are 


II  (r  3<Kr,ZK  .  92<Hr  ,Z) 
r  3r  vr  3r  ’  az2 


K24>(r,Z)  +  6(r-Zh)/D=0 


(34a) 


where  Z  is  a  unit  vector  in  the  Z  direction,  and 

J+(r , 0)  =  J_( r , L)  =  0  .  (34b) 

The  delta-function  part  of  <J>  can  be  removed  by  introducing  the  infinite 

— ►  ■+ 

medium  solution,  G(r,r  ),  where  r  =  Z  h. 

o  o 


A- 15 


,  -K  r-r 
^  1  e  '  o' 

G<r’ro)  -  7fT;  T 


(35) 


It  is  easy  to  show,  by  differentiation,  that  G  satisfies  Equation  (34a) 
except  where  r=rQ.  To  evaluate  the  behavior  at  rQ,  the  usual  path  of 
integration  over  a  small  sphere,  centered  at  r  ,  is  taken.  Let  T  be  the 
sphere  volume,  S  its  surface,  p=  |r-ro|,  and  dft  an  element  of  solid  angle. 


/ 


(div  grad  G-K  G)  dx 


-/ 


grad  G*do  -  K 


/ 
J  T 


G  dt 


/&  ^  {f  -  $  ^ 

/  dft 

+o  4itD  J 


e~Kp  2 

P  dddp 


If  rQ  is  an  interior  point ,  /dfi  is  4n.  Hence  (div  grad  G  -  K2G  + 

1  6(r-r  ))  dx  =  -  g  +  =  0.  Using  these  results  the  solution  for  <j>  con¬ 

sists  of  two  parts,  G  and  ip ,  where  ip  satisfies  the  homogeneous  part  of  the 
equation . 


>(r,Z)  =  G(r,Z)  +  i|»(r,Z) 


(36a) 


G(r,Z)  = 


-K  \4-2+(Z-h)2 


4ttD 


\42+(Z-h)2 

F  -h  ('  Y  -  **♦(.,*) 


=  0 


(36b) 


(36c) 
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J+(r,Z)  =  a  i(j(r,Z)  +  b 


+  aG(r,Z)  +  b  3-(^Z] 


(36d) 


As  in  the  previous  section,  this  is  solved  using  Hankel  transforms. 


<KP,Z)  =  I  \p(r,Z)  JQ(pr )  r  dr 

J  o 


(37a) 


^  -  k2  iKp.Z)  =  o 


(37b) 


.2  ^2  2 
k  =  K  +  p 


(37c) 


The  transformed  boundary  conditions  are 


a  ip  C P  .L)  +  b 


=  -  aG(p,L)  -  b 


3G(p,Z) 


(38a) 


a  iKp.O)  -  b 


=  -  aG(p,0)  +  b 


(38b) 


G(p,Z) 


r  _l 

Jo  47TI 


-K  Vr2  +  (Z-h)2 


V2  +  (Z-h)2 


JQ(pr )  rdr  . 


This  integral  can  be  evaluated  using  Equation  13.47(2)  of  Reference  4 
in  a  manner  analogous  to  the  derivation  there  of  Equation  13.47(4): 
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G(p,Z)  = 


(39a) 


e~k | Z-h | 
4itD  k 


3G(p , Z )  _  e_klZ  hl  Z-h 

3Z  ~  4-nD  |  Z-h  |  ‘ 

The  solution  desired  then  leads  to  the  result,  using  b  =  D/2: 


(39b) 


J+(P,L)  = 


16tt  abZkZ 


_ (bk-a) _ 

cosh  kL  *  1  (4  + 


sinh  kL 


x  [(a-bk)  (a  sinh  kh  +  bk  cosh  kh) 

-kT 

+  (a+bk)e  (a  sinh  k(L-h)  +  bk  cosh  k(L-h))] 
+e-k(L_h)(a  +  bk)/8irbk 


(40a) 


J  (p,0)  = 


(bk-a) 


16tt  ab2k2  |cosh  kL  +  i  +  sinh  kL) 


x  [(a+bk)  (a  sinh  kh  +  bk  cosh  kh)  e 


-kL 


+  (a-bk)(a  sinh  k(L-h)  +  bk  cosh  k(L-h))]j 
— k  h 

+e  n(a  +  bk)/8iTbk  .  (40b) 

The  evaluation  of  the  actual  transmitted  fluxes  is  achieved  as  in 
Section  A-2.1.3.  If  desired  the  second  term  in  each  of  the  expressions 
can  be  explicitly  integrated.  Using  Equations  3.71(13)  and  13.47(2)  of 
Reference  4, 
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f‘  e-h 

J  o 


2+K2  J0(pr)pdp  e_K  |/h2+r2 


p2+K2 


'.2 .  2 
h  +r 


(41a) 


oo 

°-h 

O 


i  h  -K  f  h2+r2 

'p2+K2  J0(pr)pdp  =  ^ - 2—2 - 

h  +r 


(Kt7ts) 


(41b) 


However,  from  a  computational  viewpoint,  the  convergence  of  the  integral 
over  p  should  be  improved  by  not  separating  these  terms.  A  cursory 
examination  indicates  convergence  for  J+  to  be  as  e-k(L-h)  if  separation 
is  done,  as  e  if  no  separation  is  done. 

As  was  done  in  Section  A-2.1.3,  the  integrals  over  a  strip 
can  be  found. 


j  o 


sin  px 


-  J+(p,L)dp 


(42a) 


=  2/ 

«'  n 


00  sin  px 


2  J  (P.O)dp  . 


(42b) 


Similarly,  for  total  fluxes  out. 


yJ+(r,L)  2irr  dr  =  2ir  J+(0,L) 
o 


J_(r,0)  2nr  dr  =  2tt  J  (0,0) 


(43a) 


(43b) 


A-2.2  Time-Dependent  Solutions 

The  only  time-dependent  solution  that  is  practical  from  a 
computational  viewpoint  is  the  pure  diffusion  problem  of  a  plane  slab 
illuminated  by  a  uniform  external  source.  The  method  of  Laplace  trans¬ 
forms  is  used.  The  relevant  transforms  are: 
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f 

J  o 


-St 


<Kz,  s)  =  /  4>(z,t)e  dt 

To 


J+(z,s)  =  f  J+(z,t) 
—  J  o  — 


-st  ... 
e  dt 


(44a) 


(44b) 


F(s) 


f (t)e  stdt 


(44c) 


where  f(t)  is  the  incident  flux.  For  convolution  purposes,  the  case  of 
f(t)  =  6 (t ) ,  F(s)  =  1,  is  of  special  interest.  The  equations  and  boundary 
conditions  satisfied  by  the  transforms  are: 


32  =  (K2  +  ^)  $  (z,s) 

3  Z 


I  -  -  I  3bs\  7  ,  _  .  x  -  3<f>(z,s) 

Vz.s)  =  (a  “  TT)  *  (Z,S)  b  3^ 


JJL,s)  =  0;  J+(0,s)  =  F(s) 


Proceeding  in  a  straightforward  way. 


J+(L,s)  =  F(s)/G(X) 


J_(0,  s)  =  F(s)  i  '  IT)  siDh  XL/G(A) 

G( X )  =  cosh  XL  +  i  (£\~  ir)  sinh  XL 


where 


X2  =  K2  +  s/Dv  , 


3  bs 

a  =  a  "  2  ~ 


(45a) 

(45b) 

(45c) 

(46a) 

(46b) 

(46c) 
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The  relation  to  the  steady-state  solution  is  apparent,  since  the  integral 
of  the  time-dependent  solutions  over  time  is  just  the  value  of  the 
transform  for  s=0. 

The  actual  time-dependent  solution  is  obtained  by  the  usual 
inverse  Laplace  transform, 

i  f ioo+e  -  st 

J+(z>t)  =  27T  /  J+(z , s )es  as  .  (47) 

—  J  -i°°+e  — 

This  is  evaluated  by  contour  integration,  giving  a  sum  of  the  residues 
at  the  poles  of  J+(z,s).  These  are  at  the  zeroes  of  G(A).  The  exact 
evaluation  of  the  zeroes  is  rather  intractable,  involving  a  simultaneous 
evaluation  of  the  real  and  imaginary  parts  of  A.  For  simplicity,  an 
approximate  evaluation  was  used.  This  approximation  neglected  the  dif¬ 
ference  between  a  and  a.  In  terms  of  the  boundary  conditions,  this  is 

* 

equivalent  to  neglecting  the  3()>/3t  term,  compared  to  41,  in  the  currents, 

J  .  Using  this  approximation,  the  values  of  A  are  pure  imaginary.  For 

^  -th 
the  j  root 

A  =  ip.  (48a) 

J  J 

tan  njL  =  2abnj(b2nj2-a2)-1,  n?0 , j=l ,2 . . .  (48b) 

s  =  Dv(-n .2-K2)  .  (48c) 

J 

For  the  special  case,  f(t)  =  6(t),  the  inverse  Laplace  transforms  then 
yield : 

* 

Fluxes,  for  consistency. 


(49a) 


®  e‘  j  Dv  n 

J*(L-t)  ‘  &  -i  cdnjj 


t  /  bn  .  > 

-  e_sJ  I  fe.  *  -r1)  ‘“"V'  Dv"j 

-  £  - J-i  G'dij) -  (49b> 


-iG'(inj) 


L  sin  n  -L 
J 


cos 


njL 


(49c) 


+ 


V 


The  error  in  the  roots  is  of  the  order  of  D/L  times  the  value 
of  the  root.  Thus  values  of  the  root  near  zero,  which  dominate  the  late¬ 
time  behavior,  are  reasonably  well  determined.  The  early-time  behavior  is 
less  well  determined.  This  does  not  seem  a  serious  deficiency,  since  the 
diffusion  model  is  poor  for  early  times  in  any  case,  e.g.,  giving 
instantaneous  transmission  independent  of  thickness. 

A-3  STEADY-STATE  DIFFUSION  FOR  CONVEX  BODIES 

In  order  to  evaluate  the  reflection  and  transmission  by  more 
general  shapes,  a  numerical  method  has  been  developed  which  is  applicable 
to  general  convex  shapes.  The  restriction  to  convex  bodies  is  to  prevent 
one  part  of  the  body  serving  as  a  source  for  another  part  by  an  exterior 
path.  The  equations  to  be  solved  are 

V2  <|>(r)  -  K2<)> (r )  =  0  (50a) 

J  (r)  =  a  <(>(?)  +  b  grad  <J>( r)  •  n(r) 


(50b) 


where,  on  the  surface,  J+(r)  is  the  outward  flux  and  J_(r)  is  the  inward 

flux,  with  fi(r)  the  outward  unit  vector  normal  to  the  surface. 

The  usual  Green's  function  solution  of  this  type  of  equation 

involves  finding  a  function,  G(r,r'),  which  satisfies  the  Equation  (50a) 
■> 

except  at  r  =  r',  where  it  has  a  delta  function  type  singularity.  The 
usual  development  proceeds  as  follows.  For  a  volume,  T,  with  surface  E , 

[G(r,r')(V2f(r)-K2*(r))  -  <f>  ( r )  ( V2G(  r  ,  r  ’  )-K2G(  r ,  r  ’ )  ]  dx 

r,r')grad  ip( r)*n  -  4>(r)  grad  G(r,r')*n]  da  .  (51) 


Since  <j>(r)  satisfies  Equation  (50a)  throughout  the  volume  and  G(r,r') 

satisfies  it,  except  at  r  =  r',  the  volume  integral  can  be  shrunk  to  be 

over  a  sphere  of  arbitrarily  small  radius  about  r',  volume  T  ,  surface  Z  . 

s  s 

<j>(r )  (V2G(r  ,r ' )  -  K2G(r,r’))  dT 


f  «Kr)  (V2G(r  ,r  ’ )  -  K2G(r,r'))  dT 

•/  rn 


/  [div(c()(r)  gradG(r,r'))  -  grad  <j>  (r)  gradG(r,r') 

J  T 


K2G(r ,r' )]  dx 


-/ 


<J> ( r )  grad  G(r,r')  •  n(r)  do 


[grad  <(i  (r)  »gradG(r  ,r  ' ) 


+  K2G(r ,r ’ )]  dx  . 
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2  2 

Since  V  G-K  G  has  an  integrable  singularity  at  r=r',  the  second  integral 
vanishes  as  the  sphere  radius  goes  to  zero.  (Typically  G  is  like 
| r— r ' | —  ,  grad  G  like  | r— r ' | ~  and  dx  like  | r-r ' |  d | r-r ' | . )  In  the  sur¬ 
face  integral,  since  <K?)  is  continuous,  it  may  be  adequately  approximated 
by  $(?')  and  removed  from  the  integral.  Equation  (51)  can  then  be  written 
as 


-<Kr '  ) 


“►  A 

grad  G(r , r ' )• n ( r ) 


da 


[G(r,r'  )  gradi|)(r )  *n(r)-ij)(r) 


gradG(r , r 1 ) *n(r ) ]  da 


(52) 


The  limit-integral  multiplying  <t>(r')  is  here  left  explicit  since  its  value 
depends  critically  on  whether  the  point  r'  is  inside  the  surface  or  on 
the  surface.  In  the  former  case,  the  sphere  radius  will  be  small  enough, 
in  the  limiting  process,  that  the  sphere  surrounds  the  point,  r' .  If  r' 
is  on  the  surface  this  will  not  happen. 

In  the  usual  textbook  formulation,  G(r,r')  is  chosen  to  match 
suitable  boundary  conditions,  e.g.,  if  n*grad  <j>  is  given  on  the  surface, 
n*grad  G  is  zero  on  the  surface;  if  <)>  is  given,  G  is  set  to  zero.  In  the 
method  proposed  here,  following  Banaugh,^^  the  function  G  is  taken  to  be 
the  infinite  medium  function.  The  boundary  conditions  specify  a  linear 
relation  between  <t>  and  n*grad  <p  on  the  boundary.  Then,  if  r'  is  taken 
to  be  on  the  surface,  Equation  (52)  is  an  integral  equation  for  the 
unknown  values  of  <j>  or  n*grad  4>,  as  desired,  on  the  surface. 

For  the  problem  under  consideration,  where  p =  r-r'  and  dm  is 
element  of  solid  angle,  choose 


G(r , r ' )  =  e-Klr-r’l|?_?.|-l 


(53) 


/ 


lim  /  grad  G(r ,r ' )*n(r)do 

'l 


lim  /  (e  p  ^)p2  dw 

'y  o  P 
^s 


Jf 

L 


=  lim  /  e-Kp  (-kp-1-p-2)p2  dw 


k 


dw 
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Equation  (52)  then  becomes 


ft<t>(r')  =  /  [G(r,r' )grad<J>(r)*n(r)-(Kr)gradG(r ,r' )»n(r)]  da 


where  ft  is  the  solid  angle  subtended  by  the  surface  at  the  point  r', 
e.g.,  4tt  if  r'  is  inside  the  surface,  2ir  if  r'  is  on  the  surface  at  a 
point  not  a  "corner". 

For  the  problem  at  hand,  only  surface  values  are  of  interest 
and  it  will  be  assumed  the  surface  is  "smooth."  The  quantity  of  actual 
interest  is  J+(r)  with  J_(r)  given.  On  the  surface 

<f>(r)  =  (J+(r)  +  J_(r ) ) 

gradc|>(r)*n  =  ^  (J_(r)  -  J+(r)) 


^  (J+(r')  +  J_(r ' ) )  =  -  (J+(r)  +  J_(r))  grad  G(r,r')*n 

-  2E  "  J+(?))  G(?,?')j  da 


J+(?'  )  =  -  J_(r '  )  -  J ^  J_(r)  grad  G(r ,r '  )*n 

"  2fb  da 

~  fz  J+^>  6rad  G(r,r')*n 


+  G(r,r')  da  . 


The  first  two  terms  are,  in  principle,  known,  although  the  actual 
evaluation  of  the  integral  may  present  some  problems.  Thus  Equation 
(55)  is  an  integral  equation  for  the  values  of  the  emergent  flux,  J+ , 
on  the  surface. 
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Computationally,  the  procedure  for  solving  the  equation  con¬ 
sists  of  replacing  the  integral  by  a  set  of  algebraic  equations.  A  grid 
of  points,  r j ,  is  established,  each  having  an  associated  region,  £  ^  ,  and 
a  local  vector,  p j . 

V?j)  =  -  h  f.  J-(W  [grad  p.  G(?k+pk-?j)-« 

K 

-  S  G(VvV]  dok 


-  h  J  VW  [grad  pkG(?k+pk-?j)M 


+  b  G( rk+ 


W] 


da. 


-  h  £  J-(W  [« 


J_(?j+P3)  [grad  p  + 

J  J 


2tt 


£  J+<W  [eradP,G(?J^j'?j>-S  *  I  0<W*J)]  ddj 


ri . 
j 


(56) 


(gradp  means  a  gradient  with  respect  to  Pj  coordinates). 

J  The  first  term  is  known.  In  the  sums,  k^ j ,  the  integrals  may 
be  replaced  by  suitable  averages  times  the  area,  £k.  Only  the  last  two 
terms  require  special  care  because  of  the  singularity  of  the  function  G 
and  its  gradient.  A  demonstration  that  these  integrals  are  finite  and  an 
evaluation,  neglecting  variation  in  the  J's  over  the  area,  is  straight¬ 
forward. 
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/ 


G(r  +P  r.)  do  =  /  (e~Kpj  p.  2)  p  .dp .  d6 

Z  .  J  J  J  J  y  J  J  J 


-/ 


e  Kp  j  dp  .  de 
I  .  J 

J 


-  2ir  p . 


/z 

"  7T 


-  2tt  F  (where  irp.  *  =  I  . ) 

J  J 


“  2  VttJ 


/ 


gradp  G(r  +p 

Z  .  J  J  J 

J 


d» '  fz  -h  (^) 


^  pdpd0 


/. 


-Kp  /-K  1  \  p.n  .  « - 

e  (T  '  ~  pdpde 


-/ 


e-Kp(-K  -  I)  £1“  dpd6 


■+  /v 

£1“  ~  JP_ 

p  2R 


where  R  is  the  local  radius  of  curvature. 


f  Srad  Pj  G(^j+Pj  )  *nda  «  -^e-Kp(Kp+l)  R_1  dpde=-|v^T:  . 


Using  these  evaluations  of  the  integrals  then  yields  a  set  of  simul 
taneous  linear  algebraic  equations  for  the  J  (?.). 


(57) 


(58) 
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An  alternate  approximate  scheme  can  be  developed  which  does 
not  involve  the  solution  of  simultaneous  equations,  thus  avoiding  time 
consuming  procedures  (in  the  case  of  many  points).  It  can  best  be 
expressed  by  reverting  to  Equation  (55).  Let  S(r')  denote  the  known  first 
two  terms  on  the  right  side  and  G+(r ,r' )  the  multiplier  of  J+  in  the 
integral . 

J+(r')  =  S(r')  -  J  J+(r)  G+(r,r’ )  do 
J+ ( r ' )  ^1  +  G+ ( r , r ' )  da  \  =  S(r')  - J [J+(r)-J+(r ' )]  G+(r,r')  da  . 

•  /  e  c\  \ 


This  form  shows  some  immediate  advantage  since  the  integral  on  the  right 
is  less  singular  at  r  -  r'  due  to  the  zero  of  J+(r)  -  J+(r').  The 
approximate  procedure  solves  this  in  two  successive  approximations. 


J,  +(r ’ )  =  S(r ')/ 


|l  +  G+ ( r , r ' )  da 


(60a) 


J2.+(r’>  =  Jl.+(r’>  - 


J*  Jj  +(r)-J1  +(r ' )  G+(r ,r ' )  da 


(60b) 


r '  )  do 


The  actual  numerical  implementation  of  this  involves  the  same  replacement 
of  the  integrals  by  sums,  but  each  step  is  now  explicit. 

A-4  DIFFUSION  THROUGH  MULTIPLE  SLABS 

The  problem  of  transmission  by  multiple  contiguous  plane  slabs, 
in  the  pure  diffusion  case,  can  be  derived  from  the  single  slab  case. 

The  case  of  plane  source  and  point  source  follow  similar  development. 

The  procedure  consists  in  using  the  emerging  flux  from  the  top  (bottom) 
of  a  layer  as  the  entering  flux  on  the  bottom  (top)  of  the  next  higher 
(lower)  layer  and  deriving  a  set  of  relations  between  the  surfaces. 
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Planar  Diffusion  Through  Multiple  Slabs 

Assume  N  contiguous  layers  with  bottoms  at  d-^dg  <# 
thickness  ,  L ^  ...  L^.  In  each  layer, 


V2<fi(Z)  -  K24>.(Z)  =  0 


fi,±(Z) 


ai*i  +  bi 


3  4>  . 

1z 


The  boundary  conditions  are 


•  d. 


(61a) 


(61b) 


Jl,+(dl>  =  f>  JN,-(dN  +  V  =  0 

(62a) 

Ji+1,±<di+l>  =  Ji,±(di  +  V 

..  N-l  .  (62b) 

The  desired  quantities  are  JN  +(dn  +  LN)  and  J-^  (d^), 

In  the  ith  layer 

4>i(z)  =  AieKiz  +  Bie_Kiz 

Ji,+(z)  =  AieKiZ  (ai  +  biKi)  +  B.e"Kiz 

(a.  +  b.K. )  . 

l  —  l  i' 

In  matrix  notation, 

[Ji(di)]  =  [PjlCA.] 

(63a) 

[J.(d.+L.)]  =  [P.]  [Q.]  [A.] 

(63b) 

where 
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[Jj^x)]  = 


Ji,+(x) 


J  _(x) 


(63c) 


[P,]  = 


ai~biKi  ai+biKi 


a.+b.K.  a.-b.K. 

1  11  l  ii, 


(63d) 


[Qj  = 


eKiLi  0 


Lo 


-KiLiJ 


(63e ) 


[A,]  = 


K.d-n 

eii 


1_Bi  e"Kidi| 


(63f ) 


The  problem  has  been  formulated  in  terms  of  a  different  set  of  matrices: 


[T.]  =  [Ri ]  [A.] 


(64a) 


where 


[Ri]  = 


a.  +  biK1 


a.-b.K. 
1  11 


(ai-biK.)  eKiLi 


(aj.+biKi)  e  KiLi 


(64b) 


[t.]  = 


J.  (d.) 

i,~  i 


Lji)+(di+Li)J 


A  recursion  relation  is  easily  developed. 


[Tj  =  [Et]  [A.]  =  [Rj  [P±]  [Ji(di)] 
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By  the  boundary  conditions 


[J^d.)]  =  , 

[Ti]  =  [RiHPi]"1  [Q^]  [Ri_1] 

"  CMJ  tTi-l^ 

[tn]  =  [MN]  [mn_1]  ...  [M2]  [Tl] 

=  [>J]  [T\] 

[tn]  =  [Rn]  [Qn]  [Pn]  [JN(dN+LN>] 


JN,-<V 


LJN,+(dN+LN)J 


-  [«„]  C^]-1  CPN] 


-1 


JN,+  (dN+LN 
0 


[Tx]  =  [Rx]  [P^-1  [J1(d1)] 


■Jl,-(dl> 

f 

=  [Rx]  [P^-1 

-Jl,  +  <dl+V- 

Defining 


A- 31 


=  cosh  K 


iLi  -  2  (vr  *  -if) 


sinh  K.L. 
1  1 


1  t± i_  Vi\ 

4  -  2  \b.K.  -  a.  ) 


sinh 


ij]  -  Wd/Vi 


X.  ,  D.D.  ,-X.X.  , 

l-l  l  l-l  l  l-l 


[7H]  =  [Mn]  [M^]  ...  [M2]  , 


straightforward  matrix  operations  lead  to 


V*‘W  ■  ("i=i  si)  f/  I  [1'-V  [‘ij!  ■  (65a> 


A  similar  expression  can  be  derived  for  _(d^), 


Jl,-<dl>  = 


f  Dj  [1,-XN]  [7^]  “ 


[l,-xN]  [U/  ] 


(65b) 


Only  Equation  (65a)  has  been  introduced  into  a  computer  program. 

It  should  be  noted  that  the  requirement  of  contiguous 
layers  is  overly  restrictive  in  this  one-dimensional  case.  Actually 
the  requirement  is  for  non-overlapping  layers, 


d1<d1  +  Lx<d2<d2  +  L2  ...  <dN<dN  +  Ln  . 
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A-4.2  Point  Source  Diffusion  Through  Multiple  Slabs 

The  slab  geometry  is  the  same  as  in  Section  A-4.1.  It  is 
essential,  here,  that  the  slabs  be  contiguous  because  of  cylindrical 
divergence.  The  source  is  a  point  source  at  dQ(£d1).  In  the  ith  region, 
the  equations  to  be  satisfied  are 


1  _9_ 
r  3r 


0 


(66a) 


J. 


a* 


(66b) 


The  boundary  conditions  are  that  J+  and  J_  are  continuous  at  the 
boundaries  and 


Jl,+(rl’dl)  2i r 


dl'd0 


Cr2  ♦  (drdo)2]3/2 


JN,-(rl’VLN)  =  0  ‘ 

As  in  Section  A-2.1.3,  Hankel  transforms  are  used.  Then 
in  the  ith  slab 


92  4'i(P,z) 


-  <j>i (p  , z )  =  0 


(67a) 


ki 


2 


+  K. 

l 


(67b) 


—  3^  (P.z) 

Ji  ±  (P.z)  =  <t>i(p,z)  +b.  - ^ -  .  (67c) 
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J-i  .  (p,z)  and  Ji  _(p,z)  are  continuous  at  the  interfaces  and 
*7l,  +  (p’dl)  =  A  e_p(dl~dO)  =  f(p) 


(68d) 


jn,-(p’Vln)  =  0 


(68e) 


JL+1 ,  +  <p’di+l>  =  Ji)+  (P.VLi>-  i=1>2  •••  N-X 


(68f ) 


The  functions  of  interest  are  JN  +(p,dN+LN)  and  J-^  _(p,d1).  These 
equations  are  formally  the  same  as  Equations  (61)  and  (62)  in  Section 
A-4.1,  with  replaced  by  and  p  appearing  throughout  as  a  parameter. 


■  f<p> PJi 


(69) 


with 


I(P)  =  &  e“p(dl-do> 


[M^  = 


cosh 

k.L. 
l  i 

4  ( 

a . 

l 

b.k. 
l  l 

♦  p2 

[Ms] 

£“n-1 

'1 

-X. 

1 

.Xi 

/\  A 

D.D 

ih  -  *  (b-i-  *  -if) 


sinh  k.L. 

l  l 


sinh  k.L. 

l  l 
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The  actual  value  of  the  emergent  flux  must  be  determined  by  an  inversion 
of  the  Hankel  transform  as  in  Section  A-2.1.3. 


v+^-w 


ym  oo 

+(P-dN+LN)  Jo(pr)  pdp 

o’ 


(70) 


A-4.3  Diffusion  Through  Overlapping  or  Non-Contiguous  Slabs 

The  restriction  to  contiguous  layers  of  Section  A-4.2  and  to 
non-overlapping  layers  of  Sections  A-4.1  and  A-4.2  can  be  removed  in 
a  crude  approximation.  The  recipe  proposed  consists  in  fixing  the  tops 
of  all  14  cloud  types  and  the  bottoms  of  the  seven  low-altitude  cloud 
types,  adjusting  the  bottoms  of  the  remaining  layers  (i.e.,  the  mid- 
and  high-altitude  cloud  types)  to  give  contiguous  layers  and  adjusting 
the  density  of  each  layer  to  give  the  original  optical  thickness. 

This  procedure  is  admittedly  crude  but  seems  the  simplest  one 
that  even  approximately  simulates  the  actual  problem. 
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APPENDIX  B 
ROSCOE-IR 

NATURAL  CLOUD  PROGRAM  LISTING 
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C  THIS  PRESCRIPTION  NORKS  ONLY  FOR  A  PARTICULAR  SET  OF  INPUTS  ROSCUER  169 
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CL J  DTR  =  CONVERSION  FACTOR  FROM  DEGREES  TO  RADIANS  SCLOUO 

CLJ  PTE  CONHUN  SCLOUD 
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APPENDIX  C 

LISTING  OF  UPDATE  CARDS  NECESSARY  FOR  SAMPLE  PROBLEMS 


•IDENT  T?R0Bl/4 
•/  FOR  PROGRAM  ROSCOER 
‘INSERT  ROSCOER. 13 

COMMON/  T?M4  /  IFLAC4 
DIMENSION  KRR(161)/  RRTR(161,4) 

‘DELETE  ROSCOER. 94 
KKK  =  0 
DO  9  J=l,4 
DO  9  1=1,151 
RRTR(I,J)  =  0.0 
9  CONTINUE 
7  CONTINUE 
KKK  =  KKK‘1 
IFLAC4  =  -2 
•DELETE  ROSCOER. 176 

PIDX2  =  PI‘DXX/2.0  $  JR  =  161 
DO  301  1=1, JR 
‘DELETE  ROSCOER. 182 
•INSERT  ROSCOER. 188 
RRR(I)  =  X 
RRTR(1,KKK)  =  T4 

‘DELETE  ROSCOER. 192, ROSCOER. 194 

•DELETE  ROSC3ER. 206, ROSCOER. 208 

•INSERT  ROSCOER. 211 
NRITE(6,103) 

103  FORMAT  (181, ‘  RANGE  *,»  SC  AC  Cl  ‘, 

$  *  SC-AC-CI  •,*  INT(  SC  )  INT(  AC  )  INT(  Cl  )  •, 

$  •  INT(TOTAL)  •) 

RSM1  =0.0  $  RSM2  =  0.0  $  RSM3  =  0.0  $  RSH4  =0.0 

WRITE( 6,105)  RRR(1 ),{ RRTR( 1, J), J=l,4) ,RSM1, RSM2, RSH3, RSM4 
105  FORMAT  (1X,F8.2,4E12.4,4E14.6) 

DO  110  1=2, JR 

RSM1  =  RSM1  ♦  PIDX2*<  RRR( I)‘RRTR( 1,1 )  ♦  RRR( 1*1 ) *RRTR (1-1,1)  ) 

RSM2  =  RSM2  ♦  PIOX2 *(  RRR( I )‘RRTR( 1,2)  ♦  RRR( 1-1 ) *RRT«{ 1-1, 2 )  ) 

RSN3  =  RSN3  ♦  PIDX2*(  RRR(I)*RRTR(I,3)  ♦  RRR(I-1)‘RRTR(I-1, 3)  ) 

RSM4  =  RSM4  ♦  PIDX2*(  RRR( I )‘RRTR( 1,4)  ♦  RRk ( 1-1 )‘RRTR (1-1, 4)  ) 

HRITE(6,105)  RRR(I),(RRTR(I,J),J=1,4),RSM1, RSM2, RSM3, RSM4 

110  CONTINUE 

•/  FOR  SUBROUTINE  OIFPftN 
•INSERT  DIFPRM.81 

DATA  TF1234  /  -1  / 

•DELETE  DIFPRM. 105 

IF (  IF1234  .GT.  0  )  GO  TO  111 
•DELETE  DIF?RM.124 

IF(  IF1234  .GT.  0  )  GO  TO  27 
•DELETE  DIFPRM. 149 

IF(  IF1234  .GT.  0  )  GO  TO  29 
•INSERT  DIFPRM. 156 
IF1234  =  *2 

•/  FOR  SUBROUTINE  SCLOUD 
•INSERT  SCLOUD. 157 
E4  =  EMISS(l) 

I F (  ISUN.E3.0  )  T4  =  THANS(l) 

IF(  ISUN. E3.1  )  S4  =  TRANS(l) 

HSORCE(l)  =  HSAVE 


Figure  C-l.  Update  IDENT  for  Problems  1  to  4 


GO  TO  30 

*/  FOR  FUNCTION  TRANSF 
•BEFORE  TRANS. 9 

CUNM0.M/  TPM4  /  IFLAG4 
•INSERT  TRANS. 11 

OATA  KCLDZ  /  0  / 

•INSERT  TRANS. 101 

I F (  KCLOUD. EQ. KCLDZ  )  GO  TO  6 
KCLOZ  =  KCLUUD 

PRODKL  =  DKAPPA(KCLOUD)  •  CLDTHK(XCLQUD) 

ADIVKtf  =  OALF  A( KCLOUD )  /  ( DKAP PA ( KCLOUD) *D  BE  TA( K CLOUD) ) 

T  =  1.0  /  (COSH(PRODKL)  ♦  0.5*(  ADI  VKB  ♦  1.0/  AD  IV  KD)*  51  tall  (PRODKL)) 
WRITEI6, 13)  KCLOUD,  T 

13  FORMAT  (*0  TOTAL  TRANSMISSION  FOR  KCLOUD  =  *,14,»  IS  *,E12.4) 

6  CONTINUE 
•INSERT  TRANS. 385 

CLJ  ThE  MAIN  PROGRAM  INITIALIZES  1FLAG4  TO  A  NEG.  VALUE. 

IP(  If LAG4.GT. 0  )  GO  TO  1301 

CLJ  FUR  SAMPLE  PROBLEM  NO.  4  d E  WANT  TO  EVALUATE  THE  ANALYTIC 

CLJ  EXPRESSION  FOR  THE  TOTAL  ( SP ATI ALL Y- INTEGRATED )  TRANSMITTANCE 

CLJ  THRUUGH  THREE  LAYERS  ABUVE  A  POINT  SOURCE. 

CLJ  FIRST,  EVALUATE  THE  TRANSMITTANCE  S  AND  REFLECTANCE  R  FOR  EACH 

CLJ  Or  THE  THREE  LAYERS  SEPARATELY,  PER  EOS.  (31A)  AND  (31U)  IN 

CLJ  CLOUD  REPORT  APPENDIX  A. 

DKIT1  =  OKI *Tl 
DK2T2  =  DK2*T2 
DK3T3  =  DK3*T3 
HSKT1  =  SINH(DKITI) 

HSKT2  =  SINHIDK2T2) 

HSKT3  =  SINHIDK3T3) 


SI 

= 

1.0/(  COSH(DKITI) 

♦  0.5*11. / AL1  ♦ 

AL1 )*HSKT1 

) 

S2 

= 

1.0/(  COSHIDK2T2) 

♦  0 . 5*( 1. / AL2  ♦ 

AL2)*HSKT2 

) 

S3 

= 

1 . 0/ <  COSH( DK3T3) 

♦  0.5*( 1./AL3  ♦ 

AL3 )*  HSKT3 

) 

R1 

= 

Sl*0.5*(l./AL1  - 

AL1) *HSKT1 

R2 

3 

S2*0.5*( 1./AL2  - 

AL2)*HSKT2 

R3 

= 

S3*0 . 5* ( 1 . / AL3  - 

AL3) *HSKT3 

CLJ  NON  EVALUATE  THE  TOTAL  TRANSMITTANCE  FOR  THE  THREE-LAYER 

CLJ  CASE,  S123. 

S123  =  SI *S2*S3/(  (l.-Rl»R2)*(l.-R2*R3)  -  R1*R3*S2*S2  ) 

IFLAG4  =  *2 

WRITE! 6, 1300)  S1,S2,S3,R1,R2,R3,S123 

1300  FORMAT  {*0  TRANSMITTANCES  AND  REFLECTANCES  FOR  EACH  OF  THREE*, 

1  •  CLOUDS  ARE...*/*  T1,T2,I3  =*,3E14.5/*  R1,R2,R3  =*,3E14.5/ 

2  *0  TTANSMI TTANCE  FOR  THREE-LAYER  COMBINATION  IS...*, 

3  *  T123  =*,E1 4.5) 

1301  CONTINUE 
•DELETE  TRANS. 390 
•DELETE  TRANS. 413 
•DELETE  TRANS. 478 

TRANSf  =  fLUX  /  PI 
•DELETE  TRANS. 480 
•DELETE  TRANS. 493 


Figure  C-l.  (Continued) 
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•IDENT  TPROB/5 
*/  FOR  PROGRAM  ROSCOER 

•DELETE  ROSCOER. 176 
DO  301  1=1,1 
•DELETE  ROSCOER. 201 
•/  FOR  SUBROUTINE  CLD4T 

•INSERT  CLOUT. 7 

DIMENSION  LLNH( 160,3),ZCFLOS( 160, 4) 

•INSERT  CLDUT.l 23 

LLHH ( IDX, 1 )  =  LL  $  LLMH(IDX,2)  =  MM  $  LLMH(IDX,3)  =  HH 

ZCFLOS(IOX,II-1)  =  PIJKL»(1.0-CFL0S) 

•INSERT  CLDUT.l 39 

LLMH ( IDX, 1 )  =  LL  $  LLMH(IDX,2)  =  MM  $  LLMH(IDX,3)  =  HH 

ZCFLOS(IDX,II-l>  =  PIJKL»(1.0-CFLOS) 

•INSERT  CLOUT. 153 

LLMH( IDX, 1)  =  LL  $  LLMH(IDX,2)  =  MM  $  LLMH(IDX,3)  =  HH 

ZCFLOS(IDX,II-l)  =  PIJKL*(1.0-CFLOS) 

•INSERT  CLOUT. 173 

LLMH( IDX, 1 )  =  LL  $  LLMH( IDX, 2)  =  MM  $  LLMH(1DX,3)  =  HH 

ZC FLOS (IDX, II-l)  =  PIJXL*(1.0-CFLOS) 

•INSERT  CLOUT. 191 

LLMH( IDX, 1 )  =  LL  $  LLHH(IDX,2)  =  MM  $  LLMH(IDX,3)  =  HH 

ZCFLOS(IOX, il-l)  =  PIJKL*(1.0-CFLOS) 

•INSERT  CLDUT.209 

LLMH ( I  OX, 1 )  =  LL  $  LLMH(I0X,2)  =  MM  $  LLMH(IDX,3)  =  HH 

ZCFL0S(IDX,1I-1)  =  PIJEL*(1.0-CFLOS) 

•INSERT  CLDUT.228 

LLMH(IOX,l)  =  LL  $  LLHH(IDX,2)  =  MM  S  LLMH(IDX,3)  =  HH 

ZCFLOSdDX,  11-1)  =  PIJXL*(l.0-CFLOS) 

•BEFORE  CLOUT. 233 
URITE(6,303) 

303  FORMAT  (2H1  ,*  IDX  LL  MM  HH  *,*  TRANS  EMISS  *, 

$  •  CC=  2  CC=3  CC=4  CC=5  UT 

$  *  HT(SUM)  *) 

RSUM  =0.0 
DO  310  IX=1, 160 
RSUM  =  RSUM  ♦  UT(IX) 

URITE(6,305)  IX,(LLMH(IX, I),I=1,3),TRANS(IX),EMISS(IX), 

$  (ZCFLOS(IX,I),I=l,4),MT(IX),RSUK 

305  FORMAT  (2X,4I4,2X,8£12.4) 

310  CONTINUE 

•/  FOR  SUBROUTINE  DIFPRM 

•DELETE  DIFPRM. 105 
•DELETE  DIFPRM. 124 
•DELETE  DIFPRM. 149 
•/  FOR  SUBROUTINE  SCLQUO 

•INSERT  SCLOUD. 254 

WRIT£(6,103)  IDX,  TR  ANSI  160 ), EMISS  (  1 60  ),  HT  ( 1 60) 

103  FORMAT  (2 HO  ,4X,»  IDX  =  *,I5,5X,*  TRANS(160)  =  *,E12.4,5X, 

S  *  EMISSC 160 )  =  *, El 2. 4,5X,*  UT  (160)  =  *,E12.4) 

•INSERT  SCLOUD. 260 

MRITE(6,1 11 )  SLGCLD,SLCCLR,SFR,TUI,T12,T2TP,RDIST,SDIST0, 

S  THETI,THETE, PHIE 

111  FORMAT  (»0  SLG:lD,SLGCLR,SFR,THI,T12,T2TP,RDIST,SDIST0,»/ 

S  *THETI,THETE,PHIE*,/(5X,6E12. 4) ) 

WRITE(6,l 13)  (TRANS(I ),UT  (I), 1=1,11) 

URITE(6,115)  (EMISS(I),WT1(I),I=1,II) 

113  FORMAT  (»0  TRANS, UT  =•/(  5X,  8E1 2.  4)  ) 

115  FORMAT  (*0  EMISS, UT1  =•/( 5X,8E1 2.4) ) 

•INSERT  SCLOUD. 266 

MRITE(6,1 17)  (WT(I),UT1(1),I=1,II) 

117  FORMAT  (»0  UT,UT1  = */(5X, 8E1 2. 4 ) ) 


Figure  C-2.  Update  IDENT  for  Problem  5. 


MOENT  TPROB/6 

•/  FOR  PROGRAM  ROSCOER 

•INSERT  ROSCOER. 13 

COPMON/XKOOOO/  XCC ! 64), YCCI64 ), ZCC! 64 ) 

•DELETE  ROSCOER. 17 

DATA  SOLLAT,  SOLLON  /  0.0,  -0.7B5398  / 

•DELETE  RQSC3ER. 139, ROSCOER. 140 

IF (  I.EO.l  )  MRITEI6, 123) 

123  FORMAT  !*1CL0UD  ALAM  TCOEF  EHISS  XCC  *, 

1  •  YCC  ZCC  •) 

WRITE( 6,1 25)  K CLOUD( K), ALAH, TCOEF, TEHISS, XCC ( I ), YCC( I ),ZCC! I) 
125  FORMAT  !1X, I5,F8.2,E12.4,E12.4,3E12.4) 

*/  FOR  SUBROUTINE  CLOUD2 

•DELETE  CLOUD2.72 

•DELETE  CLOUD2.84 

•DELETE  CLOUD2. 96 

•/  FOR  SUBROUTINE  CLOUD3 

•DELETE  CL0UD3. 185 

•DELETE  CLOUD3. 216,CLQUD3. 217 

•/  FOR  SUBROUTINE  DIFFUS 

•DELETE  DIFFUS. 143 

•DELETE  DIFFUS. 165 

•DELETE  DIFFUS. 222 

•/  FOR  SUBROUTINE  LOS 

•BEFORE  LOSE. 67 

COHMON/XXOOOO/  XCC! 64), YCCI64 ), ZCC! 64 ) 

•INSERT  LOSE. 118 

XCC!L)  =  0.0 
YCC!L)  *  0.0 
ZCC!L)  =  0.0 
•INSERT  LOSE. 153 

XCC! LL)  =  XX!IC,JC) 

YCCiLL)  =  YY!IC,JC) 

ZCC! LL )  =  ZZ! IC, JC) 

•INSERT  LOSE. 203 

XCC!LL)  =  XX!IC, JC) 

YCCILL)  =  YY! IC, JC) 

ZCC! LL)  =  ZZ! IC, JC) 

•INSERT  LOSE. 245 

XCC! LL)  =  XX! IC, JC) 

YCCILL)  =  YY! IC, JC) 

ZCCILL)  »  ZZ!IC,JC) 


Figure  C-3.  Update  IDENT  for  Problem  6. 
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APPENDIX  D 

TRANSMITTANCE  AND  REFLECTANCE  OF  TWO-  AND  THREE-LAYER  PLANAR  SLABS 


To  help  validate  the  code  it  is  valuable  to  compare  the  code's 
results  with  information  obtained  by  other  means.  One  type  of  calculation 
performed  by  the  code  is  the  computation  of  the  photon  flux  emerging 
from  one-  two-,  and  three-layer  slabs  at  a  given  lateral  displacement 
from  the  vertical  through  the  point  source  below  the  slabs. 

Sample  Problems  1  to  3  in  Section  4-3  are  for  one-layer  slabs 
and  Problem  4  is  for  a  three-layer  slab.  To  provide  a  check  on  the 
code's  operation  for  Problems  1  to  3 ,  we  used  Equation  (32a)  in 
Appendix  A  for  the  spatially-integrated  transmittance  of  a  planar  slab. 

By  performing  a  spatial  integration  of  the  code-computed  flux,  we  com¬ 
pared  the  code's  result  with  the  analytic  result  in  Equation  (32a). 

Problem  4,  with  its  three  layers,  is  more  complicated.  Hence,  we  need 
another  exact,  analytic  expression  for  the  transmittance  of  three  layers, 
given  the  transmittance  and  reflectance  of  each  of  the  three  separate 
layers.  We  derive  such  a  formula  here;  it  is  applied  in  Section  4-3.1. 

The  approach  is  to  derive  formulas  for  the  transmittance  and 
reflectance  of  two  layers  and  then  to  apply  the  two-layer  formulas  to 
the  three-layer  configuration.  The  derivation  is  greatly  facilitated 
by  drawing  an  appropriate  diagram,  such  as  that  in  Figure  D-l.  The 
total  transmitted  flux  is  seen  to  be  a  geometric  series: 

T12  =  T2T1  +  R2R1  +  ^R2R1^  +  (R2R1^  + 


Vi 

l-r2ri 


(1) 


Similarly,  the  reflectance  from  the  lower  layer  includes  a  geometric 
series : 


R12  =  ri  +  Ti  r2  t1  +  R2Ri  +  <R2R1)2  +  (R2R1)3  +  •••] 


T12r2 
R1  +  1-R„R 


2  1 


(2) 
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2  2  3  3 

T.  T,  R„R,  T.T,  R_R,  T  T  R  R 
2  121  2121  2121 


/  /  /  / 


T1R2R21 


T1R2R1 


\  \  \  \ 

2  2  2  32  24. 

T1R2  T1R2R1  T1R2R1  T1R2R1 


Figure  D-l.  Illustration  of  contributions  to 
transmitted  and  reflected  fluxes 
for  a  two-layer  configuration 


D-2 
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Note  that  Equation  (1)  for  the  transmittance  is  symmetric  in  the  indices 
1  and  2  whereas  Equation  (2)  for  the  reflectance  is  not  symmetric. 

If  one  makes  the  assumption  of  conservative  scattering  (i.e., 
no  absorption),  so  that 

Tx  +  R1  +  1 
T2  +  R2  "  1  » 

then  one  can  show  that  the  reflectance  is  symmetric,  i.e.. 


R12  R21 


R1  +  R2  ~  2R1R2 
■*‘-RlR2 


and  that 


T12  +  R12  +  1 


For  three  layers,  consider  Figure  D-2  where  we  denote  the  pair 
of  Layers  1  and  2  as  Layer  A  and  Layer  3  as  Layer  B.  The  symbol  RA+ 
denotes  the  net  reflectance  for  flux  incident  on  the  Layer-1  side  of 
Layer  A  and  f?A_  denotes  the  net  reflectance  for  flux  incident  on  the 
Layer-2  side  of  Layer  A. 

The  transmittance  for  flux  incident  on  the  Layer-1  side  is 

tab  =  tbta  Ci  +  rbra-  +  <rbra->2  + 


1 “RBRA- 


The  reflectance  for  flux  incident  on  the  Layer-1  side  is 

rab  =  ra+  +  ta2rb  Cl  +  RbRa_  +  (Rbra-)2  +  •••  1 

-R  +  ^ 

-  ra+  i-rbra_  • 

By  using  the  expressions 

T  _  T2T1 
A  l-R2Rl 

TB  =  T3 

RA.  s  ■  *1  - 


Figure  D-2 .  Illustration  of  contributions  to  transmitted 
and  reflected  fluxes  for  a  three-layer  con¬ 
figuration  expressed  in  terms  of  those  for 
a  two-layer  configuration. 
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A-  "  21  "  "2  T  i-R1R2 

RB  =  R3 


we  find 


T  T  T 

T  -  T  =  3  2  1 

123  "  aAB  1~R2Ri 


T22R1 
R2  +  1-R. R 


12 


T  =  _ 

1 123 


T3T2TI 


<1-»2B1>  d-iyij)  -  T2  R3»l 


which  is  symmetric  with  respect  to  an  interchange  of  the  1  and  3  and 
hence  is  valid  for  the  source  irradiating  either  Layer  1  or  Layer  3. 
The  reflectance  R^23  ^or  a  source  below  Layer  1  is 


T 

p  —  d  -  r,  .  A  KB 

R123  -  RAB  RA+  +  1-R  R 

B  A- 


T12r2 

'  ri  +  i-r2r1  + 


/VlYr 

V-Wi)  2 


T22r1 

R3  [R2  +  l-R^ 


Rj_  (l-R^)  +  T12R2 


i-r2ri 


(T2T1> 


1-R2Ri 


(1-R„R2)  (1-R„R  )  -  T_^R  R 


+ 


whereas  the  reflectance  ^^21  ^or  a  source  above  Layer  3  is 


tb2ra- 

p  —  R  —  R  4-  ** 

n321  ”  nBA  kB  l-RgRA_ 


t!|„  ,  T228l 

*3  \k2  l-R-jRg 


R3  + 


1  -  R„ 


R  *  ^ 

2  1-R^j 


One  can  show  that  for  conservative  scattering,  i.e., 


that 


TV  +  Ri  =  1  for  i  =  1,2,3 


R123  R321 


and 

T123  +  R123  =  1  • 

Finally,  in  Figure  D-3  we  illustrate  the  enhanced  effective 
transmission  of  a  single  layer  due  to  a  reflecting  surface  located  below 
the  source.  Let  T  be  the  transmittance  of  the  layer  without  the  reflecting 
surface.  With  the  reflecting  surface  present,  there  are  two  contributions 
to  the  transmittance: 


Tx  =  T  +  TR2  +  TR4  +  . . .  =  T  / (1-R2) 


T9  =  TR  +  TR  +  TR  +  ...  =  TR/Cl-R'6) 


Tx+T2  =  T(l+R)/(1-R  )  =  T/ ( 1-R) 
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i 


The  factor  of  increased  transmittance  due  to  the  reflecting  surface 


1 

1-R  ’ 


a  result  stated  in  Section  2-2. 
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